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Russian ‘Scientific School’ System Explained 


964D0133A Moscow POISK in Russian 
No 44 (338), 2& Oct 9S-4 Nov 95S p { 


[Arucle by Yuriy Tretyakov, academician, Russian 
Academy of Sciences, and Igor Melikhov, academician, 
Russian Federation Academy of Natural Sciences under 
the “Portrait of the Problem’ rubric: “We Are Not 
Saving That in Which We Take Pride’ | 


{[FBIS Translated Text} In a recent issue of our 
newspaper (POISK No. 41), we published the Rus- 
sian Federation governmental decree ‘Concerning 
Governmental Support of the Russian Federation's 
Leading Scientific Schools." The program approved 
in that decree calls for the distribution of special 
grants among the country’s leading scientists on a 
competitive basis. Obviously, the measure will not 
solve all the pressing problems that have developed. 
The authors of the article that we are publishing to- 
day have analyzed the paths of development of Rus- 
sia’s scientific schools and presented their own ver- 
sions of supporting them. 


Modern science may be represented as a system of two 
types of scientific teams, namely, basic research labora- 
tories and base laboratories. Basic research laboratories 
are geared toward short-term study of individual man- 
ifestations of natural and social processes and are the 
most important type of laboratory at the present time. 
Base laboratories, on the other hand, are tailored to mul- 
liyear systematic study of natural phenomena. Base lab- 
oratories are generally assembled into groups, with each 
using approaches, ideas, and methods developed by one 
of the scientific schools. 


A scientific school is a team of very highly trained 
researchers that works on solving a large-scale problem 
and that during the course of its work, trains scientific 
workers in the art of research. The problem is generally 
of so large a scale that it requires the efforts of severa! 
generations of scientists and is solved on the basis of 
original ideas and methodologies developed by the main 
school. Schools are generally organized in the form of 
an institute or department or in the form of individual 
base laboratories or informal associations. 


Scientific schools are an important factor in society's 
development. It is precisely schools that determine 
how necessary it is to obtain new information and 
the pace and quality with which information must be 
obtained. It 1s schools that make expert evaluations of 
the information obtained in basic research and base 
laboratories. Without the official stamp of scientific 
schools, new information would be denied entry into 
mankind's store of knowledge. By interacting with one 
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another, the schools compare information trom different 
branches of science. We owe our harmonious world 
ouulook to scientific schools, Were it not for the schools 
Our understanding of nature would be as tar trom 
the truth as a heap of stones is from the Ta) Mahal 
Furthermore, each school is 4 source of creative impulse 
and all that is generally referred to as spirituality 


By our estimates, the world now has approximately 
§,000 natural science schools encompassing all dire 
lions of scientific exploration, 


In the nineties, Russia has had about SOO functioning 
natural science schools that have been responsible tor 
the country’s big scientific and technical achievements 
Many of these schools have by and large shaped their 
times. Natural science 1s now literally being transformed 
as a result of comprehension of the role of nonlinear 
processes and self-organizauon in nature. “Nonlinear 
thinking” has become a new paradigm tor researchers 
Indeed, the schools of L. Mandelshtam and A. Andronoy 
analyzed the role of nonlinearity both in general and in 
detail back in the forties to the extent that little can 
be added to their analysis today. And N. Kurnakov’'s 
school gave science a methodology of physicochemical 
analysis that has long dictated chemistry’s success tn 
the investigation of complex natural systems. 


Russian’s scientific schools have worked so effectively 
that they have created the impression that Russian sc! 
ence as a whole 1s in a state of well-being. In actual tact, 
Russian science began lagging behind Western science 
back in the seventies. In the eighties, the development 
of Western science accelerated sharply while Russian 
science continued advancing at its previous pace. From 
1980 through 1990, the Russian scientist's labor pro- 
ductivity increased by a factor of 4-4 while that of sci- 
entists in the United States increased by a factor of 20 
By 1990 our laboratories’ equipment and methods were 
four to six umes worse that those in the United States 
New ideas appeared with increasing rarity. Our labo- 
ratories’ cognitive abilities became significantly lower 
than those in the West. Russian science could not sup- 
ply the information support required for the pace of eco- 
nomic development required of leading countries. The 
financial support of science had to be increased. But 
the reverse occurred: In recent years investments in sci 

ence have been cut by two-thirds. As a result, updating 
of laboratories’ instrument park ceased, unul hy 1994 
the “average age” of most (60 percent) measuring 1n- 
struments was in excess of 15 years versus less than § 
years in the United States and Japan. Existing instru- 
ments began malfunctioning because of a lack of spare 
parts. Scientific associates’ wages fell catastrophically 

Some researchers left for abroad. Others entered differ 

ent spheres of activity. Researchers do not duplicate one 
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another's work in schools; each has his own inimitable 
task. For that reason, each researcher's departure is a 
blow to his school, And most sciences have lost more 
than a third of their researchers. The influx of youth 
into the schools has stopped, and directors have learned 
to hide their tears when talented individuals leave tor 
commerce 


The amount of available scientific literature has been 
reduced, which has inevitably resulted in a decrease in 
schools’ intellectual potential. The intensity of research 
has deceased by a factor of two to eight. The quality 
of scientific production has decreasea as well: Obsolete 
and worn-out equip,snent no longer provides the mea- 
surement precision that practice requires. Over the past 
5 years virtually all Russia’s schools have degraded at 
one pace or another. Some have even demonstrated bril 
lant achievements. Against the background of general 
collapse, however, such achievements were perceived 
as miracles 


Five years is an enormous period in the life of a 
scienufic school. In science, the measure of time ts the 
ume during which the volume of information doubles 
In the natural sciences the said period is between 5 
and 15 years long. In chemistry in particular it is 8.3 
years. This means that 75 percent of the knowledge 
on which chemistry will be based in 16.6 years is not 
yet known today. New knowledge will lead chemists to 
think in a new manner. Chemistry will be transformed 
into a science unknown to us today, and schools that 
have lagged behind scientific progress will not find a 
place in it. Although a 1- to 2-year delay in a school’s 
development is correctable in the future, a shutdown 
for 8 to 10 years is lethal: the school cannot be reborn 
regardless of the amount of financial investment. Those 
politicians at various levels who talk and write so 
much about science obviously do not agree with this 
statement. As a rule, they express guarded optimism and 
certainty that Russian scientists will be able to endure 
everything and quickly forget these trying times. In this 
case, however, optimism is unjustified. If the situation 
does not change radically within the next 2 or 3 years, 
some scientists may perhaps endure but there will be no 
schools. 


Without a system of scientific schools, Russian cannot 
be independent, much less a great power. A country’s 
independence and greatness is, along with everything 
else, determined by its information potential, which 1s in 
turn directly linked to the activity of its schools. There- 
fore, if the schools are destroyed, they will have to be 
created anew, which is a very complicated matter. Af- 
ter having analyzed the history of the development of 
the world’s largest schools, one can conclude that each 
school develops in three stages. In the first stage. the 
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potential head of the school develops hypotheses and 
methodological approaches to solving a problem. Next 
the research methods are assimilated and personnel are 
trained. In other words, the school is established. In the 
second stage, very simple model systems are studied and 
particular laws are established. In the third stage, com 
plex systems are studied and general laws of their be 
havior are deduced. In the said process, the information 
output from the school ts insignificant untl the middle 
of the second stage, after which it snowballs. A schoo! 
amasses most of its Knowledge in the third stage. The 
first stage lasts between 2 and S$ years, whereas the sec- 
ond stage lasts between 5 and 10 years. Consequently, 
during the course of that ume, a school produces virtu 
ally no information but rather simply absorbs resources 
that may later be recouped many times over. If, for one 
reason or another, a school collapses before the end of 
the second stage, the total yield from its work wall turn 
out to be low. The incubation period in a school’s de- 
velopment is an objective inevitability. The creation of 
a system of scientific schools therefore requires uncom 
pensated expenditures for many years. Obviously, those 
expenditures will be many times greater than what 1s 
needed to keep schools from collapsing. Only the ex- 
treme impoverishment of a country can justify the col- 
lapse of its scientific schools. Russia, which imports 
more expensive automobiles trom abroad than any other 
country throughout the world, can hardly be considered 
impoverished. In Russia, the development of science ts 
obviously considered a secondary task whose financing 
can be indefinitely suspended in this difficult period. But 
that view ts erroneous. Limiting the financing of science 
for 8 to 10 years is equivalent to killing it 


A rational amount of financing of scientific schools can 
be estmated by proceeding trom the expenditures on 
science in the United States. Russia does not have her 
Own experience in determining the reasonableness ot 
expenditures, whereas in the United States science ts fi 
nanced very economically. According to publi hed sta 
tistics, approximately $19 billion was spent on basic 
research in the United States in 1990. Fifty-four per- 
cent of that amount went to universities and research 
centers where scientific schools are concentrated. In the 
United States, approximately 400,000 individuals were 
involved in basic research, and the material and informa- 
lion support of their work (not Counting their wages) as 
eraged approximately $10,000 per individual researcher 
per year. In Russia, an average of five base laborato 
ries belongs to each school. If it is assumed that each 
of our 500 natural science laboratories has an average 
of 20 individuals, that means that approximately $0,000 
individuals are involved [in hasic research]. Given the 
same yearly funding to support the research of one in 
dividual (excluding his wages) as in the United States 
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this would mean the equivalent of $450 million each 
year. We are certain that that level of expenditure ts 
objecuvely necessary to provide the currently required 
intensity and quality of scientific invesugation. If that 
level is not reached, the effectiveness of the investiga- 
tion falls sharply and the work becomes uncompettive 
or even useless. The opumum level of expenditure does 
not depend on the will of some legislator or civil servant 
but is instead determined by the level of development 
of world science. 


Basic research must be financed by the government be- 
cause knowledge is not a commodity and the likelihood 
that the accumulation of knowledge be financed by mar- 
ket structures is low. In the United States, for example, 
firms contribute no more than 21 percent to the financing 
of basic science. Calling upon base laboratories to earn 
their own money by filling orders trom market struc- 
tures is therefore the equivalent of calling for a shift 
from basic research to applied research. Many scientific 
schools are now being forced to fill such orders, and 
the disastrous effects of the practice on the solution ot 
basic problems are already clear. 


Financing of scientific schools should be a priority. As 
a tule, schools have postgraduate students and doctoral 
candidates, organize conterences and seminars, create 
knowledge bases and expert systems, and publish mono- 
graphs and textbooks. One such example is G. Devy- 


atyy’s school, which not only developed the theory of 


deep purification and many methods of obtaining high- 
purity materials but which also trained 90 candidates 
and 15 doctors of sciences, published several mono- 
graphs and textbooks, conducted a university lecture 
course, created a data bank and exhibition/collection 
of high-purity materials, organized an all-union confer- 
ence, and worked on the publication of a special journal. 
All that required a great deal of funding over and above 
the amount provided by its research work. 


Because schools conduct lengthy scientific investiga- 
tions, their financing should be predetermined 5 to 8 
years in advance. A governmental program for devel- 
oping national schools that would reflect all directions 
of investigation for which Russia assumes responsibility 


before the world community may become the form of 


long-range financial planning. It 1s widely believed that 
a system of long-range financial planning that nearly 
excludes competition among researchers cannot provide 
the necessary pace of work and is therefore inferior to 
a system of short-term financial support on a compet- 
itive basis (a grant system). The fact that most scien- 
tific funds throughout the world use the grant system 
is confirmation of the grant system's superiority. It ts 
clear, however, that the nature of financing must corre- 
spond to the nature of the research being financed. If 5 
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to & years is required to solve a research problem, fi- 
nancing of the research work should be planned for the 
same ume period. The grant system is suitable tor ba- 
sic research laboratones and tor those base laboratories 
whose problems may be broken down into subproblems 
effecuve enough to win a competion. The system ts 
not suitable for scienufic schools conducting long-term, 
front-line scirenufic exploration. Ot Course, any school 
can aruficially single out an effective chain of tasks 
from its problems; however, this generally destroys the 
logic of their work. Both the pace and quality of long- 
term scientific invesugation should also be regulated by 
a scientometric analysis of laboratones’ activity. 


Such analysis makes it possible to identify scientific 
schools, determine which base laboratones are inca- 
pable of providing the required pace o1 work, and rank 
projects from the standpoint of the amount of financing. 
A general method of scientometric analysis has recently 
been developed. Once it has been adapted tor the respec- 
tive branches of science, it may be used to esumate the 
capabilities of all the country’s teams. Given the coun- 
try’s strained material resources, this type of “inventory” 
of laboratories seems necessary. After a base labora- 
tory has received long-term financial support, it must 
be subjected to monitoring by independent experts who 
evaluate the pace and quality of its work annually. An 
expert evaluation should be conducted by many experts 
so that it can be objective. It 1s therefore a good idea 
to combine expert evaluations conducted by individuals 
with expert evaluations conducted by groups, specifi- 
cally an evaluation that involves discussion of works at 
seminars and conferences. The said discussion may be 
conducted analogously to a dissertation defense, which 
is to say by using a secret ballot procedure. Seminars 
and conferences will then become stmult for accelerat- 
ing works with the most democratic form of monitoring 
of their course. We suggest that centralized financing of 
scientific schools within the framework of a governmen- 
tal program with democrauc monitoring of the course 
of the research based on a scientometric analysis of lab- 
oratories’ work will result in a pace of development of 
basic science consistent with Russia’s place among the 
world’s countries. 


Creation of S&T Lobby Urged 


Q64D0131A Moscow POISK in Russian 
No 45(339). 4-10 Nov 9S p ] 


[Article by Igor Goryunov: 
of Science?” 


Where Are You, the Party 


[FBIS Translated Text} Russian science. as a result of 
conducting market reforms, has become a demilitanzed, 
deideologized” and open world. So stated Minister 
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Boris Saltykov, speaking in the Palace of Scientists 
al a meeting dedicated to formulation of the strategy 
and tacucs of the scienufic community on the eve 


of the parliamentary elections. But in the opinion of 


Lyudmila Vartazarova, chairman of the Socialist Party 


of Workers, the ‘deideologization” has led to a loss of 


purposefulness,” that is, a loss of clear guideposts in 
the scientific and technical policy of the state. 


The threat which this carries was set forth in the 
Manifesto in Defense of Culture,” proclaimed at the 
meeting, which was signed by scientists and cultural 
workers of St. Petersburg. 


The St. Petersburg specialists feel that the course 
being pursued by the state has led to the development 


of exceptionally dangerous trends. The breakdown of 


Russian science 1s taking place, morals are going 
downhill and cultural values are being lost. 


The manifesto stated: “Specific measures for correcting 
the situation must be undertaken immediately.” In 
the announced document it is emphasized that the 
authorities of today do not comprehend what is going 
on. Repeated appeals of scientists to higher legislative 
and executive agencies and to the president “have had 
a zero effect.” 


However, for a true scienust a negative result also 1s 
a result. The authors of such messages were convinced 
that the appeals of scientists, even the most meritori- 
ous, to the pertinent authorities will scarcely find any 
response. In the Soviet years one could count on a re- 
action from the authoritative bodies. Now no heed is 
given to a cry from the soul even from the most emi- 
nent persons in the country. 


Accordingly, it is necessary to seek other methods for 
influencing organizations responsible for the fate of 
Russian science, education and culture. There must be a 
new concept of interrelationships between the scientific 
community and the authorities. Its principal provisions 
must be shared by scientists of the most different 
political persuasions: from the “Demvybor Rossi” to 
the Russian Federation Communist Party. 


However, the problem is not only the formulation of 


such a concept. This ts a necessary but by no means 
an adequate conditic.. It is necessary to organize a 
scientific lobby to exert an influence on the authorities 
and public opinion. The scientific community also must 
be concerned with its organization if it actually wishes 
to participate in the struggle f drawing up the budget 
and its expenditure, for an appropriate tax policy and 
science legislation 


Without a “brain center,” working out a= scientific 
ideology for the development of society, it will be 
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impossible for example, for an agranan or masked 
“criminal” lobby to prevail 


But for the ume being only some forms for conduct 
of the preelecuion struggle have been defined it has 
been decided to organize work on finding candidates 
in single-mandate districts who could, regardless of 
their party affiliation, represent the interests of science 
in the State Duma and Federation Council. After a 
corresponding analysis a list of such candidates must 
be published in the press. And it only remains for the 
sclenufic Community to vote tor “their people 


At the meeting no preference was given to any one of 
the political currents. The essence of the position of sci- 
enusts on this score was expressed by the adoption of a 
document calling on different political parties to assume 
specific obligations in ‘ormulating a state strategy in 
the scienufic and technical sphere. The RF Communist 
Party, “Zhenshchiny Rossu,’ Demokraticneskaya par- 
liya, Demvybor Rossi and Kongress russkikh ob- 
shchin” are even now ready to sign this document 


In general the background against which the discussion 
of the position of the scientific community transpired 
on the eve of the elections was created by an address of 
the American president, received by Russian screnusts 
via Internet and read from the pocium. Bill Clinton 
expresses disagreement with the congress, planning to 
reduce expenditures on science ‘for the sake of financial 
stabilization.” “The future does not have its voters. But 
we all at the same time are electors of the future 
And we must resist limitations on investments in those 
spheres which will be crucial for economic growth of 
the country in the 21st century. 


So it also is in Russia. Each of us makes his choice. 
But at the same time we are choosing the future of the 
country and its science 


Russian State Science Centers Form Lobbying 
Association 

964D0131B Moscow NEZAVISIMAYA GAZETA 
in Russian 11] Nov 9S p2 


[Article by Marina Kalashnikova: “A New Lobbyist 
Has Appeared in Russia. Russian Science For the Time 
Being Is Receiving Lip Service Rather Than Real 
Assistance’ } 


([FBIS Translated Text] Russian science 1s manifesting 
inventiveness in order to survive. On the eve of the 
elections it must push all possible buttons: state, pri- 
vate, public and emotional. An Association of State Sci- 
ence Centers (ASCC), whose introduction took place 
the other day at the All-Russian Aviation Materials Sci- 
enufic Research Institute (VIAM), has been organized 
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This is something like a Trade Union of Directors of 
State Science Centers (SSC) or a ‘volunteer tire com- 
pany tor rescuing children from a fire.” The associa- 
tion intends to lobby in the State Duma, the Federation 
Council, ministries and departments capable of supnly- 
ing money and affording advantages for saving the in- 
stitutes which are the leaders of Russian science. Can a 
structure tradiuonally appealing to state institutes be et- 
fecuve under market conditions and a change in national 
priorities? Wall it be capable of convincing the new 
holders of wealth and resources in Russia — the tuel- 
energy complex and the banks — to allocate a ‘chunk 
for the support of science? 


All of science cannot be saved if for no other reason 
than that Russia inherited it from the USSR in a greatly 
bloated condition. In 1993, by a presidential deciee, and 
then by a resoluuion of the RF government, the interim 
status of State Scienufic Center (SSC) was authorized 
for enterprises, scientific organizations and institutions 
of higher education. Those selected received the mght 
to special forms of state support for a period of 2 
years. The SSC designation was contended tor by 400 
scientific organizations among the 5,000 existing in 
Russia and 61 won this ttle. Regionally the leaders 
were Moscow and Moscow Oblast, where the SSC 
designauon was awarded to 40 institutions; there were 
12 in St. Petersburg, 3 in Novosibirsk, 2 1n Kaluga 
Oblast and | each in Ulyanovo, Orlovo and Tomsk 
Oblasts and Stavropol Kray. Eleven SSC are located 
in ¥ scientific cities — Zhukovskiy, Obninsk, Obolensk, 
Troitsk, Protvino, Mendeleyevo, Dmitrovgrad, Koltsovo 
and Nizhnyy Arkhyz. But for the time being there 
is no SSC in a single secret administrative-territorial 
formation 


For example, the Kurchatov Institute in Moscow and 
the All-Russian Scientific Research Institute of Grain- 
Legume and Groat Crops in Orlovo Oblast are included 
in the same list. The selection criteria did not pro- 
vide for proportional representation for either regions 
or branches. According to the members of the selecting 
commussion, the world rating of achievements was taken 
into account, but the principle of “who needs it more 
dominated. For example, space scientific research insti- 
tutes were not included because they were characterized 
by a relatively favorable situation and are enjoying fund- 
ing under other state programs. The conferring of the 
SSC designation does not change the forms of owner- 
ship. According to Nikolay Gorbatenko, responsible for 
work with the SSC in the Minnauka (Ministry of Sc- 
ence and Technical Policy), this is something like the 
‘boutonmere of the Legion of Honor, giving the nght 
on an interim hasis to enjoy the special tavor of the state 
The label is applied in front of the name of institution, 
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A) 


even if its name already contains such a word combina- 
lion. For example, the Vektor Virology and Biotechnol- 
ogy Center in Koltsovo village in Novosibirsk Oblast 
carnes a long ttle twice SSC. But what will one not 
bear tor the sake of one’s own salvation! 


However, for the me being such a granting of utles 
has brought no significant reliet. Intually the SSC idea 
assumed full maintenance by the state. A presidential 
decree and a government resolution only provided tor 
advantages. In actuality, however, great privileges were 
not received other than treeing trom taxation on land 
and real property. But this does not solve the principal 
problem because YO percent of the debts of the SSC — 
natural Monopolists are now primarily for energy, 
communicaucn services, water supply and sewerage 


The tdea of interbranch scientific and technical com- 
plexes, talling through in the late 1980's, again loomed 
behind the intention to award the SSC designation not 
only to one, but at the same time to several related inst- 
tuuons. At Zhukovskiy it was proposed that § organiza- 
tions be combined under this tle, but only two received 
it — the Flight Testing Institute and the Cential Aerohy- 
drodynamics Insutute — and each separately. The Rus- 
sian Academy of Sciences for a long time objected to 
the SSC idea, proposing its own versions. The opposi- 
tron was dropped when the advantage of gaming access 
to state investments became evident. The key question 

the mechanism by which the SSC relate to the state 

remains unclear to any and all of the “dramatis per- 
sonae the SSC directors, as well as the representa- 
tives of the Minnauka and the Department of Science 
and Higher Education of the RF government. While the 
Minnauka, Russian Academy of Sciences, State Com- 
mittee for the Detense Industry, State Committee for 
Property and other departments have sought a compro- 
mise, the SCC have had their water, electricity, heating 
and communication cut off and 6 of them today are on 
the brink of bankruptcy 


Among the immediate undertakings, German Zagaynov, 
general director of the association, announced that 
during the period 14-24 November there would be an 
exhibition-report of SSC activity during the last 2 years, 
to be held at VIAM facilities. This will determine their 
right to he voted into the list of SSC for the next period; 
protection of the interests of the SSC betore arbitration 
tnbunals and lobbying in legislative instituuions wall be 
discussed. Our political parties and movements do not 
always reflect the interests of narrow and specialized 
groups. A lobby may possibly help 


The awarding of the SSC status, in the opinion of 
Andrey Kulagin, head of a Minnauka administration, ‘18 
more likely a tnbute to past services than a policy for 
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determining future priorities.” There also is a directly 
opposite opimon. German Zagaynov, ASSC general 
director, asserts that the choice of SSC must lay the 
foundation tor a new state scienufic and technical 
policy. For the ume being it is being worked out, to 
use the expression of Mikhail Kirpichnikov, director 
of the Department of Science and Higher Education of 
the government, in a “Brownian movement” tashion 
Today i 1s possibly easier to survive on one’s own, 
betting on and playing up to some one of the natural 
monopolies or that narrow part of the military industrial 
complex which has already been able to find its way 
into the internauonal market (space and the exportation 
of competitive weapons). In the changing situation the 
greatest influence and the money for which science must 
contend is not in the bands of the defense branch, but 
with the oil and gas interests and the banks. Their need 
for scientific and technological development work 1s 
incomparably less than was true for the former detense 
industry and even the present-day defense sector and 
they differ in structure from the traditional activity 
of Soviet science. They will not die from starvation 
who are able to propose a cheap method for finding 
oil more rapidly and more simply expelling it from a 
well. But right now it is difficult to find sponsors tor 
a drug, invented by our science, which 1s capable ot 
replacing an anubiotic. The importation of drugs will 
bring incomparably greater profits to those who control 
the pharmaceuticals market. Will the ASSC or the state 
he able to straighten this out? 


In those matters within its competence, the ASSC has 
done plenty. Its initiative, in the persona of Mister Za- 
gaynov, is vigorously penetrating many fields and ot- 
fices. But there also are fears: how will the association 
distribute the money collected from one source or an- 
other tor support of the SSC and will the activity of the 
association not diminish the attention of the Minnauka 
and other departments which now fund the SSC” Yev- 
geniy Fedosov, director of the SSC Institute of Aviation 
Systems, feels that the principal problem is support ot 
the national experimental base. The Minnauka 1s in no 
position to assume this burden in its entirety. The State 
Committee for the Defense Industry for the most part 1s 
funding the development of specific technical systems 
Scientists are receiving support from grants. Will the 
ASSC get the necessary funds for developing the ex- 
perimental hase’ 
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At the NPO Aerodinamika there exists the bold idea 
of going tarther than the ASSC and establishing a 
scienufic exchange where the technical proposals ot 
the science centers and the orders of enterprises will 
find one another. A databank will be accessible on a 
commercial basis to both Russian and foreign capital, 
opening the way to the internauonal labor market 
An entire aircraft will not be ordered trom us, it 
is far more likely that the order will be for some 
part for at, the development of which will be more 
perfect and advantageous. assures Vladimir Sidenko, 
the general designer at Aerodinamika. This will increase 
the chances for small business, which will be more 
flexible and inexpensive under market conditions. It is 
easier (oO Close a smal! laboratory after it has performed 
its mission, Whereas the giants are regenerated into a 
self- maintaining system and a burden to the budget. The 
high economic etficacy of such a scheme possibly would 
make investment in science more attracuve for new 
Russian business and banks. And science also would 
turn out to be more profitable than when in the grip 
of the state. An exchange without fail must have tts 
own development bank, from which without additional 
bureaucratic steps there would be funding tor both basic 
and applied research. The exchange, bringing together 
the most diverse technical finds, will be able to become 
a place tor smoothing out the relations between the 
science lobby and the branch interests, as well as 4 tool 
of pohtcal influence 


For the ume being such projects sound utopian. But 1 
is essential to support both them and the ASSC if sci- 
ence is ever to break away from its long-time practice 
of begging the state to give’ and energetically involve 
itself into the “Brownian movement 
same equal basis as its other numerous present-day par- 
ticipants. The ASSC, according to Mister Kirpichnikovy, 
Is a great boon for the transition period, an attempt 
under new conditions to establish a new structure” In 
his opinion the Russian Academy of Sciences could not 
perform the function of a scientific lobby because ‘for 
too long a ume it has occupied an almost suicidal posi- 
tion and has greatly discredited itself. In the Minnauka 
it is hoped that on the eve of the elections the ASSC 
will convince thinking people to keep Russian science 
in mind. Mister Kirpichnikoy feels that ‘scientific spe- 
clalists must again feel that they have a sound root over 
their heads 


Of politics on the 
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Space Closed Gas Turbine Power Plants 
YS7AIIIOA Moscow VESTNIK MOSKOVSKOGO 
GOSUDARSITVENNOGO TEKHNICHESKOGO 
UNIVERSITETA: SERIYA MASHINOSTROYENIYE 
in Russian No 1, Jan-Feb 95 (manuscript received 
27 Aug 92) pp 25-32 


[Article by Vladimir Leonidovich Samsonov, candidate 
of technical sciences, head, Section for Space Closed 
Gas Turbine Power Plants, EM Scientific Research 
Insutute, Mcscow State Technical University iment N. 
E. Bauman, the first paragraph is an abstract; UDC 
629.7.54] 


{FBIS Translated Text] The results of multiyear research 
and development work on space closed gas turbine 
power plants carried out at the Moscow State Technical 
University iment N. E. Bauman under the direction of the 
author are presented. The fields of application in which 
such power plants are most advantageous are examined. 
The results of testing of experimental copies of the 
power plants, their principal assemblies and elements 
are given. Aspects of use of the developed power plants 
in conversion work are discussed. 


Space closed gas turbine power plants (CGTPP), which 
in the literature are often simply cailed power plants 
(PP), operating in a Brayton cycle, have a number ot 
advantages. In contrast to a closed piston-driven PP, 
operating in a Surling cycle, they are characterized by a 
higher specific power and a long continuous-operation 
useful life. In contrast to closed steam turbine PP, oper- 
ating in a Rankine cycle, CGTPP are characterized by a 
high efficiency and high reliability, which ts attributable 
to the use of a working medium in the form of a mix- 
ture of neutral gases having excellent compatibility with 
construction materials. The enumerated qualities result 
in a preference for space CGTPP among other PP of the 
generator class. 


Due to the use of gas bearings in CGTPP turbogener- 
ators there 1s virtually no wear in these critical units. 
These PP therefore have a long continuous-operation 
useful life attaining 30 or more years. As indicated by 
research, in the range 0.5...5,000 kW space CGTPP, not 
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yielding in reliability to PP based on solar cells, and 
also based on thermoelectric and thermoemission con 
verters, considerably surpass them with respect to such 
Indices as Compactness, Mass and economy 


Figure | shows the typical layout of a space CGTPP 
{1}. 
t 















































Fig. 1. Layout of space CGTPP: 1 — radiator; 2 — 
payload; 3 — heat source; 4 — recuperator. 


Key: (1) Gen; (2) HRHE 


The power conversion system [T includes 4° units 
turbogenerator-compressor unit: (TGC), recuperator, 
source heat exchangers (SHE) and heat removal heat 
exchanges (HRHE). These units are connected by lines 
thereby forming a closed gas circuitin which a helium- 
xenon mixture with a molecular weight opumum for 
the PP power level circulates. The principal unit of this 
system the TGC unit for increasing reliability 
has a single shaft with a common rotor resting on two 
supporting and one two-sided supporting gas-dynamical 
hearings 


In the heat source system it is possible to use a solar 
concentrator with a heat pickup, nuclear reactor of 
radioisotope chamber 
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Fiy. 2. Layout of stand for tests of CGTPP units: | — compressor unit; Il — turbine part; III — 


recuperator part; IS 


— heater part; V — dynamic simulator of TGC and electric subsystem of CGTPP; 1 


— vacuum valves; 2 — throttle; 3 — recuperator. 


Key sit Mip +) Water 4 KP/S, (5 HRHE: (6) I; 
T} SSIDa fh system operates 


either | e direct release of heat trom the heltum- 


xenon mixture through the wall into space or with an 
ntermedi at juid circuit, depending on the makeup 


t the CGTPP on the spacecraft. It 1s proposed that a 


droplet radiator with the direct radiation of heat into 
space by hot droplets of a liquid heat transfer agent 
irculating in the heat dissipation system [2] be used 


tor high-power CGTPP, the area of whose radiator 


surf ac ‘ SS y of their placement aboard 

‘f we er 

The elec system TV is used tor conversion, rectifica- 
r distribution of the electric cu cnt produced by 

tne general the TG Hoot the r conversion 

system, dependiag on the n iKeup Of the pacecralt pay- 


The gas regulation system V is intended tor compen- 


? 7. yy ‘ 7 los i . 7 en) | 5 > 7 . 
sating for muicroleakage from the closed circuit of the 


power conversion system, as well as for regulating the 
guantty of gas if when there are large power con- 
sumption changes 

Using such a system American specialists have con- 
structed a surtace CGTPP with a power N. = 2-15 kW 
with an electric simulator of the heat source (prototype 


Aspace Power pian in the course of the tests accu- 


mulating f $0,000 hours of continuous fault- 
free oper “ legradation of the output char- 
acteristics {4}. The efficiency of this plant exceeds 30 
percent with an electric power of 7 kW. As a result, 
Americ ecialists preferred a solar CGTPP over so- 


cells for new-generation orbital stations with a power 


(7) Gen 


consumption 20-300 kW. According to their data, such 
a choice will make it possible to save 3-4 billion dol- 
lars during operation of the Freedom orhital stauon over 
the course of 30 years [4]. We obtained similar results 
in a study of the possibility of using a solar CGTPP 
for new-generation Mir orbital stations and a general- 
purpose space platform. 


Due to high economy and reliability the use of space 
CGTPP with a power N_ = 0.5-10 kW with a radioiso- 
tope chamber as a heat source 1s extremely promising 
The high efficiency of this plant (20-35 percent) not only 
reduces expenditures on loading ?“Pu into the radioiso- 
tope chamber by tens of millions of rubles (in compars- 
ison with the case of use of traditional thermoelectric 
generators), but also ensures good layout properties of 
the plant in the spacecraft. Such plants have no alterna- 
tives for spacecraft with high power consumption during 
flights to Mars, to distant planets and to outer space 


The combination of a CGTPP with a nuclear reactor for 
spacecraft with a high power consumption (more than 
100 kW), as a result of the high efficiency of power 
conversion (up to SO percent), provides substantial 
advantages over thermal emission PP with respect to 
size and weight characteristics. A high efficiency of the 
CGTPP can be attained by use of ceramix materials, as 
well as materials based on carbon fibers, which ensures 
a temperature increase at the entrance into the turbine 
up to 2,400 K 


Space CGTPP, in combination with a turbine gas 
expander, are extremely promising tor producing a great 
guantity of cold at the cryogenic temperature level 
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necessary for the cooling of infrared transducers and 
containers for the storage of cryogenic fluids under 
orbital conditions. When using such solar plants it 1s 
possible to reduce the overall dimensions of the plant 
by 1-1.5 orders of magnitude in comparison with the 
dimensions of solar cells in combination with traditional 
piston- driven cooling plants. In this case the TGC unit 
in the CGTPP 1s replaced by a turbocompressor-gas 
expansion unit 


And finally, space CGTPP are without an alternative 
in a case when for the spacecraft payload there must be 
delivery not only of electric power, but also the pumping 
of a great quantity of gas through it. This applies, in 
particular, to such users as supersonic electroionization 
CO, and CO lasers, as well as a Cu vapor transonic 
laser, intended for power transmission in space [5]. 


It can therefore be concluded that CGTPP are some of 
the principal power sources for space vehicles used for 
different purposes. 


The Moscow State Technical University imeni N. E. 
Bauman, with the support of the aerospace industry, has 
for a long time been investgating and developing space 
CGTPP in a broad range of electric powers. However, 
the experimental study was made of the PP as a whole 
and their assemblies for an electric power level 1...20 
kW 


Due to the uniqueness of the gas-dynamical and heat ex- 
change processes in the assemblies and elements of the 
centrifugal compressors, centripetal turbines, gas bear- 
ings of the TGC block rotors, recuperator, heat ex- 
changer for heat dissipation, radiator, etc., transpiring 
at low Re numbers (< 10°) and numbers Pr = 0.2...0.5, 
characterisuc of helrum-xenon mixtures, stand equip- 
inent was developed for their autonomous investigation 
and testing. 


Figure 2 1s a diagram of this stand. It consists of 


compressor, recuperator, heating and turbine parts, as 
well as the power plant electric system intended for 
investigation and testing of the supporting systems of 
the TGC unit rotors. 


The basis for the compressor part of the stand 1s a single- 
shaft low- temperature turbocompressor which includes 
the centnfugal compressor K1 and the drive air turbine 
TI to be tested. The unit rotor rests on gas bearings. The 
working medium of the tested centnfugal compressor 
(He- Xe, Kr, Ar and others) is separated from the air by 
means of a cutoff device 


When testing other CGTPP elements the compressor 
part of the stand serves as the source of the working 
medium 
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In invesugations and tests of the recuperator and cen 
tripetal turbine (CT) all the components 
shown in Fig. 2 were in operation. The turbine part o 
the stand is the most complex. Its basis is a single- 
shaft high-temperature turbocompressor with the CT 
to be tested, braking compressor and rotor on high- 
temperature gas-dynamical bearings 


of the stand 
t 
: 


The stand (part V) 1s used in invesugat.ng and tesung 
the gas- dynamical support systems of the rotors for the 
CGTPP single-shaft TGC units. The tests are carned out 
using a dynamic simulator of the TGC uni with heating 
of its hot part by quartz lamps 


The hot assemblies and elements of assemblies of all 
parts of this stand are heat insulated. The electric heater 
of the heating part of the stand also 1s carefully heat 
insulated 


The stand was outfitted with thermocouples, pressure 
transducers, current eddy devices for measuring smal! 
movements in elements of the gas bearings and wheels 
of the CC and CT, flow meters, chromatographs, tre 
quency meters, etc. Highly economical stages of the 
CC with an efficiency 80 percent, CT with an efficiencs 
90 percent, high-speed a-c generators (n = 24,{%))..45 
rpm), supporting gas-dynamical systems of the TGC 
unit rotors in the electric power range 1...20 kW and 
compact recuperators with a degree of recuperation up 
to 96 percent were investigated and tested on this stand 


A surface prototype of a space CGTPP with a power N_ 
= 3 kW was constructed on the basis of the developed 
technology. Figure 3 1s a diagram of this stand. The 
role of a space heat source (solar concentrator heat 
pickup) 1s played by an electric heater, whereas the role 
of the radiator 1s simulated in the form of a liquid-water 
heat exchanger. The gas in the HRHE 1s cooled oy a 
silicone fluid which is then cooled by water from the 
main network. 


The specificauons of the surface copy of the space 
power plant are as follows: nominal electric power 
4.0 kW, current voltage 220 V, degree of increase 
IN gas pressure in Compressor 1.75, frequency 1,433 
Hz, rotation of TGC unit rotor 40,000 rpm, degree of 
recuperation 95 percent, gas temperature at entrance into 
turbine 1,130 K, at entrance into compressor 408 
K, working medium — helium-xenon mixture with a 
molecular weight 83.8. The TGC 1s designed in the form 
of a single-shaft unit whose rotor rests on thin-film gas 
hearings of the hinged type. The outer diameter of the 
CC wheel is 90 mm, that of the CT wheel is 105 mm 
and that of the generator rotor 1s $8 mm 


The body of the generator stator 1s cooled by a silicone 
fluid. The temperature of the hottest point of the stator 
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of this quadrupole generator with a massive rotor does 
not exceed 200°C. The wheels of the turbodynamo were 
cast, the wheel of the CC was made from tutanium alloy 
and the CT wheel was made from a heat-resistant alloy. 
The bimetallic generator rotor is fabricated by heat- 
diffusion welding. The working surface of the thin-film 
gas bearings has a polyamide coating. The mass of the 
TGC unit (Fig. 4) is 43 kg and the mass of its rotor is 
4.2 kg; the overall dimensions of the unit are: maximum 
diameter (along the connecting pipes) 330 mm, length 
580 mm. 


The recuperator is a countercurrent lamellated-ribbed 
heat exchange unit with a matrix formed by heat 
exchange surfaces with short straight offset ribs. The 
construction material is corrosion-resistant steel. The 
thickness of the ribbed surface is 0.1 mm and that of the 
separation plates is 0.2 mm. The overall dimensions ot 
the recuperator matrix are: length SOO, width 300 and 
height 400 mm. The recuperator weight is 80 kg. 


The HRHE is made in a stand variant. Its lamellated- 
nbbed construcuon ts fabricated from aluminum alloys. 
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Fig. 3. Diagram of surface prototype of space 
CGTPP with power N. = 3 kW. I — radiator 


simulator; [1 — ballast load; II] — electric 
regulation unit; [V — payload; 1 — electric heater; 
2 — recuperator. 


Key: (1) TGC; (2) MPa; (3) Gen; (4) Water; (5) HRHE; 
(6) kg/s 
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The electric heater has Nichrome tubular heat-releasing 
elements. It is constructed in the dimensions of the solar 
pickup-accumulator tor a CGTPP with a power N. = 
20 kW. The power plant is startled up by means of a 
generator operating in a motor mode. The synchronous 
frequency of rotor rotation in the TGC unit is 24,000 


rpm. 


The total ume which the power plant funcuoned ex- 
ceeded S00 hours with 100 “start-stop" cycles. The re- 
sults of the tests indicate a potential possibility of at- 
taining a prolonged (years) Contunuous-operation useful 
lite. The resulting plant efficiency attains 24 percent. 


Detailed information on a surface prototype of a space 
CGTPP with N. = 3 KW is given in [6]. 


Applicable to solar CGTPP of promising orbital stations 
and space platforms, the Moscow State Technical Uni- 
versity iment N. E. Bauman is investigating designs of 
a space solar pickup-accumulator (Fig. $). In particular, 
experiments are being carried out for determining the 
operability of niobium capsules with LiF the princi- 
pal elements of such pickup-accumulators. In order to 
determine voids accompanying LiF hardening, one of 
the capsules making up this element was first subjected 
to tomographic examination with the simultaneous pro- 
cessing of the results on a computer after thermal! cy- 
cling in vertical and horizontal positions. Intormation 
on these experiments was presented in [7]. 


Research on the conversion applications of space 
CGTPP indicated that they may play an important role 
in solving ecologic and resource-conserving problems 
on the Earth [8]. 


In particular, it seems promising to use low-power 
CGTPP operating on organic tuel for cathode protec- 
tion Of main gas pipelines, radio relay Communication 
systems and small farms in inaccessible regions of the 
planet instead of low-efficiency organic-fuel steam tur 
bine PP employed at the present time 
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Fig. 5. Tubular element with niobium capsules filled 
with LiF of space solar pickup-accumulator. 
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Fig. 6. Schematic diagram of CGTPP with power N. = 3 kW operating on natural gas: 1 — radiator; 2 - - 
TGC unit; 3 — electric regulation unit; 4 — gas heater; 5 — electric fan; 6 — final user of released heat; 7 
— air recuperator. 


Key: (1) kW; (2) V; (3) kg/s; (4) MPa; (5) W; (6) air; (7) g/s; (8) Generator; (9) Natural gas; (10) To atmosphere; 
(11) Recuperator; (12) HRHE; (13) To heat user; (14) From heat user 


Figure 6 is a schematic diagram of a CGTPP with a 
power N, = 3 kW operating on natural gas for line users 
On main gas pipelines. For making maximum heat use 
in the heat source system of this CGTPP there is an air 
recuperator and a terminal gas-fluid heat exchanger of 
a system for heating communication system equipment. 
Independent tests of this system with a total duration of 
47 hours with more than 100 “on-off' cycles indicated 
its reliable and efficient operation [9]. Tests are now 
also being carried out of experimental copies of CGTPP 
with a power N. = 3 kW operating on natural gas which 
include the proposed heat source system. 


The conclusion can therefore be drawn that the use of 
low-power CGTPP both in space and on the ground is 
promising. 
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Gas Turbine and Steam Power Plants With 
Ceramic Modules 
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[Article by A. A. Kulandin, doctor of technical sciences, 
professor, academician, Russian Academy of Trans- 
portation, V. Kh. Dubershteyn and A. V. Sudarev; the 
first paragraph is an abstract; UDC 621.438] 


[FBIS Translated Text] A study was made of the ap- 
plication of ceramic construction materials (CCM) for 
the high-temperature uncooled feed-through parts of gas 
turbine power plants (GTPP) in the case of their au- 
tonomous operation and in combination with a recov- 
ery unt’. The characteristics of a thermodynamic analy- 
sis are given with optimization of the unit-modular ap- 
proach to the finalization of design schemes for GTPP 
and the desirability of using a cycle with intermediate air 
cooling in the compression process is pointed out. The 
possibility of increasing the efficiency of thermal electric 
power plants by 15-18 percent is established. 


In the technical-economic assessment of thermal elec- 
tric power plants the most important factors taken into 
account are the thermal economy or total efficiency, de- 
termining the specific expenditures of fuel and the spe- 
cific thermal load on ecologic conditions in the location 
region, as well as the specific (per unit of produced 
electric power) discharge of harmful substances (NO,, 
CO and others). Among the different types of stationary 
thermal power equipment it is gas-vapor power plants 
with a recovery unit (GVPP/RU) which have the high- 
est indices 


The total efficiency (7, of GVPP/RU is dependent 
primarily on the efficiency of the gas turbine (GTPP) 
and steam turbine (STPP) plants (1), and n,), the degree 
of binanity of the thermodynamic cycle (4, and the 
recovery unit (7,. With an accuracy to 2-3 percent the 
efficiency of the GVPP/RU is defined by the equation 


m21-(1-nm) Cl - 9) (1) 


Naturally, it 1s desirable that all the components entering 
into ths equation be as high as possible. But they are 
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interrelated (an increase in some results in a decrease 
in others) and these relationships determine opumization 
of the parameters of the working process and the GVPP/ 
RU cycle with respect to ensunng a maximum total! 
efficiency. There is no need for optimizing solely the 
Maximum temperature cycle (T_ [cc = closed cycle], 
which must be as high as possible when using any 
GTPP cycles. As indicated by an analysis, an increase 
in T. for the GVPP/RU has a greater valuc than for an 
autonomous GTPP. 


} 
I 








2 4 


Fig. 1. Schematic diagram of GTPP with ceramic 
modules: 1) air entry into modules with ceramic 
turbines; 2) modules with ceramic turbines; 3) 
metal power turbine; 4) exhaust; 5) variant of 
common low-pressure compressor; 6) air entry into 
common low- pressure compressor. 


In the traditional embodiment of elements of the feed- 
through part of a metal turbine the T_ increase 1s limited 
by both the increase in expenditures of compressed air 
for covling and the uniqueness of the employed ma- 
terials and coatings. A complicauon of the technology 
results in a sharp increase in the cost of the items. The 
traditional materials for GTPP turbines are heat-resistant 
alloys, those close to the limit of reserves with respect 
to an increase in the maximum cycle temperature 


The use of ceramic construction materials (CCM) tor 
the fabrication of high-temperature uncooled elements 
of the feed-through parts of a GTPP makes possibie a 
substantial increase in efficiency of both an autonomous 
GTPP, but especially a GVPP/RU. With an increase 
in T. in this case there is assurance of a continuous 
increase in 7,, it becomes feasible to realize a GVPP/ 
RU with n, = 1 with a simultaneous increase in 7, and 
n,. Thus, all four coefficients entering into equation (1) 
on which the total efficiency 7, 1s dependent, increase 


However, taking into account exisung technologies. a 
nuriber of restrictions are imposed on items made of 
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CCM: the mean mass temperature T7 must not ex- 
ceed 1,300°C; the size of the loaded parts is limited 
(the height of the blades of the axial turbines must not 
be greater than 50-60 mm, the outer diameter of the 
working wheel of the axiradial turbines must not be 
greater than 180-200 mm); the turbine disks, consid- 
ering strength requirements, must be solid, which gives 
nse to definite design and technologic difficulties in con- 
necting the shaft to the disk, and also makes impossible 
the successive placement of several turbine stages on 
the same shaft; an increased negative influence of heat 
cycling with high levels oj the ume temperature deriva- 
tive 6/Ar predetermines use of aGVPP/RU primarily for 
bearing base loads. 


The mechanical properties of the CCM exert a signif- 
icant influence both on the design-layout solutions and 
on the specifics of a thermodynamic analysis for opt- 
mizing the GTPP heat circuits and cycles. 


The need for obtaining a high unit power of a GTPP as 
part of a GVPP/RU with a restricted discharge of the 
working medium through sections of the CCM feed- 
through parts leads to a unit-modular solution for a gas 
generator (Fig. 1). The working medium is directed from 
several modules with ceramic turbines (CTM) into a 
single metal pressure turbine (PT). In this case variants 
are possible either with a tree PT (solid lines in Fig. 
1) or with a general low-pressure compressor (LPC), 
driven by the necessary number of PT stages, which 
depending on the required frequencies of rotation may 
have either a common or different shafts with stages ot 
the free PT (dashed lines in Fig. 1). 


In turn, the structural layout of the CTM 1s two- or three- 
shaft with individual shafts for each turbine stage. The 
number of CTM shafts is determined by the admissible 
gas temperature before the PT, with an increase in which 
the number of shafts 1s reduced, but the design and 
technology of the PT blade crowns is more complicated. 
Thus, a simplification of one structural unit results in a 
complication of another 


In a thermodynamic analysis a study was made of four 
GTPP cycles: Brayton (B), regenerative with final re- 
generation (FR) with intermediate regeneration (IR) and 
with intermediate air cooling in the Compression process 
(IC). In contrast to the traditional thermodynamic anal- 
ysis of cycles, where as efficiency indices use is made 
of the efficiency 7, and the disposable work ||, when 
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using the CTM it is necessary to replace |, by the CTM 
disposable unit power N,.. This is attributable to the 
fact that it is impossible to ensure the supulated CTM 
power by the necessary discharge of working medium, 
as is possible in ordinary GTPP. Accordingly, cycles 
which without a loss of efficiency have higher values 
of the product of the degree of pressure increase and 
work |, are advantageous. In a comparison of the ther- 
modynamic efficiency of the cycles the unit power N.. 
must be regarded as conditional, since in its determina- 
lion use is made of some maximum admissible area of 
the critical section, which is determined after detailed 
design work. 


The dependencies n, = {( N.., derived with the partic- 
ipation of S. L. Yegorov, are illustrated in Fig. 2 for 
the cycles illustrated in the figure. When T_ = const the 
indices of the B cycle are dependent on only a single 
parameter: the degree of pressure increase 7, . The FR 
and IR cycles are dependent not only on 7, , but also 
on the degree of regeneration, assumed to be 0.85. The 
IR cycle is shown for the case of the most feasible posi- 
tioning of regenerative bleeding on the distribution line 
between the last stage of the ceramic turbine and the 
PT. 


The indices of the IC cycle are dependent on two 
parameters of the working process — the total degree 
of pressure increase 7, and its distribution in the 
compression process in the LPC and in the high-pressure 
compressor (HPC) (7, ,. and 7,4,.). For each 7, level on 
the lines 7, = f(N,, there are two characteristic points 
corresponding to such relations between 7, ,,. and 7. 
with which in one case 1, max 'S attained, and in the other 
case — N, _, Is observed. 


The principal advantage of a cycle with FR 1s a high 
efficiency. But with relatively low 7, values about 6- 
8 in the region 7, 4 low CTM unit power (N, = 
0.5 MW) is obtained, which increases their number 
Thus, the advantage of the FR cycle and a GTPP with 
a metal feed-through part (small 7, 1s transformed into 
an inadequacy for a GTPP with ceramic turbines as part 
of the GVPP/RU. A cycle with IR has a considerably 
lesser maximum efficiency than a cycle with FR, but 
higher levels of reasonable degrees of pressure increase 
(7, = 14-18), which substantially increases N, with 7, 
ow 10 both cycles. 
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Fig. 2. Dependencies 7, = {(N, of GTPP with T, = 1,300°C applicable for Brayton cycles (B) with final 
regenei ation (FR), intermediate regeneration (IR) and with intermediate air cooling (IC). 


Key: (1) B; (2) IR; (3) FR; (4) AC; (5S) m,.; (6) m; (7) MW 








Fig. 3 


Cycle B is used in gas turbine construction both 
autonomously and as part of a GVPP/RU. In a design 
with CTM, due to an absence of losses in cooling at T 
= 1,300°C, there is an appreciable increase in efficiency 
(up to 36.5-37 percent) in comparison with the best 
modern GTPP. GTPP with CTM operating in cycle B 
as part of a GVPP/RU make it possible to achieve an 
efficiency of 54-55 percent. 


A cycle with AC for GTPP with CTM is most attractive 
both in an autonomous variant and as part of GVPP/ 
RU units with different thermal and design schemes. 
The principal advantages of this cycle are: a high 
efficiency (41-42 percent), high degrees of pressure 
increase, which realistically may attain 35-50, which 
reduces the required number of CTM; a low-temperature 
heat exchanger for intermediate AC 1s far simpler to 
construct than high-temperature heat exchangers for FR 
and IR. 


When using a “classical” scheme (Fig. 3) in the makeup 
of the GVPP/RU a moderate 7, value was adopted 
(point a in Fig. 2) with degrees of pressure increase 
in centrifugal LPC, IPC and HPC equal to 3.7, 3.5 and 
2.5 respectively. The moderate value 7, = 32 with T. = 
1 ,300°C was governed by the need to have an acceptably 
high gas temperature at the entrance to the recovery unit 
for steam generation. The gas temperature in the exhaust 
is 515°C and the temperature of the live steam is 460- 
480°C. 
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Fig. 4. Thermal scheme of GVPP/RU with AC, with 
LPC and IPC drive by steam turbine: 1) LPC; 2) 
IPC; 3) condenser; 4) steam turbine with 
intermediate steam heating; 5) recovery unit; 6) 
metal PT; 7) two-shaft CTM unit; 8) intermediate 
air cooler. 


Key: (1) LPC, IPC, HPC (left to right); (2) CTM 


We emphasize that the use of a GTPP with intermediate 
AC in the GVPP/RU is possible only with a high gas 
temperature beyond the combustion chamber. Otherwise 
the temperature of the live steam will be low, which 
evidently requires additional conveyance of heat at the 
input into the recovery unit, a decrease in the degree 
of binarity (yn, < 1) and an efficiency n, in accordance 
with equation (1). In the case of autonomous use, in 
order to increase m, there are no restrictions on reduction 
of exhaust temperature. Accordingly, even when T, = 
1,300°C it is possible to achieve x, = 50 (for example 
Tp = 5.8-6.0 in the common axial compressor, 7,,. 
= 3.5 and 7,,,. = 2.5 in the IPC and HPC of a two- 
shaft CTM). The temperature before the common metal 
turbine is entirely acceptable and is about 1,040°C. In 
this case a 5 ‘te high efficiency will be attained n, = 
41.8 percent (point c in Fig. 2). 


We note sore peculiarities of the GVPP/RU heat 
system shown in Fig. 3. 


A steam turbine with intermediate steam heating, but 
without regenerative takeoffs; with a temperature of the 
combustion products at the PT output 515°C, n, about 
0.37-0.38. 


The flow of cold feed water at the condenser output is 
divided into two parts. Part of the discharge is directed 
into an air cooler, from which after heating it is directed 
to a corresponding heat exchange section of the recovery 
unit; another part is directed to the cold output (in the 
direction of gas movement) heat exchange section with 
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ceramic heat exchange surfaces not subjected to low- 
temperature Corrosion, 


The diagram shows a recovery unit of he feed-through 
type (the boiler also may of the drum type). 


The total efficiency of the GVPP/RU, designed in 
accordance with the diagram in Fig. 3, is 7, = 0.6, but 
the power of a gas turbine with intermediate AC is more 
than double the power of the steam turbine. The latter 
makes possible some modification of this scheme by 
using the PT not for the generation of electric power, 
but for driving the LPC. The output electric power and 
the efficiency in this case do not change, but there is 
an increase in the power of the GTPP PT by a factor 
of almost 1.5. A merit of the modified scheme 1s the 
installation of only one generator, but a shortcoming 1s 
that autonomous GTPP operation is impossible 


The following goals are pursued in another of the 
considered GVPP/RU heat systems (Fig. 4) 


Simplification of design by a changeover to a single- 
shaft CTM scheme, for which the gas from the ceramic 
turbine stage enters at a temperature about },150°C into 
the high-temperature section of the RU 


Reduce the te:nperature of the gas before the metal 
power turbine to 800-850°C 


Increase the temperature of the live steam upon entry 
into the GTPP to that adopted in modern condensation 
SPP — a temperature 540-560°C. This increases the 
efficiency , (without regenerative takeoffs) from 37-3 
to 40-41 percent, and also SPP power 


In Fig. 4 the heat system is shown in a variant of a LPC 
and IPC drive by a steam turbine, although it also can be 
designed in a variant similar to that shown in Fig. 3, that 
is, with a drive of the compressors (or one of them) by 
a gas turbine and drive of the electric generator proper 
by a steam turbine. 


A problem for thermal engineering equally important 
as an increase in the efficiency of the power units is, 
as noted above, the problem of reducing the specific 
discharge of harmful substances into the atmosphere 
(NO,, CO). Although there are no fundamental differ- 
ences in measures for reducing the content of NO, and 
CO in exhaust gases in high-temperature GTPP with 
ceramic feed-through parts in comparison with ordinary 
GTPP and SPP, this question requires independent de- 
tailed analysis. 
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Fig. 5. Heat system for GVPP/RU with 
stoichiometric fuel combustion: 1) LPC; 2) IPC; 3) 
two-shaft CTM unit; 4) combustion chamber; 5) 
steam heater; 6) metal PT; 7) recovery unit; 8) 
steam turbine; 9) condenser; _— intermediate air 
cooler. 


Key: (1( LPS, IPC, HPC (left to right) 


We note only that there is basis for expecting a 
possibility of reducing the content of harmful impurities 
by a specific organization of the combustion process 
in chambers constructed of CCM in comparison with 
ordinary metal chambers. 


It is possible to use heat systems in which there is 
stoichiometric combustion in the absence (or with a low 
content) of oxygen in the combustion products (Fig. 5). 
This also makes it possible to use “dry” purification 
methods, including conversion to inexpensive catalytic 
devices not requiring use of precious metals. 


Conclusions: 


1. The use of CCM for the fabrication of high- 
temperature uncooled elements of the feed-through 
parts of a GTPP will make possible a substantial in- 
crease in the efficiency of thermal electric power plants 
(by 15-18 percent), both on the basis of autonomous 
GTPP and especially in GTPP/RU. 


2. Due to known technologic restrictions on the size 
of loaded parts fabricated from CCM it is desirable to 
employ a block-modular approach to organization of the 
feed-through part of a GTPP. These restrictions also lead 
to a specific approach when optimizing the parameters 
of the working process of the thermodynamic cycles and 
heating schemes. 
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3. For GTPP it is desirable to use a thermodynamic 
cycle with intermediate air cooling in the compression 
process. 


4. GTPP/RU heating schemes with the stoichiometric 
combustion of fuel and using catalytic purification 
devices for exhaust gases, applying the conversion 
principle, are of interest but require further research. 


Method for Nonlinear Correction of Dynamic 
Characteristics of Rudder Drive of Spaceplane 
Using Standard Model 


957A1132A Moscow VESTNIK MOSKOVSKOGO 
GOSUDARSTVENNOGO TEKHNICHESKOGO 
UNIVERSITETA: SERIYA PRIBOROSTROYENIYE 
in Russian No 1, Jan-Mar 95 pp 38-44 


[Article by A. G. Borovkov, candidate of technical sci- 
ences, assistant professor, I. V. Misyuk, graduate stu- 
dent, A. G. Kalashnikov and M. N. Demenkov, Au- 
tomatic Control Systems Department, Moscow State 
Technical University imeni N. E. Bauman; the first para- 
graph is an abstract; UDC 62-82-585.2+62-52.001.5] 


[FBIS Translated Text] The efficiency of correction 
of the dynamic characteristics of a rudder drive are 
examined using information on the rate of the pressure 
drop in the pressure cylinder working cavities and a 
refined model of the drive load is presented. 


It is proposed that the frequency characteristic of the 
rudder drive of a space vehicle be corrected, and 
specifically that the resonance peak at frequencies 
corresponding to the frequency of fluid oscillations in 
the fuel tanks be reduced. Clearly expressed resonance 
phenomena arise as a result of the weak characteristic 
damping of the rudder drive (RD), as well as the 
low damping constant of the kinematic circuit in the 
neighborhood of the characteristic frequency of the 
elastic circuit. The presence of low-frequency resonance 
in the rudder channel is inadmissible because in this 
case the characteristic frequencies of oscillations of the 
elastic contour and the elastic oscillations of the rocket 
body converge. 


A special feature of the dynamics of a rudder drive for 
a class of flightcraft, control of whose thrust vector is 
accomplished by means of the jet engine drive, is the 
great moment of inertia of the rotating engine; limited 
rigidity of kinematic circuit elements; low damping 
constant; nonlinearity of the friction charactcnstics. 


Accordingly, for stable RD operation it is necessary, on 
the one hand, to ensure a stipulated phase lag which 
is attained for the most part by an inciease in the 
drive quality factor, and on the other hand, to ensure 
suppression of resonance at the frequency of the elastic 
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circuit, which is amplified with an increase in the quality 
factor. It follows from an analysis of the mentioned 
requirements that without additional correction devices 
the simultaneous satisfaction of these requirements is 
impossible. 


One of the methods for solving this problem is the 
introduction of a correcting feedback in the rudder 
drive scheme for the derivative of the pressure drop 
in the cavities of the rudder drive pressure cylinder — 
hydromechanical correction (HMC) [1]. The essence 
of HMC is based on passage of the working fluid 
from the pressure cylinder high-pressure cavity into 
the low-pressure cavity at the time of appearance 
of oscillatory processes. For a digital RD with a 
mechanical feedback the use of HMC makes it possible 
to use the positive properties of the analog internal 
feedback for correction of the RD-load circuit. A 
shortcoming of hydromechanical devices of such a type 
is their operation in the LF range of the amplitude- 
phase-frequency characteristic, an increase in structural 
mass and size. 


The formulation of the problem involves solving the 
problem of correcting the rudder drive-motor elastic 
circuit by digital correction (DC) devices. A special 
feature of this system (Fig. 1,a) is the absence, due 
to design complexities, of a technical possibility for 
measuring the load coordinate, that is, the angle of 
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rotation of the longitudinal axis of the engine jet relative 
to the rocket body. In the formation of the DC signal use 
is made of information on the intermediate variables of 
State of the system, on the pressure drop in the pressure 
cylinder cavities and movement of the pressure cylinder 
rod (Gray code). 


There are different methods for solving problems of 
this class. A method for organizing correction by 
use of a standard model, taking into account the 
acuion of dissipative forces in the kinematic scheme 
and the influence of the nonlinear segment of the 
characteristic curves of frictional forces in the bearings, 
was proposed in [2]. The use of a standard model makes 
it possible to discriminate load information criteria. 
An algorithm for control of the system on the basis 
of the output coordinate was described in [3], where 
the necessary accuracy was attained by curve fitting 
of the control algorithm coefficients. On the basis of 
the collected information an algorithm was written for 
digital correction of the rudder drive and this produces 
the correcting signa! fed to a 12-digit summator, thus 
performing system signal correction. 


Mathematical model «* loaded RD with digital correc- 
tion. Figure | is schematic diagram of the RD-load elas- 
tic system. The rudder drive 1s a digital electrohydrome- 
chanical tracking drive with three amplification stages. 
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Fig. 1. Schematic (a) and structural (b) diagrams of rudder drive with correction. 


Key: 1. U,; 2. summator; 3. U__; 4. DAC; 5. Injection; 6. Discharge; 7. Microcontroller; &. Digital manometer, 9. 
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A 12-digit digital-to-analog converter (DAC) is used as 
the first amplification stage; a 7-digit binary code (the 
lowest 5 digits are used for correction) is fed to the DAC 
input and the mechanical movement of the working rod 
is read from the DAC output. The next two amplification 
Stages are represented by a slide valve hydrodistributor, 
in each case with a rigid mechanical coupling. 


Information on movement is read from the pressure 
cylinder rod by means of « system of digital sensors of 
linear movements and together with information on the 
pressure drop in the pressure cylinder cavities, received 
from a digital manometer, is fed to a computer which 
constructs a mathematical model of the load, on the 
basis of which a nonlinear feedback signal is shaped 
which is proportional to the power gradient, which is 
fed to the DAC input. 


The dynamics of the considered system is described by 
differential equations of mouon with the following as- 
sumptions: the DAC is a net lag component with a time 
constant T,, = 0.015 s; the injection and discharge pres- 
sure is constant; the elastic modulus (hydraulic stiff- 
ness) and viscosity of the working fluid are constant 
and are not dependent on number and temperature; hy- 
drodynamic forces in the second and third amplification 
Stages are not taken into account; the masses of the slide 
valves, pistons and pressure cylinders are neglected due 
to their smallness in comparison with the mass of the 
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load; RD operation in the neighborhood of the cavitation 
inodes and the phenomenon of supulated rod position 
drift are not described. 


With allowance for these assumptions we write the DAC 
equation X_, = I[t- 7,JT,,.(U,/O,,)A,,, where r,, 
is the lag lume constant; A. is pac quantization — 
DAC rod movement in response to the sending of a 
single control signal; T,. is a function for discarding 
a fractional part; the equation describing the second 
amplification stage 1s 


Kemp = Xi > Ky Xanp/K,, OF TX ump + amp = KX. 
where K,, is the amplification factor for the speed of 
Stage 2; K, , 1» the kinematic feedback factor; T. = 
; K, = I/K, , (here and in the text which 
follows the asterisk indicates that the preceding letter 
is surmounted by a dot and two asterisks represent two 
dots]; the equation for the kinematic circuit relating the 
coordinates of the piston and load is 


CgAx = Sip - F,, sign (x*), where Ax = x, - x; 
S, is the area of the pressure cylinder piston; Ap is the 
pressure drop in the pressure cylinder cavities; F,_ is the 
dry friction force on the rod and piston seals; C,, is the 
reduced stiffness of the working piston. 


In turn, 


(1/ Cy +1/ Coop +1 / Cop +1/ Cy +1/ Cy): 


1/ Cy +1/ Cyogg +1 / Cop 1 / C0, 


subscripts end superscripts recur in th sussian 
ollew. Their renditi~-n in the English tres ation is 


co = fb} BX = ins wp red; T= po HTTP = fr/ti 


a 
coat © bedy; MILT = rod; # = load; K = att; 3 = valve; Ty = fr) He £1) 


here C.. C._, are the rigidities of attachment of the 
RD to the support and the load; C_,, C_,, C,, are the 
rigidities of the RD body, loaded part of the rod and 
ihe hydraulic spring, governed by the compressibility 
of the working fluid in the cavities of the RD pressure 
cylinder; x. is the coordinate of the valve of the third 
amplification stage; A. is the positive valve coverage 
of the principal hydroamplifier. 


cor; CXC @ sys. The meaning of these and other 
ev 


The equation of motion of the load m x**. =F’ - 
Fir (FP, x* uy) Po = Clan = KO, - F., where 
C., is the reduced rigidity of the RD, K. is the 
position force coefficient; F, is static force, F,, is a 
nonlinear function describing the dependence of the dry 
and viscous friction levels on the total (external, elastic 
and static) force F° imparted to the load and the velocity 
of load movement. 





rr report may contain copyrighted material. Copying and dissemination 


is prohibited without permission of the copynght owners 

















f 


FLF <\.3F py iq =O: 


br, +13 Fp - K.,arctg(K .x, ). 





Here b is the viscous friction damping factor; F, is dry 
friction force; F,, = (2/m)0.3F, is the dry friction factor; 
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if 


in ™ rf - = sign (23 
0 


where K, is the RD velocity amplification factor; p, 
is the injection-discharge pressure difference; x |... is 


the effective coordinate of the third amplification stage 
valve. 





C. 
— a 2 p= Sem ay | ~ SEo +——r—v,, 
C. +C, Py 


op 


where v,_, is the velocity of load movement, 


The follow ag restrictions are imposed on the coordinate 
x 


valve ce’ 


‘0, x4] < 43: 
Xyq * }%3 — 4; sign (x3), A; < [xs| < 23 max: 
(x5 maz ~ 45) sign (x3), [25] 2 25 max 
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K, is a coefficient characterizing the dependence of dry 
friction on the velocity of load movement. 


The velocity equation for the RD third amplification 
Stage is written in the form 


sign(t - sign (x34)) - “2 ap? 


Py r 


Using the equations for RD-load kinematic coupling, 1t 
is possible to transform this equation. Since 


Cc F 
Ap = —~ Ax + —P sign (x,), To 
Sy Sy 


and neglecting the derivative of the second term, we 
obtain Ap* = (c_/s,)Ax*. Then 


Py r+ Cay 


Here x... 18 Maximum valve motion. 


The equation for the third stage feedback has the form 
Kae = (1 + Kx, - K,x,. The x, coordinate can be 
found from the x, coordinate: x, = (c_/c, )Ax + xX... 
where c, = (I/c., + Wc) '. 


The system of sensors of linear movements is described 
by the equation x, = T, (x /O,)A.. where A. is 
sensor discreteness. 
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In the mathematical model the load will be described 
by a simplified differential equation 
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ME + Ding, + Ky Zing + Fo + Fp S190 (Xoyy) = Capy (Farr Xm) 


{In the above equation HM = simulated load] 


Jumplike deviations of the x,_, coordinate, having a ran- 
dom character, may be observed under the influence of 
aerodynamic perturbations and nonuniform fuel com- 
bustion. It is therefore desirable to obtain this coordi- 
nate not by mathematical simulation, but directly from 
the system. This problem can be solved by installing a 
digital manometer connected to both pressure cylinder 
cavities. Then 





» ow _ AP ty — Frege 81D (Fee) 
tm * or ~ Com 


where r._ is the lag constant for the correcting device 
(integration time of the mathematical model equations); 


0.| Kea - Fx| > Ay: 


@ 


Usop 





( te | Ket - Xie )/ Kea 


where Ap is information on the pressure drop in 
the pressure cylinder cavities, read trom the digital 
manometer. 


It must be noted that in this case there is no need to de- 
termine the velocity of the rod — only its direcuon need 
be determined. Since we will obtain x*,,,, directly from 
the system in real time, the order of the mathematical 
model equation will be reduced to |, and this reduces its 
integration time in the correction process. The equation 
for the principal correcting relationship has the form 


= Trunc Upon / Su )4y It - tx 


here A. is the magnitude of the insensiuvity zone within 
which no correction is made, K,, is the static system 
transfer factor, K,, is the feedback coefficient. Figure 
1,6 is a block diagram. 


Study of the equations of motion using digital correcting 
device. The system of equations was studied by tts in- 
tegration by the fourth- fifth order Runge-Kutta method 
using a computer with the following RD-load system 
parameters: 


m = 579 Hc? cu; b = 260 Hem”; Fy = 0; K, = 1890 Hem’: 
Koc, 0,5; Kee = 0,166; Ky, = 475 c°';  p, = 2100 Hem”: 


K, =81,3c"'; 43 = 0,005 cm; 


Z3 max = 0,154 cm: K, = 500: 


Fy = 2000H: c, =18-10° Hem; cy =11-10° Hem: 


Cop 2 40:10° Hew; cc, = 25-105 Hem: c 


F 


and the digital correction unit parameters: 4. = 0.0001 
em; K., = 50,000; 7, = 0.02 s; A, = 0.01 cm. 


= 44.10° Hem”. 


pn @ 22000; S, = 157.6 cu?; 4, = 0.000625 cm; t, = 0,015 ¢ 


The equation for the load model was integrated with 
the interval 7. and the remaining equations with the 
interval TS 7. 





is report may contain copyrighted material. Copying and dissemination 
is prohibited without permission of the copyright owners. 














22 SPACE SCIENCE AND ENGINEERING 


The system was studied with the sending of a command 
Signal in the form of a step constituting 4 percent of the 
maximum rod stroke (0.05 cm), which corresponds to 
the amplitude of the command signal during in- flight 
rocket stabilization. 


Figure 2 shows the transient process of the RD-!oad 
system with correction "on" and “off.” The presence 
of a scatter of steady values between the corrected 
and uncorrected system is attributable to the existence 
of a high level of dry friction, which gives rise to 
an insensitivity zone. Figure 2 also shows a corrected 
command signal with allowance for lag of the DAC and 
correcting device. The graph shows that the control is in 
antiphase with the oscillations of the uncorrected system 
and has some lead relative to it. 























6,/ ~ 


























| w 
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Fig. 2. Time dependence of velocity v,_, with 
correction (2) and without correction (1); (3) shape 
of correction signal. 


Key: 1) s 


It must be noted that with the sending of sinusoidal 
command signals to the input this correction exerts a 
weak influence on the phase and amplitude of the output 
signal up to 1.5 Hz and therefore there will be assurance 
of the phase lag necessary for stable rocket motion at 
the frequencies of escillations of the fuel components 
in the tanks (0.5...1 5 Hz). 


Possible field of applicability of correction. The pro- 
posed correction variant can be regarded as an alterna- 
tive when correcting the drive of the Energiya booster 
described in [1]. 


Nonlinear correction with a feedback, proportional to 
the rate of change in kinetic energy, can be used in 
both digital and analog forms and makes it possible to 
solve problems in the dynamics of actuating systems 
with allowance for requirements on the smoothness of 
motion of the controlled object when there is a lag in 
the direct circuit. 
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[FBIS Translated Text} Designs, equations of motion 
and an analysis of errors, speed and noise immunity 
parameie-s are proposed for surface gyrocompasses 
constructed on the basis of gyrotachometers, finding 
wide use in systems for azimuthal orientation of surface 
objects employed in many branches of the economy 


The modern development of technology raises problems 
in the development of systems for the azimuthal on- 
entation of surface objects (SO) for both defense (ar- 
mored and reconnaissance vehicles, “mobile” missiles, 
etc.) and civilian (tunneling equipment for shafts and 
tunnels, bottom work, objects situated within geophys- 
ical, oi! and gas wells, etc.) purposes. The solution of 
such problems requires gyrocompasses (GC) having a 
quite high accuracy and requiring only a short time for 
seeking the meridian (up to 8 minutes), small volume 
(up to 3. dm’ and a low mass, labor input and cost. These 
requirements are not met by the gyrocompasses used 
on seagoing vessels and spacecraft and the azimuthal 
orientation systems employed in aviation [1}. Accord- 
ingly, work has been undertaken in proposing, analyzing 
and developing other methods for gyrocompass orien- 
tation based, for example, on measurement of apparent 
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deviations of gyroplatforms and free gyroscopes, mea- _ horizontal reduction of platforms and gyroverticals and 
surement of information signals from the channels for — use of gyrotachometers (GT) [2]. 



































. Diagram of gyrocompass on a horizontal platform: A — 1, 2) accelerometers; 3) axis of azimuthal rotation 
of platform with GT; 4,7) actuating motors of system; 5) moving GT component; 6) GT body; 8, 9) 
mechanical device for reading course angles; 10) actuating motor with tachogenerator of system for 

azimuthal rotation of platform; 11) electromechanical device for reading course angles; B — switch of 
gyrocompass orientation signals: (at bottom — geometric, at top — analytic). 


Key: 1. computer; 2. V_.; 3. V_; 4. A; 5. moment sensor; 6. angle sensor 





This report may contain copyrighted material. Copying and dissemination 
is prohibited without permission of the copyright owners 














u SPACE SCIENCE AND ENGINEERING 


Systems for azimuthal orientation of SO constructed on 
the basis of GT have proven to be the most universal. 
We will examine them in greater detail because infor- 
mation published on them is of a fragmentary character 
(3, 4). 


In order to be specific, by the term GT we will mean 
a precision gyrounit (GU) with two degrees of freedom 
whose moving component has a feedback. Such a GT, 
installed on a base fixed relative to the Earth, strives 
to match the vector of its kinetic moment (virtually 
coinciding with the rotor axis (RA) of the gyromotor) 
with the vector of the horizontal component of the 
velocity of the Earth’s rotation, lying in the plane 
of the meridian. Theoretically, the GT also may use 
any other type of gyroscopic device for measuring the 
angular velocity (laser gyroscope, dynamically tunable 
gyroscope, etc.). 


The functional-kinematic layout of a GC constructed 
on the basis of a GT for the most part is determined by 
two methods: use of information read from the GT or 
minimizing of the GC error caused by nonhorizontality 
of the GT input axis. This error is generated by the 
vertical component of the velocity of the Earth's rotation 
and is already known from the theory of a Foucault GC 
as the inclination error. The figure illustrates an example 
of the functional-kinematic layout of a GC constructed 
on the basis of a GT. 


With respect to the method for using information, GC 
schemes are divided into geometric and analytic. In 
geometric schemes the GT output signal is used in 
reducing the GT RA to the plane of the meridian, after 
which the RA, like a magnetic needle, indicates the 
direction to the north. In analytic schemes the GT output 
signal is fed to the input of a computer, which using 
some algorithm computes the SO course angle. 


The inclination error also can be minimized geometri- 
cally and analytically. In geometric minimization the GT 
input index is reduced tc the plane of the horizon and 
thereafter is held there. This can be realized by mount- 
ing the GT either on a platform which is held in the 
horizontal plane by means of tracking systems whose 
sensing elements are accelerometers or on a physical 
pendulum. In analytic minimizing of inclination error a 
platform or pendulum is not required, but the GT body 
is attached directly to the SO body and accordingly du- 
plicates all its inclinations. In this case the GC has large 
inclination errors. In order to lessen them, acceierome- 
ters are mounted on the GT body. The signals of these 
accelerometers, containing information on the angles of 
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SO inclination, are used for compensating for the incli- 
nation errors. We note that the attachment of the GT 
body to the SO body, simplifying GC design, reduces 
the noise immunity of GC readings because the GT also 
senses the oscillations and vibration of the body. 


The examined GC schemes can be constructed using 
standard-produced precision GU, accelerometers and 
computers, which is a substanuial merit of such GC. 


The principal operating characteristics of GC con- 
structed using a GT are as tollows: time required for 
reaching a working mode, that is, the time expended 
On performing servicing Operations; ime required for 
determining the course angle; error in determining the 
course angle; noise immunity of GC readings. 


The time required for reaching a working mode 1s de- 
termined for the most part by the setting of the nominal 
GT thermal mode. Computations and tests of a GC con- 
structed using a GT indicated that with a negative am- 
bient temperature the power of the built-in system for 
thermostating a precision GU is inadequate for bring- 
ing the GT to a nominal thermal mode. Accordingly, it 
is necessary either to develop an additional system for 
thermal protection of the GT or a system for algorithmic 
compensation of GC temperature errors. The problems 
arising in this Case were examined in [5]. 


The remaining operational characteristics are deter- 
mined on the basis of a GC mathematical model. We 
will construct such a model for the GC scheme illus- 
trated in the figure on the assumption that the GT input 
axis 18 geometrically held in the plane of the horizon, 
the GT has reached a nominal thermal mode and the 
dynamic processes in the electric circuits and electronic 
units can be neglected in comparison with mechanical 
systems. 


We introduce into the examination a base system of 
geographic coordinates (CS) O&n C referenced to the 
Earth and right-handed CS Oxyz, referenced to the body 
of the SO, GC body, platform, GT body and GT moving 
component. We stipulate the relative angular position of 
these CS as a result of angular movement of the first 
kind, that is, Euler-Krylov angles are used. 


Using the introduced CS, we obtain (with allowance 
for quantities only of the first infinitesimal order) the 
following proyections of the absolute angular velocity 
onto the axis of the GT moving component and onto 
the platform z-axis: 
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@, = @, sin(y + ¥,) + (ay + w)(a, sin ¥_ - By COS 7, ~ By - B,) + 
+0, COW + Ya Yen + Yun + Yu + Yr) + (DNCOB 7p — (Yun + Yun + Yu + 
Ye) 8in Yq) + Risin Yq +! uy + Yun + Yq + Vp) C08 Yq) + Gy + (7, + ¥,) By - 

-(6, + B,) Yu+ a, 


y 


@, = Wy COS(W + Yq) + (We + WB, SIN 7, + Aq COB Yq + Oy * a) a em + 


+ Ye + Yun + Y,) sin( yw + ¥_,) - d(sin 7, + (Yun * Yun * Vx + Y,) cos Yo ” (1) 

+ teos ¥_ — (Yen + Yun + Yq + Yr) SIN Yq) — (7y + Yr) @q + By + (a, + @,) ¥_ + B,: 

w, = w,((B, + B,)sin(y + 7,) - (a, + a,)cos(w + y,) + B, sin yw - a, cos w] + 
+ + w + O1(B, + B,)c0s ¥q + (a, + @,)siny_ - (x - By) + mB, + 

+B,)sin ¥_ - (a, + @,)cos y, +(0-a,)]+ (A, +B,)a, -(a, +a,)8,+ Yq 7, 
w! = @, (A, sin y - a, cos y) + 0(8, - K) + K(d - a,) + we + Wt Yq, 


Bote: The following subscripts and 


ipts recur in the Russiar 


supersecr 
formulas which follow. Their rendition in the English translation is 


indicated: 


fr @=gs we pl; Ke co; KW @= fin pl; HT = in pl; 2M = ms; Ay = as 


Rea Be ;TT = Ty He TT; 


these and 


where 


®, = U cos 9, @ = Using 


are the projections of the velocity U of the Earth's 
rotation onto the axes of the geographic CS; ¢ is 
the latitudinal position of the GC; 4, «, ~ are SO 
pitching, banking and course angles; o, /3 are the angles 
of deviation of the GC units (to which the CS are 
referenced) from the plane of the horizon and 7 — in 
the plane of the horizon; 4*, «*, * are the angular 


K ry2, _ 
= J, = Dy ¥_% Kyl — Hoos - yo, + (Jy - Js Ml@y) 7 (o, 


+ wr@,|+M; = 0; 


where J, J.*, J,‘ are the axial and equatorial moments of 
inertia of the moving GT component, D, is the damping 
coefficient of the GT moving component; K_, is the GT 


moment sensor transfer coefficient; i is current strength 


pH = in dG; H = in. The meaning of 
abbreviations is evident from the text. 


velocities of the SO body caused by force impacts acting 
upon it (wind loads, ground tremors, ground subsidence, 
etc.) [here and elsewhere in the text the asterisk indicates 
that each of these Greek letters is surmounted by a dot); 

a B*.. Ty» a*, PP, T". are the angular velocities 
of the platform and the ‘GT moving component caused 
by external and internal force impacts acting upon them. 


The GC mathematical model includes three groups of 
equations. 


The equation for the GT moving component (with 
allowance only for quantities of the first infinitesimal 
order) is 


ry + 


(2) 


in the GT feedback circuit; H is the kinetic moment 
of the gyromotor; M.* is the total perturbing moment 
Operative around the GT output axis. 
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The equations for the GC electric circuits are 


| “lp . ; 
Upy * Udy ~ Kye Uys Ur + KyUyyi f= Uy Re i Re = Re + Ry’ 
U, #iR;U, = Uy_-Ky (WU, +U%U,, = KU, U, =U,+ (3) 
+ Kopf Yq: 
where ator; K., K,,, K,,, Ke K,,, are the transfer coefficients 


for the angle sensor, amplifiers, strength of the feedback 
regulator for platform reduction and tachogenerator; R., 
Uy, Up, Uyy Use: Ven Ury: Uy. Usa. Ur R,, are active electric resistances of the GT output re- 
sistor and the control winding of the moment sensor; n 

is the reducer transmission ratio. 


Uy. 


are the output and spurious signals of the angle sensor, The equation of motion for a platform with a GT is: 


amplifiers Al and A2, gyrotachometer and tachogener- 


-J, @'+nM,-M,=0; S,ny.+M, = K,Uy). (4) 


where J, is the moment of platform inertia with  §, is the motor rotor damping factor; K_ is the motor 
allowance for the reduced moment of inertia of the transfer coefficient. 

motor rotor; w," is the absolute angular acceleration of 

the platform, M,, is the moment created by the motor; Differentiating expression (1) in time for w,” and w5, 
M,, is the moment of resistance of platform supports; — we obtain 


a= +7, + Of, - K)+t(0-a,)+ 0(6,- «)+ 
+ i(O - @,) + @,[(8, + a, W)sin w - (a, - fA, y)cos wi (S) 
o;=y+7,+7, + HB, +B,)cosy, + (a, + @,)sin y, - (x - B,)|+ 
+ (By + B,) sin Yq ~ (ay + c)08 75; + (0 - coy)]+ Gy By + By) - 
~ By (a, + a) + Ha, sin 7, +B, cos 7, ~ (k-B,) + Fully + &,)c08 7, - 
~(B, + B,) sin 7,1} + (A, sin y, - &, 008 ¥_ + (0- Gy) + 7qllay + a,) x 
x sin 7, + (A, + B,)c08 yn) + Oy By - By & + @y(, sinly + ¥_)- 
~ Gr coslw + ¥,) + (w + 7, M(B, + A,)cosly + 7.) + (a, + (6) 
+ a@,)sin(y + 7,)]+ ia, sin y +B, cos w) + B, sin y - a, cos y). 
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Substituting expressions (1), (5) and (6) into (2) and (4) 
and excluding the intermediate variables in equations 
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(2)-(4), we obtain the equation of motion for the GC 
platform: 


(Jy + Kink, KyRJ,) igt US, + Ky KyaKyKep)n” ~ Kink, Kyg x 


x RH(B, +B.) ¥, + Kj AK, KyR, Ha, sin(y + 7,) = M,, 


where M_,, is the periurbing moment. 


In a general case the analytic expression for M_ is 
unwieldy because it is written on the basis of the 
unwieldy expressions (1), (5), (6). We will examine the 
principal special cases. 


Ideal case when M.. Then for a steady platform 
position we have 


sin(y + ¥,) = 9. (8) 


Equation (8) has two solutions 


(9) 
(10) 


V+7, = 9% WM =~ Va: 
Ww +7, =180°; ww, = 180° - 7, 


Solution (9) corresponds to a stable equilibrium of the 
platform and solution (10) corresponds to an unstable 
equilibrium ("spurious zero”). In GC schemes with con- 
tinuous platform reduction the spurious zero problem is 
not pertinent because the slightest perturbation of the 
GT moving component (as is always observed in prac- 
tical work) causes a triggering of the tracking system 
and platform reduction to a position of stable equilib- 
rium. 


An analysis of the coefficient on 7*, shows that 
depending on the value and signs on the random angles 
8, and f, it may be both positive and negative. In the 
latter case the motion of the platform will be unstable. 
The reason for this is the possible angle of deviation 
8, + 8, of the GT output axis from the platform axis 


l 
sinly + 7,) = w+ 7, = ——(M, - 
Ho, 


(7) 


of rotation. In this case the component of angular 
velocity of platform reduction (1. + § )y*,, is projected 
onto the GT input axis. As a result the reacts not 
only to the useful angular velocity w, but also to the 

urbation (3, + J, )y*,. The nature of platform motion 


is dependent on the ratio of the values w, and (/3, and 
BY: 


In experimental studies of GC an unstable platform mo- 
tion was detected in the form of increasing amplitudes 
of platform oscillation about the plane of the merid- 
ian, passing into continuous platform rotation with an 
increasing velocity. This phenomenon has been called 


“gyroscope capture.” 


The coefficient on -y*.,, and therefore also the nature 
of platform motion, may change during the time of 
GC operation due to a change in the latitude ¢ of 
GC position and change in the angle 9. + J. A 
differentiating link must be introduced in order to ensure 
the stability and quality of platform reduction system 
operation. A DG-1 actuating motor-generator was used 
as such in the constructed GC variants. 


On the basis of solution of the nonlinear differential 
equation (7), found, for example, by the so-called 
“alinement” method, it is possible to find the time 
required for platform reduction to the plane of the 
meridian with large initial angles ¢ + 7,, < 180°, that is, 
the time required for determining the course angle. The 
minimum experimental time for determining the course 
angle with an initial angle of deviation of the platform 
from the plane of the meridian by 179° was about 3 
minutes. 


In the first approximation we find the error in determin- 
ing the course angle. For this purpose we will assume 
that the GC is in a steady position (that is, a*, = /* 
=7,=a* = = *+* = ()), but tremors of the SO 
body are absent (that is, § = « = u* = 0). Under these 
conditions the expression for M,, is simplified and on 
the basis of equations (1)-(4) we obtain 





> ast 
Kan o KU y2 o Kase re | 
OKA Kak Kak 


+ (PB, +B, — Gq SIM 7_ +B, 605 ¥n 1G — (Yen + Yun + Me + Ve COSIY + Yq). Ni 
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An analysis of expression (11) reveals that a GC with 
continuous platform reduction has no methodical errors: 
all the error components are instrumental and have a 
quite clear physical sense. We note that the error is 
dependent on the latitude of the place of GC operation 
and increases with an increase in latitude. 


The GC readings are read when its platform is reduced 
with a smali error (11) to the plane of the meridian. In 
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equilibrium and a controlling signal is sent platform 
reducuion system, which begins to perform damping 
oscillations relative to the plane of the meridian. The 
sending of readings from the GC is interrupted. This is 
the essence of the poor noise immunity of readings of 
GC constructed using geometric schemes. 


We will find analytic expressions for possible perturbing 
influences on the GC at the time when readings are sent 


this position the directional gyroscopic moment Hw(d _ from it. For this we obtain an expression for M.. on the 
+ 7, is small. As a result, even weak perturbing effects assumption that the GC is ideally constructed 
drive the movable GT component from a position of 
IM, = M, = 0. 
US = US = 0,0,= Bp = te = Oe = Be = = Yan = Man & OD) 
the platform is in the plane of the horizon (a, = 4, = With allowance for these conditions we find w,", w,§, 


0) and that the banking angles « and pitching angles 4 
of the SO are equal to zero. 


Jn Wnt (S, n? +n°K aKyaK, Ka) int Ky nk, Ky,R.Ha,(y + y 
a+ @ 7(B, sin y - a, cos y)}- K7! 


ond alW + O(8, - &)+ K(p- 


w,*, w* and we substitute them into equations (2) and 
(4), AS a result, we obtain 


le 
nk, + 


-[ 


. K RAJ, [w+ 75+ ¥-+ iis: sin Yq +B. cosy, -K + B,) + X(8, si siny, - 


- & cosy, + 0 - 


a,) + @, B,- @, B,- w,(a,- B,, sin y + (12) 


+ a, cos w)! + D, i, + H(8cos 7, ¢tdny, + @, + a.)- 
(Jy - JEM Bcos 7, + K sin y, + G, + @Miccos /,. - P sin 7, + B+ B.D) 


The right-hand side of equation (12) characterizes the 
perturbing influences, each of which may drive the 
GC platform from a stable position, and shows that 
the perturbing influences on the GT and the platform 
were caused by the angular acceleration of any of the 
moving GC components (7**,, 7** ) and SO angular 
acceleration (y**), as well as virtually any combination 
of SO angular accelerations (6*, «*) and moving GC 
components (a*,, A, a*, is {here the two asterisks 
indicate that two dots surmount. the accompanying 
Greek letter]. 


The impact of high-frequency harmonic effects can be 
reduced somewhat by careful choice of parameters for 
the GT feedback and platform reduction system, but it 
is virtually impossible to lessen the influence of low- 


frequency and pulsed perturbing influences. 


GC schemes with time separation of GT signal measure- 
ment and platform reduction processes were proposed 
and studied for the purpose of lessening the shortcom- 
ings characteristic for GC schemes with continuous plat- 


form (P) reduction. We will examine one of these as an 
example. The GC scheme shown in the figure will be 
supplemented by an integrator (1) and controlling unit 
whose purpose is setting of the time intervals and per- 
forming switching of the electric circuits. The cyclo- 
gram of GC operation contains three subcycles (SC). 
The first subcycle is a pause necessary for the damping 
of all transient processes in the GC systems. In the sec- 
ond SC the GT output is connected to the I input, the 
GT output signal “is accumulated” at the | output, no 
controlling signal is fed to the actuating motor and the 
platform is fixed in position; the second SC ends upon 
the elapsing of the time T, of the second SC. Then be- 
gins the third SC — the GT-I connection is broken, the 
I output is connected to the control winding of the ac- 
tuating motor, I is discharged and the actuating motor 
rotates the platform, bringing it close to the plane of 
the meridian. The third SC continues for the duration 
of one cycle. Then the cycles are repeated, as a result 
of which the platform is discontinuously reduced to the 
plane of the meridian. 
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The motion of the GC platform in one of the cycles is 
described by the system of equations: 
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-J,@8+nM,-M,=0; S,nj,+M, = KUy Uy= 


= Kh + Uz 


where U,**, U ,*, U,™ are the | initial, current and 


spurious discharge voltages; T,* is the discharge time 
constant; K, is the I transfer coefficient. 


J 9, D, = = 
K,yKy, K 


where N ., is the deflecting moment operative around the 
GT output axis and in a general case containing dynamic 
perturbations and inertial moments) and 


(gyroscopic 
quasistatic perturbations of the weighting moments type 
(Ms. 
Assuming that N., = 0, using formulas (13) and (14) 
we obtain the equation of motion for the platform in 
one cycle: 


Js ¥_* n’s, ) - aK, Kyu —"/% - M,. (15) 


By solving equation (15) we find the angle and the time 
required for reducing the platform in one cycle: 


yoo AT | -aets (16) 


n’s, 


(17) 


opin 


: UP umes, 


T, 
UM =U; + Ky {U, at (13) 
0 


We find the GT-U, output signal by using expressions 
(2) and (3) and in the final expression retaining only 
terms of the first infinitesimal order: 


U,+ KU, = RLHa, sini + y,)+ Nel, (14) 


Using expressions (16) and (17) we will construct a 
graph of platform motion from which it is possible to 
judge the stability of GC operation and the time for 
determining the course angle. 


We find the error in determining the course angle on 
the assumption that the platform has reached a steady 
position (7**, = 0, 7*,, = 0). Then from equation (15) 
we have: 


nK Kye Me 


"-M, =0, (18) 


where t° is the specific moment in time for onset of the 
equality (18). 
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Expanding equality (18) with allowance for expressions 
(13) and (14), we obtain 
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y+ In - J Ka Me Kalle at |, 
Ho, | aK, Ky,K, RT, K,RT, j 


+ (B, + B, - Gy sin 7, +B, 008 ¥,)t9 O— (Yen + Yun * Yq * Yr) COS(W + Yq ).(19) 


It follows from a comparison of formulas (19) and (11) 
that with a correct choice of the parameters K, and T, the 
error in GC schemes with the separation of processes 
is somewhat less than in GC schemes with continuous 
platform reduction. 


Schemes with discontinuous platform (P) reduction 
make ible a considerable increase in the noise 
immunity of GC readings. The interference signal 
accumulated at the I output is not always capable of 
causing the triggering of the platform reduction circuit. 
We will show this for the case of a pulsed impact on 
the gy WY component, most probable in actual 
practice. solution of equation (14) in the case of 


a pulsed impact has the form 


UP = Ute sin(ant + ?,.). 


(20) 


T, 

' K.u" | ?.’ \ 
Uy = K,, fura = ~ 1o@ "(1 snar, sonar} 
+o; \ or, 

0 f . 


In a correctly designed GC T, >> T,, w >>\ 1/T,, then 


Us - (Ky, / w)UP. (21) 


Expression (21) shows that with proper choice of the 
K, and w quantities the noise signal at the I output can 
be considerably reduced, and accordingly, the noise 
immunity is increased. 

In addition, the GC schemes with discontinuous P 
reduction give a fic possibility in these schemes 


for cutting off the P reduction system at the end of this 


cycle when the I output signal is equal to zero. In this 
case P does not have a possibility for moving relative to 
the SO body and its noise immunity is virtually ideal. 


where U.™, U,” are the current and initial levels of the 
GT interlerence signal caused by pulsed influences; w 
is the frequency of the characteristic oscillations of the 
GT moving component; ¢, is the initial phase of these 
oscillations; T, is the time constant of the interference 
signal damping exponent. 


The U.™ signal is fed to the input of I, at whose 
output we will have a spurious controlling signal U,™. 
Assuming that U" = O, ¢,.. 4 0. omain 


Fr 


However, when the reduction system is cut off errors 
may arise in the GC readings due to both to rotation 
of the SO body in the azimuthal plane and a random 
compensation for useful signal interference at the | 
output. Fundamental solutions for the detection of these 
errors and their elimination are simple. 


Among the shortcomings of GC schemes with discontin- 
uous P reduction it is necessary to include an increase in 
the time for determining the course angle due to "down 
time” of P in the course of the “pause” and “signal ac- 
cumulation” SC. In order to lessen this shortcoming it is 
proposed that prior to the onset of standard discontinu- 
ous P reduction the platform be rotated in the direction 
of approach to the plane of the meridian by 90°. This 
procedure gives a decrease in the maximum time for de- 
termin + he meridian by approximately half because 
even in the worse case standard P reduction begins with 
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the initial angle of RA deviation from the plane of the 
meridian by 90°, not 180°. 


Analytic GC schemes are used for the purpose of 
reducing the time required for determining the course 
angle and simplifying the GC design. However, their 
realization has its special features and difficulties. 


Analytic GC schemes will be examined in Part II of the 
article. 
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Requirements on Navigational Support for 
Experiments From Lunar Polar Satellite Orbit 


9S7AI1117A Moscow GEODEZIYA I KARTOGRAFIYA 
in Russian Apr 95 No 4, pp 10-20 


[Article by Yevgeniy Pavlovich Aleksashin, doctor of 
technical sciences, professor, Yuriy Sergeyevich Timo- 
feyev, candidate of technical sciences, senior scientific 
specialist, and Aleksandr Mikhaylovich Shirenin, can- 
didate of technical sciences, senior scientific specialist, 
Moscow State Administration for Geodesy and Cartog- 
raphy; UDC 528.715:629.78] 


([FBIS Translated Text] The most rigorous requirements 
On navigational support for the operation of satellite- 
borne instrumentation are imposed by researchers who 
are engaged in on-line and postflight processing of the 
results of an orbital survey of the surface for the pur- 
pose of constructing a reference coordinate system and 
mapping, with det.rmination of the parameters of the 
Outer gravity field and study of the dynamics of the 
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translational-rotational motion of a celestial body [2, 8, 
11). As a rule, if these experiments are supported by ad- 
equately precise navigational data, no problems arise in 
their use for coordinate-time referencing and processing 
of the results of other experiments. Accordingly, in ex- 
amining the problems involved in navigational support 
of the instruments carried aboard a lunar polar satel- 
lite (LPS) we will first of all start from the require- 
ments which are imposed on the above-mentioned ex- 
periments. 


The coordinate-time referencing problem is formulated 
in the following way: in a coordinate system rigorously 
referenced to the moon it is necessary to determine the 
vectors of the positions 


Ry = (Xq, Yq. Ze) 


of points on the physical surface of the moon on the 
basis of the coordinates of their images measured in 
the coordinate system of the satellite-borne instrument 
and corresponding to two or more moments in time. 
A special case of this problem is determination at the 
moments of observation of the spherical coordinates of 
the traces of the principal axes of the instruments on the 
physical surface of the moon or a celestial body. We will 
use a dynamic selenocentric system, as it was defined in 
[3], as the coordinate system rigorously coupled to the 
moon. For computing R*, (Note: Here and in the text 
which follows an asterisk following a letter indicates 
that in the original text the letter has an overline; two 
asterisks — the letter has an overline with a caret above 
the overline; three asterisks — the letter carries a caret; 
four asterisks — the letter carries a tilde] we will use an 
algorithm realizing direct orbital intersection based on 
a search for the mean-weighted position of a physical 
point on the segment of the minimum distance of the 
intersecting projecting rays which correspond to the 
images of this point obtained at two moments in time. 
This algorithm is described by the formulas: 


Ras = Raa — Re Bo M75 Bie Mites: 
p) apt: Py ae et YO 
nop - pee GTA Re 2 PrRas: 
b, = Welepyepys b= W/elep,-epiak 
aR. = 0.5(k yp, +hp,). 


(1) 
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T 


= [Xsiery Yaqbeiyp Za 141] 


are the position vectors of the LPS at the times t, t,,, in 
the system rigorously referenced to the moon; 


77 = 


= [xp Y» Sh Piet @ Mew owl) 


are the position vectors of the image of the very same 
physical point, stipulated in the coordinate system of a 
videoinstrument with the focal length f at the times t,, 
bo 


tn hat M= 
* (hyp Ay Ay) Mis = (Auger divie1) Aye) 


are the matrices for transformation from the instrument 
system to the system rigorously referenced to the moon, 
stipulated at the times t,, t,,, and being functions of the 
Orientation parameters 


EF Udy Ah 
Dn - (Aas): dat): Aue) | 


It follows from (1) that the vectors R*, of the sought- 
for point is the vector function 


R, = Ally) 
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of the 17-dimensional vector of the P* parameters 


pP* = (x, Yn» x, Y, Zz, | 4, Ay, Xe: Now 


Xe Yiew 2+ Astii\s Aiitiys Ayiesys fi. 


Thus, the components of the P* vector (2) must be ac- 
cessible at the time of the experiment for the coordinate- 
time referencing of the results of on-board iments 
from a LPS orbit. They are obtained in the form of eval- 
uations together with the covariation matrices of errors 
from the results of making and processing navigational 
measurements containing random errors. The compo- 
nents of the subvectors r* and the covariation matrices 
of their errors K, are computed on the basis of measure- 
ments of the coordinates x, y of the images of physical 
points in the instrument system. The components of the 
subvectors R*, and covariation matrices K,, of their 
errors are obtained on the basis of trajectory measure- 
ments, components of the subvectors A* and covariation 
matrices of their errors K, on the basis of measurements 
of the LPS orientation parameters [1]. 


Assume that P** is the vector of evaluations of the 
parameters and that K, is the covariation matrix of errors 
of the P** vector. The coordinate-time referencing 
problem essentially involves computation of evaluations 
of the position vectors R**, of physical points and the 
covariation matrices K**,, of their errors. We have 


R, ™ FAP); R,, = GiR,G,. (3) 
where 


G, = (0,F 5) 


is a gradient 17 x 3 matrix of the vector-function F*, 
for the P parameters; 


7, =(a/ax,...8/a/) 


is a 1 7-dimensional vectorial differentiation operator for 
the P* components. 
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The components of the sought-for vectors R**, and 
the covariation matrices K*,, by virtue of (1) and (3) 
are essentially dependent on the programs, makeup and 
accuracy of the navigational measurements used for 
determining the components of the vector of the P* 
parameters, as well as on the processing algorithms 
defining the rules for computing the P** evaluations 
from the vector sample of measurements. 


The principal objective of this study was therefore 
the construction of mathematical models of the pro- 
cesses transpiring on the real flight of a LPS and di- 
rectly related to the collection of data necessary for the 
coordinate-time referencing of the results of operation of 
the onboard instruments. The constructed models were 
used in numerical computer experiments for determin- 
ing the a priori accuracy characteristics of the refer- 
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encing results. Some recommendations have been for- 
mulated on planning the LPS orbit, selecting the type, 
makeup and programs for navigational measurements 
for ensuring the necessary accuracy in the referencing 
of experiments. 


We will proceed on the assumption that the LPS car- 
ries instrumentation which can be used in navigational 
measurements: a frequency generator which fixes the 
time scale aboard the satellite and a device for regis- 
tering the moments of events transpiring on this scale; 
devices ensuring the conducting of surface radiotechni- 
cal measurements for determining the LPS orbit; gyro- 
scopic system for working out and maintaining the stip- 
ulated LPS orientation modes; topographic and sidereal 
videoinstruments; altimeter. 

















Table 1 
Experiment H, km T revolutions/day IDL, ° b,° P(%) 
I 100 15sgm 12 11 1.45 0.3(21) 
ll 200 2 08 12 1.2 2.97 0.7(23) 
Il 300 218 10 1.3 4.45 0.5(12) 























With respect to performance of the scientific missions 
of the LPS flight and supporting experiments with the 
required navigational information the parameters of the 
working orbit from which these experiments will be 
made are not a matter of indifference. In particular, 
the conditions of solar illumination of the surveyed 
surface and the geometry of its coverage by the scanning 
swaths of the topographic camera are of considerable 
importance for the problems which we examined. The 
LPS orbit must be planned in order to meet these 
or other illumination and coverage requirements. The 
usual set of parameters for the planning of orbits is: 
a — semiaxis major, e — eccentricity, i — angle of 
inclination of the orbital plane to the selenoequator, w 
— distance of the pericenter from the orbit ascending 
node, a, — selenocentric right ascension of ascending 
node, r — time of transit through the pericenter. In 
surface surveys preference is given to images of an 
approximately identical scale, which predetermines the 
choice of an orbit with an eccentricity equal to zero, 
that is, a circular orbit. Since all points in such an orbit 
are located at identical distances from the surface, it is 
unimportant which of them, instead of the pericenter, is 
taken as the reference point. Thus the parameters w and 
7 become free. Naturally, for maintaining the status of a 


polar orbit the angle i must be close to 90°. It therefore 
follows that when planning a LPS orbit it is possible to 
vary only two parameters: a, and a. 


Assume that a,, 4, are the selenocentric right ascension 
and declination of the sun at the time t. Under the term 
“illumination conditions” of some point on the lunar 
surface at the time t we wil] understand the angle of 
solar elevation 7,, above the local horizon. If a,, 4, are 
the selenocentric right ascension and declination of a 
given point on the lunar surface, then 


Yrc = arcsin(sin 6, sin 6. + eos 


+ cos 4, cos 6. cos (a, — a@,)). (4) 


It follows from (4) that for a supulated moment in ume 
a,, 4, and a stipulated declination of the point 4, the 
angle +,, is dependent only on the difference a, - a,. If 
the point belongs to the trace of the ascending branch 
of an orbit with i = 90° on the selenocentric celestial 

. then a, = a,, that is, the illumination in the 
circle of declinations 4, in this case 1s dependent only on 
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the difference Aa, , of the selenocentric right ascensions 
of the node a, and sun a,. By stipulating the maximum 
admissible 7,, in the selected circle of declinations 4, 
from (4), we obtain the values 


dag, dag 


corresponding to the “morning” and “evening” angles of 
elevation + ,.. 


Now from the relation Aa, = a, - a, we determine 
@,, ensuring the stipulated illumination conditions. The 
realization of the required a, is accomplished by proper 
choice of the moment of transfer of the LPS from a 
standby orbit to a working orbit. Due to the orbital 
motions of the Earth relative to the sun and moon a, 
increases as time passes, whereas a, for i = 90° remains 
unchanged. In order to ensure the required illumination 


nv 


a4 


nboumens (1) 
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conditions as long as possible it is necessary that the 
difference Aa,, correspond to the “morning” case, 
that is, be positive. In addition, it makes sense to 
change somewhat the LPS orbital angle of inclination 
i somewhat in order that due to precession the orbital 
node be displaced away from the sun. Node precession 
is determined using the formula [4] 


aan(Ry/e) Vs 


429 = x1-e')' “ 


cos i(t-4), (5) 
n= M 


where C,,,, = -0.206 x 10° is the zonal harmonic in 
expansion of the force function of lunar attraction; R,, 
= 1,738 km is the mean radius of the moon; p = 4,902.90 
km? /s* is the lunar gravitational parameter. 








6, 
"% ' 6UMOX 





nt aumor (1) 


Fig. 1. Geometry of coverage of lunar equator by satellite transits as function of spacecraft altitude, 
scanning swaths and end laps. 


1. Satellite transit 


It follows from (5) that Ala, > 0 when cos i < 0, that is, 
when 90° < i < 180°. In order to retain the possibility 
of observing the pole from orbit it is admissible that 
the camera axis be deflected from the local vertical in a 
range 5°, that is, the status of the polar satellite will be 
conserved if i = 95°. 


The LPS orbit semiaxis major a entering into (5) is 
the most significant parameter determining the geom- 
etry of coverage of the lunar surface by topographic 
camera scanning swaths. The change AL, of the se: 
lenographic longitude of the LPS orbit ascending node 
during a satellite transit (transit interval) is governed by 
two effects: by the precession 4 during the LPS 
revolution period T = t -t, and the displacement of the 
node during the time T due to lunar rotation, that is, 


Alo i dag -e,7, (6) 


where w,, = 0.2661699489 x 10° p/s is the angular 
velocity of lunar rotation; 


| 7 = 2aVe/Vp 


is the period of LPS revolution. 


Using (5) and (6) and varying the a parameter, it is 
possible to obtain different AL, (a). However, for the 
considered problems there is interest only in those 
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AL,(a) which are found in a definite relation to the 
width of the topographic videocamera scanning swath 


b, = 2[arein| "5 sin arcs 39 - 
- arcig 7]. <i (7) 


where | is the length of the frame edge; f is the 
videocamera focal length; H is LPS flight altitude. 


Figure | shows the geometry of coverage of the 
selenoequator by satellite transits as a function of the 
altitude of the circular orbit H, = a - R,. For a 
Caiucta with the parameters | = 180 mm, f = 400 mm 
the figure also shows the scanning swaths (7) and their 
end laps. Table 1 gives the periods of revolution T, the 
number of rotations per terrestrial day and interval of 
satellite transit AL,. If the survey scale is limited to 
not less than 1:800,000 and a minimum side lap P, as 
follows from Fig. 1 and Table | the optimum orbit (in 
the sense of coverage geometry) will be an orbit with 
an altitude H__., = 300 km, a = 2,038 km. 


In Fig. 2, in a projection onto a selenocentric celestial 
sphere with the radius Ry, are plotted a spot of 
illumination of the moon by the sun with isolines of 
angles of equal solar elevation +, and 10-day traces 
of the LPS orbit showing the relative change in the 
positions of the orbital traces in the solar spot as a 
function of time. The initial position of the ascending 
orbital node relative to the sun, that is, the positive 
(“morning”) difference tygDla,, at the initial moment 
was selected under the condition that at the point of 
the trace 6, = 40° the solar angle of elevation is 4, 
= 30°. In Fig. 2 it is possible to trace the change in 
illumination conditions both along the flight trajectory 
and as a function of the change in the difference of the 
selenocentric right ascensions of the ascending node of 
the LPS and the sun. The convergence of the node and 
sun after 10 days is about 9°, that is, Ala, during each 
10 days of flight changes by about 9°. 


It also is of interest to ascertain the required number 
of days for continuous mapping of the moon from the 
selected orbit by a topographic videocamera (f = 400 
m, | = 180 mm). Figure 1 and Table | show that when 
working with a side lap of 23 percent the program for 
camera activation on the ascending revolution must be 
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structured in such a way thal a survey is made after 
two transits on the third, that is, on the n+3, n+6, etc. 
transits. The displacement of the survey transits over 
the physical surface of the moon due to its rotation and 
precession of the orbital node using formulas (5) and 
(6) even in this case is AL, = -3.9°. 


The time between survey transits is 0.2875* (24,840 = 
6.9"). Proceeding on the basis of AL,,, it is necessary to 
have 92 such transits in order to cover the selenoequator 
completely by scanning swaths with an overlap 23 
percent. Hence the required number of days for a 
continuous lunar survey is 26.45*. During this time the 
difference Aa,,, changes by -23.8°. In the middle of the 
survey period it is desirable that Ala, = 0°. Then it 
is necessary to make the initial Ala, equal to +11.9°. 
Figure 2 shows the initial and final positions of the 
survey transits in the solar spot. 
























































Fig. 2. Relative change in positions of traces of LPS 
orbit in solar spot on the moon as a function of 
time. 


1. Days 


Thus, a LPS orbit with the following parameters was 
selected for numerical modeling: — circular orbit (e = 0, 
w = 0, a = 2,038 km, H = 300 km); — inclination angle 
i = 95°; — initial difference Ala, of the selenocentric 
right ascensions of the ascending orbital node and the 
sun +11.9°, that is a, = a, + 11.9°; — the initial time 7 
of transit of the node is equal to the time during which 
a, is computed. 


For a comparative assessment of different sources of 
Navigational information when determining the vector 
P* a computer was used in making experiments simu- 
lating the real processes of navigational determinations 
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of the orbit and orientation of the LPS and in computing 
the selenocentric coordinate vectors R*, of the observed 
points. The referencing efficiency criterion used was the 
Euclidean norms of the traces of the covariation matri- 
ces K,, of the errors of the vectors R*, 


wa 
Ser = (Ox + Fy +9) 


where a,, 7,, 7, are the dispersions of the errors of the 
components of the vectors R*, standing on the diagonal 
of the matrix K,,, determined from expression (3). 


A LPS orbit, selected by the method described above, 
was adopted as the reference orbit for all experiments, 
whereas the orbital orientation mode in which the +z 
axis of the topographic camera was directed along the 
LPS radius-vector, the +x axis was directed along the 
transversal in the direction of motion and the +y axis 
was directed along the vector of the kinetic moment 
was used as the nominal mode of LPS orientation. 
Coaxially with the topographic camera system there also 
is the coordinate system of the sidereal camera, but their 
optical axes are directed in opposite directions (-z and 


+z respectively). 
A system of first-order nonlinear differential equations 


[2] was used for predicting the orbital motion of the 
LPS relative to the physical surface of the moon. 


R, = ¥,; nex 


, - a= ne Bnac—2@a V,) - 
5 


where g*,. , , is the vector of the total perturbations 
from the acentrality of the lunar gravitational field 
and the attraction of the Earth and sun; w*,, is the 
vector of angular velocity of the moon; €*, is the 
random vector of unsimulated accelerations with the 
mathematical expectation E[gy*, . . and the stipulated 
covariation matrix 


Ky = EE EY} 
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Since the components of the vectors R*,, V*, are 
determined on the basis of trajectory measurements, 
they are known with errors. Accordingly, it also is 
necessary to predict the covariation matrix of errors of 
the vectors R*, along the trajectory. We introduce the 
LPS vector 


Tr -_ 


T 
a: =(R,.Vs5) 


and its covariation matrix of errors K,,. Then, in 
accordance with the linear theory of stochastic processes 
[5], the covariation matrix K es satisfies the Riccati linear 
differential equation 


K = OKs + (GK) + Kew (9) 
é 


where G., is the gradient matrix of the right-hand sides 
of equations (8) for the q*, components. This equation 
also was used in the simulation for predicting K,, along 
the orbit. 





tatu aket 
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Fig. 3. Trace of optical axis of Zond-8 spacecraft 
camera on celestial sphere in inertial orientation 
mode during lunar survey. 
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The LPS phase vector q*, and the covariation matrix 
K,, at the time of measurements are evaluated using a 
Kalman filter [5S] employing the following algorithm 


(GK.,0,+K)W = 1; U = K-G,W; 
a «7 +ul(d-d} - 6,(% - a) |: 


Ky = (1-U0}) Ke (1-UGl)* + U Rv", 


where G, is the gradient matrix of the measured 
functions of the components of the phase vector q*; 
ae " is the covariation matrix of errors of the vector 

, predicted on the basis of equations (9) at the 
nd of measurements; K****, is the covariation matrix 
of the measurement errors; q**, is the LPS phase 
vector predicted using equations (8) at the time of 
measurements; q**,* is the evaluation of the LPS phase 
vector obtained afler processing of the measurements; 
d** is the vector sample of measurements; K*** .* is 
the covariation matrix of errors of the vector q**.’, 
obtained after processing of measurements; d* is the 
computed value of the vector sample of measurements; 
q**,.* is the value of the vector q**,* when perf: ming 
an iteration with the number v; I is a unit matrix. 


Two approaches were used in simulating the process of 
determining LPS orientation. 


In the first it was assumed that the LPS is in an orbital 
orientation mode, supported by gyroscopic systems 
within the limits of the accuracy level admissible for 
them. An analysis of operation of actual orientation 
systems in such a mode shows that the deviations of 
the spacecraft axes from a stipulated position within the 
limits of the threshold of response of the gyroscopes are 


unpredictable. 


The space vehicle control system maintains a stipulated 
Orientation in a pulsed mode and the motion of the 

ical axis of the camera resembles Brownian motion 
Gig. 3) (2, 7, 8). Accordingly, in this study when 
using the first approach the orbital orientation mode is 
simulated by the uniform rotation of the LPS relative to 
the y-axis in such a way that the -z axis traces the lunar 
center of mass. This rotation is perturbed by random (in 
the form of white noise) moments relative to each axis. 


In the second approach the model of LPS rotation is 
not stipulated and the actual orientation at the time 
of the experiment is determined from siderea! camera 
data. The measurement information is formed from star 
images determined in the camera coordinate system, x. 
y and the corresponding equatorial coordinates a****, 
é**** of stars selected from the onboard computer 
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catalogue. In this case the following iteration algorithm 
[1] was used for determining the LPS orientation 
parameters A* and covariation matrices K, of their 
errors: 


S = BD+Ac, 6 = _(F+ Al + BAd); (1) 
SL =}, 


I = cL+Al; 3d = DE +Ad; 


where A** are the a priori orientation parameters; A*,,. 
q+ iS the vector sample of the measurements; 


T° = (x, » Xt, 4, Zh 


x**** y**** are the coordinates of the image of the 
i-th star, measured in the instrument coordinate system; 


are the components of the unit vector of the star from 
the catalogue corresponding to it; d* is the vector of the 
computed measurement values; K****, is the a prior 
covariation matrix of errors A**; A is a gradient matrix 
consisting of the partial derivatives of the limitations F* 
for the A* components; K****, is the covariation matrix 
of measurement errors; B is a gradient matrix consisting 
of the partial derivatives of the limitations F* for the d* 
components; F* is a vector function of the limitations. 


In this problem for the i-th star the limitations are 
written in the form 


F = &*- KM,, (I)? = inst in 
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K, is a proportionality factor; A**,, K**,,, is the 
sought-for vector of the orientation parameters and the 
covariation matrix of its errors. 


The evaluations A**., of the LPS orientation parameters 
computed using algorithm (11) describe its orientation 
in the star catalogue system, that is, in the geoequatorial 
system of the means of the equator and equinox 
of the epoch J 2000.0. However, for their use in 
the intersection formulas (1) they must describe LPS 
Orientation in a system rigorously referenced to the 
moon. The parameters A**, and the matrix K***,_, are 
scaled using the algorithm: 

the A** are used at the measurement time t 
in forming the matrix M,,,, ,, (t, A**,, of the “instrument 
- inertial” conversion; 

then the following are computed at the time t: M(T,, 
t) is the matrix for precession from T, = J 2000.0 to 
t and M,(t) is the “equator-ecliptic’ conversion matrix 
(3); M,,(t) is the "ecliptic-system rigorously coupled to 
the moon” conversion matrix [10]; 

the “instrument-rigorously coupled to the moon” con- 
version matrix is computed at the moment in time t 


M() = MylQMQMy(Ty WMygs (12) 


the sought-for parameters A*, are computed using the 
te 


the sought-for parameters \* are computed using the 
M(t) matrix. 


In a general case the process described above can be 
represented in the form 


f = Al). (13) 
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Then in a linear approximation 


(14) 


Ry = BKipGip 


where G,, is a gradient matrix containing the partial 
derivatives of the vector function F*, for the compo- 
nents A* and the parameters of physical libration of the 
moon 7, 7, p. This matrix is computed using the central 
finite-difference approximations of each derivative en- 
tering into it; K***, is the covariation matrix of errors 
of the parameters A** and the parameters 7, 7, p, which 
is stipulated in the form 


Ri. 0 | 
Ri, = |—wetn “|. 
Kom 








The submatrix K, , , is determined on the basis of 
studies of the accuracy of the libration models used at 
the present time (6). 


In computer simulation of the process of coordinate- 
time referencing of the results of studies from an 
LPS orbit the makeup and program of navigational 
measurements were selected on the basis of the real 
conditions for their application and were varied from 
experiment to experiment in the following way. 


Experiment R + G. The orbit is determined only 
from surface radiotechnical (R) measurements once at 
12". Orientation data are received by telemetry from 
gyroscopic systems at any moment in time during 
operation of the onboard scientific instruments. 


Experiment R + S. The orbit is determined the same 
as in the preceding case. Orientation by the stars (S) is 
Carried out at any moment in time. 


Experiment R + T + S. The orbit is determined from 
surface radiotechnical measurements onc: at 12" and 
from photogrammetric (T) measurements of onboard 
videoimages of the lunar surface on each survey transit, 
that is, once at 6°. The orientation is determined from 
the stars for any moment in time. 
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Fig. 4. Graph of " prediction-evaluation" process 
for numerical experiment R + T + S when 
predicting LPS orbit in interval 2 days. 


The software which we developed for computer-aided 
simulation of the navigational referencing process al- 
lows variation of virtually any parameters on which this 
process is functionally dependent. However, in this case 
the referencing accuracy which can be attained when 
Operating existing equipment is of interest. This also 
determined the makeup of the parameters adopted as 
the initial parameters when making the experiments: 


topographic videocamera: focal length f = 400 mm, 
length of frame edge, | = 180 mm; star camera: f, = 
70 mm, | = 10 mm; 


accuracy in measurements of the coordinates x, y in the 
videoimage field 7, = 0.2 x 10' mm; 


the surface trajectory measurements were stipulated in 
the form of LPS phase vectors q**, in a selenocentric 
dynamic system as the results of processing of sessions 
of direct measurements of range and radial velocity. The 
accuracy q**, was stipulated on the basis of data from 
the Zond-8 spacecraft [2]: o,,. = 500 mm, o,, = 0.1 m/ 
s; 


the gyroscopic system formed the orbital orientation, 
the accuracy in maintaining which was determined in 
conformity to the data in Fig. 3: 7,. = 1,800"; 


orbital photogrammetric measurements. The trajectory 
was stipulated in the form of a series of coordinates 
of the LPS centers of mass at the times of survey in a 
dynamic selenocentric system. They were obtained as 
a result of construction of a photogrammetric network 
along the survey path, as described in (2). Their 
accuracy was determined as o,. , = 150 m. 


The simulation process involves repeated repetition 
of the following steps in the time interval of active 


existence of the LPS (26°). 
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Fig. S. Errors in angular orientation of star sensor 
for narrow- and wide-angle cameras. 


The Cauchy problem is solved using the initial data 
y ” 
@s (Q, Ky (4) 


for the differential equations (8) and (9) for the purpose 
of predicting q**,, K***., at the moment t,, , of the 
next measurement. As a result the predicted evaluations 


-*» A 
Qs (howds Kolber 


are computed. 


Depending on the type of trajectory measurements (sur- 
face radiotechnical or orbital photogrammetric) at the 
time t,,, they are processed by a Kalman filter iteration 
algorithm (10). As a result the refined evaluations 


C3 (beds R alte) 


are computed. These data occupy the corresponding 
places in the vector P** (2) and the covariation matrix 
K***, (3). Figure 4 shows a graph of the “prediction- 
evaluation” process in a steady state for the case R + 
T + S for determining the LPS orbit in an interval of 2 
days. In a prediction interval 6 hours the error in LPS 
position attains |,300 m, but after determining the orbit 
at the beginning of the survey half-transit it is about 80 
m. Figure 4 shows that within the limits of a working 
half-transit (1 hour) the prediction error is about 217 m. 
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Simultaneously at the moment in time |, an evaluation 
is made of the orientation parameters A*. This evalu- 
ation is generated by a random process of a stipulaied 
intensity if use is made of data from gyroscopic systems 
or is computed using data from x, y measurements of 


star images. 


In any case the evaluation algorithm (11)-(14) is used 
for this purpose. As a result, the evaluations 


a 


are computed at the time t,,,. These are entered into 
the vector P** and the matrix K***,. Figure 5 shows 
the field of view of the star videoinstrument typical for 
this case and the errors in determining LPS orientation 
obtained from measurements of the coordinates of star 
images. 


Such images for each moment of observation are formed 
using data from the star computer catalogue and the 
parameters of the computer videoinstrument. The size 
of the field of view changes, depending on the length 
of the frame edge. Figure 5 shows narrow- and wide- 
angle fields. In both cases the greatest error is for the 
angle ¢ of turning of the image about the optical axis. 
It is virtually an order of magnitude greater than for the 
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other angles and is decisive when assessing the accuracy 
in LPS orientation. It also is interesting to note that 
the errors in the angles y and @ are in actuality not 
dependent on the size of the frame within the limits of 
the stipulated variations. 


Uniformly positioned points are selected in the zone 
of overlap of topographic camera images registered at 
the times t, t,,. The coordinate vectors R*, of the 
observed points and the covariation matrices K,, of their 
errors are computed for each of these points using the 
algorithm (1) and employing evaluations of the position 


Bra), Reclts Rte)» Kislted 
and the orientation 


Ti) Rik), Lear)» Ritter. 


positioned in the vector P**,*(t,, |, and the matrix 
K***.*. The norms of the traces of these matrices 7,, 
characterize the accuracy in referencing of the results of 
operation of the onboard instruments at the times ,, {,,, 
to the system rigorously coupled to the moon. 
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Table 2 
The Error im determining podtion, ortentation and referencing in experiment 
boars 
R+G; R +S; R+T+S 
ons ™ rs Cop ™ Ons ™ o «- Opp ™ Cas ™ 15: Cap ™ 

6 85 99 266 

12 150 Ms 1240 150 211 424 79 47 242 

»” 82 54 254 

% 167 1802 2195 167 274 612 85 116 328 
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The, Error in determining position, orientation and referencing in experiment 

=a R+G; R +5; R+T+S 

Cay ™ ors opp ™ Try = P15: opp = Crs ™ ors: Cap m 

54 77 101 255 

60 150 1804 733 150 368 460 85 185 282 

66 78 121 297 

72 156 1802 1786 156 294 533 82 106 302 

Mean 1398 492 273 


























In Table 2 we give the results of an evaluation of the 
accuracy in referencing obtained using numerical exper- 
iments on a computer. The greatest referencing accu- 
racy a = 273 m can be attained if for determining 
the orbit joint use is made of data from surface 
trajectory measurements and onboard surface videoim- 
ages, but for determining the actual LPS orientation — 
data from onboard videoinstruments containing images 
of sectors of the star sky (experiment R + T + S). In 
this case the referencing error is almost five times less 
than in determinations by the standard method (R + G). 


There is stil! another possibility for increasing the accu- 
racy of navigational referencing when using instruments 
with the adopted parameters. As was noted above, the 
worst accuracy in determining orientation is for the an- 
gle ¢ of turning of the star camera relative to the optical 
axis. If there are two star cameras mounted aboard the 
LPS at an angle about 90° the accuracy in determining 
the actual orientation increases substantially, which also 
results in an increase in referencing accuracy. Appropri- 
ate computations were made for confirming this point. 
They gave the following results: 7,,_. = 81 m, o,, 
= 13", Opneen = 63 m, that is, the accuracy increases by 
a factor of approximately 4 in comparison with the R + 
T + S method and by a factor of almost 20 relative to 
the standard R + G method. 


The cited results show that for navigational referencing 
of experiments from a LPS orbit with the required accu- 
racy it is necessary that it carry aboard instruments mak- 
ing it possible to refine the orbit from videoimages of 
the lunar surface and to determine the actual orientation 
from videoimages of sectors of the star sky. Depending 
on the purposes of the studies, in the operation of a LPS 
it is necessary to select the elements of the initial or- 
bit and the required accuracy in navigational referenc- 
ing. Precisely these data are the point of departure in 
a search for the optimum operating programs, makeup 
and parameters of onboard instrumentation which will 
ensure navigational referencing with the stipulated ac- 
curacy. In this case algorithms and modeling program 


systems similar to those described in this study may be 
useful. 
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Modeling a Pilasma-Tuned Satellite Antenna 
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received 31 Dec 93; first paragraph is IZVESTIYA 
VYSSHIKH UCHEBNYKH ZAVEDENIY. RA- 
DIOFIZIKA abstract) 


{[FBIS Translated Text} Experimental and theoretical 
studies were conducted to examine the effect that a 
radio-frequency [RF] discharge excited in rarefied gas 
by the field of a radiating antenna has on the antenna’s 
input and radiating characteristics. The conditions of the 
experiment were such tha! the inhomogeneity plasma 
was self-localized near the end of a dipole antenna 
and the parameters of plasma formation could be 
altered by varying the pressure of the gas, frequency 
of the signal, and power fed to the antenna. Excitation 
of an RF discharge results in modification of the 
antenna’s near field and, consequently, in alteration of 
its electrodynamic characteristics. 


it is a known fact that using a large-signal operating 
mode on satellites and descent vehicles may result 
in excitation of an RF discharge in the rarefied gas 
surrounding the antenna. Spontaneous development of a 
discharge (a plasma inhomogeneity) generally results in 
deterioration of the radiator’s characteristics, detuning 
of the feeder section, and (consequently) partial or 
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even complete loss of radio communications. As recent 
studies have shown, however, a controlled RF discharge 
may be self-localized in the vicinity of the end of 
the dipole. The antenna’s near field changes during 
the process of formation of the plasma inhomogeneity, 
which leads to redistribution of the antenna current and 
modification of the radiator’s input impedance. This 
article is devoted to studying the possibility of using 
an end plasma load to control the characteristics of 
a thin antenna (a << |, = A, where a and /, are the 
radius and length of the dipole and A is the vacuum 
wavelength). The results of an experimental study and 
numerical modeling of such an antenna are presented in 
the article. 


Figure 1 illustrates the laboratory model on which 
the experimental study was conducted. An asymmetric 
dipole located over a metallic shield was studied, and 
the end of the dipole was placed in an evacuatable 
tank. The antenna was excited by an RF generator that 
created and maintained a plasma inhomogeneity around 
the end of the dipole. All of the antenna’s parameters 
were measured by recording the distribution of current 
along the antenna. The mean square value of the antenna 
current was recorded by using an azimuthal magnetic 
field sensor. The sensor was based on an oscillatory 
circuit that had a microminiature incandescent-filament 
lamp connected to it. The light signal from the sensor, 
which was a function of the value of the antenna current, 
was fed along a lightguide to the photoreceiver unit and 
recorded by a voltmeter or oscillograph. The effect of 
the end plasma load was studied with a constant voltage 
at the RF generator’s input. Figures | and 2 present the 
results of the experiment. 


In the case of a short antenna (/, < /4 [Figure 1)}), 
after the end plasma load had been formed, the antenna 
current increased as a result of resonance tuning. The 
maximum antenna current was shifted to the antenna’s 
free end, which in turn increased its radiation resistance. 
Both the resonance tuning of the dipole and the increase 
in radiation resonance led to an increase in the power 
emitted by the antenna. In the case of a long antenna 
(i, 2 A [Figure 2]), both standing and traveling wave 
conditions occurred. In both cases, an increase in the 
power radiated by the antenna could be observed. 
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Figure 1. Diagram of the experiment: 1, radio-frequency generator; 2, metal screen; 3 glam docharg bab 
4, antenna; 5, current sensor; 6, capacitance divider; 7, plasma inhomogeneity; 8, oscillograph; 9, 
distribution of current along the antenna (1, without a plasma load; II, with a plasma load; f = 50 MHz; p = 


3 torr; —, 


All of the principal distinctions of the operation of an 
antenna with a plasma load may be explained within 
the framework of the method of equivalent circuits. In 
accordance with the said method, an antenna with a 


ee ee ee 


arm; Z,, = 60fin(A/xa) - 0.58] is the wave impedance 
of the line; A = 8 - ia is the wave number; 


B = 2nid; a = RJZ,; and R, is the line’s per-unit- 
length resistance as selected by proceeding from the 
requirement of equality of the radiation losses and losses 
in the line. The distribution of current in a line with an 
equivalent length and complex load Z, at its end is given 
by the expression 


ee experiment); 10, lightguide. 


31/2w -1 


UT bal fot ~ spe 7 tin - 2 


where z is the longitudinal coordinate and /, is a 
constant. Figures 1 and 2 present the results obtained 
within the framework of this method. The values of Z, 
were selected by proceeding from the requirement of the 
best consistency with the data obtained in experiments. 
Traveling wave conditions occur in the case where the 
impedance of the plasma load equals the antenna’s 
wave impedance (Figure 2, curve 2). Standing wave 
conditions occur when the said impedances are unequal 
(Figure 3, curve 3). In both cases, the input impedance 
and radiation resistance of the antenna are modified. 
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Figure 2. Distribution of current along the antenna (/, = 1.5 m; f = 180 MHz) outside the bulb. 1, without a 
3, 


A stricter analysis of the operation of a dipole antenna 
plasma load was performed within the framework 

of the integral equation method.’ An asymmetric thin 
end plasma load and is located over 

an infinite metallic plate and excited by the aperture of 


mode! of a symmetric dipole located in a free space and 
excited by a ring of magnetic current with a density of 


with a plasma load; 2, p = 0.5 torr; 3, p = 0.1 torr; ..., experiment; —, theory. 


The requirement that the total electrical field on the 
antenna’s axis equal zero may be used to obtain an 
integral equation relative to the antenna current. The 
total electrical field on the anvenna's axis is created by 
the following sources: the annular magnetic current, the 
antenna current, and the current of the plasma space. 
The effect of a plasma load on the total electrical field 
may be viewed as the field of radiating polarization 
currents j = iuXee, - ¢,)E, where D is the total electrical 


field, « = 
en1-5 tL, 
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Figure 3. Model of the antennas: a, monopole excited by a coaxial antenna; b, dipole excited by an annular 
magnetic current. 1, metal plane; 2, coaxial line; 3, antenna arm; 4, plasma space; 5, annular magnetic 
current. 


frequency, w is the frequency of 
the signal, and v is the frequency of collisions of 
electrons with neutrals. In the case where the plasma 
layer is rather thin, the electrical field has only a 
radial component that may be expressed through the 
distribution of current along the antenna conductor: 


£ 
= 
e 
i 


i-3¢ a 
4= e ds 2xg' 


where / is the antenna current and | _,z are cylindrical 
coordinates. The condition for the electrical field on the 
antenna axis makes it possible to write the following 


| 


| Bula) + Bos) + Bods) =0, (1) 


where E,(z), E,(z), and E,(z) are the z components 
of the electrical fields excited by the antenna current, 
plasma current, and annular magnetic field, respectively. 
It is known that the electrical field may be represented 
as 


° . i 1 
B= ——(VdivA + k*A) - 
iwtette ) + ) . “ead + 
where & is the wave number, 


. te -. & ate 
S & 3 dd I bere 


are vector potentials, j is the current source, and r is 
the distance between the point of the source and the 
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observation point. By using the latter two expressions, 
we can obtain the following: 2. Fa(4a) + B5(20) + Buslse) = 0, 
s ° : 
1 
E,,(z “<i. | where E_* and E_* are components of the total electrical 
rite A(2\(2r{(1 + thr) field linked to the n-th harmonic (2) and z, represents 








where a, and b, are the inner and outer radii of 
the plasma space and z, and z, are the longitudinal 
coordinates of the plasma’s boundaries. Condition | thus 
yields an integral equation relative to the known current 
I. The method of moments, which has been described 
elsewhere,‘ is used to obtain a numerical solution of the 
integral equation. 

In accordance with that method, the antenna current is 
represented in the form 


I= ¥ hen {= Ny} (2) 


. 


To find J, we replace the strict condition (1) with the 
following less strict condition: 


the coordinates of a discrete set of points spaced evenly 
along the antenna axis. By taking this condition into 
consideration, we can obtain a system of linear algebraic 
equations equivalent to the integral equation system: 


| 2h +. + Zinke = V; 
7 ‘Snh + ...+ Sake = Vy | 








Sash; +... + 3], =. 


where 
Bam) + Be (te 
Fun ars wate rae) Ven = ~Bi( Sa). 


This system was solved relative to J. The rapid con- 
vergence of the solution and small amount of required 
computer time were achieved by replacing series (2) 
when calculating E,, with the series 


r = ri, nz’). 


The values [’, were expressed through /. by proceeding 
from the requirement of mean-square consistency of / 
and [ at those segments of the antenna covered with 
plasma. Thanks to the condition |, << |, (/, being the 
length of the plasma loz” + may be much lower than 
N, which speeds up the computation significantly. 


Figure 4 presents selected results of the numerical 
modeling. The formation of a plasma load close to 
the end of the dipole results in an evening out of the 
distribution of current in the dipole and its resonance 
tuning. 
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Current distribution along the antenna 


2 are the real and imaginary parts 
parts of the antenna current with a 


er = 0.15, a/X = 0.007) obtained by numerical 


the antenna current without a plasma load, and 3 
load (w/w = 1.5, ww = 0.1, 


and 
with the length and radius of the plasma load equal to 0.02 and 0.014,, respectively, and with U = 1 V). 


The results of the experiments and theoretical exami- 
that is located near an antenna’s end and formed by the 
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{[FBIS Translated Text] The results of research car- 
ried out at the NPO [Scientific Production Associa- 
tion] Molniya are presented. Different concepts of mul- 
tiply reusable space transport systems are examined 
and the requirements on such systems are formulated. 
In particular, three types of winged orbital aircraft 
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with liquid-propellant rocket engines are discussed: an 
aerial-launched orbital aircraft with a jettisonable ex- 
ternal fuel tank (MAKS concept), an aerial-launched or- 
bital aircraft with internal tanks (HOTOL-Interim con- 
cept) and a single-stage orbital aircraft with a vertical 
takeoff and horizontal landing (SSTO concept). The de- 
pendence of payload mass on orbital altitude and incli- 
nation is determined for the three studied orbital aircraft 
variants. 

The NPO Molniya has studied possible ways to develop 
future multiply reusable space transport systems [1, 
2}. Considerable efforts are being undertaken in the 
United States for studying the possibility of replacing 
the Space Shuttle at the beginning of the next century 
[3-6]. Despite some differences in the approach to this 
problem, it can be concluded that the development of 
single- and two-stage launch systems with promising 
ramjet engines with subsonic combustion and with 
supersonic combustion is a task for the remote future. 
Different variants of multiply reusable space transport 
systems using only jet engines are regarded as the most 
preferable variants of new launch systems since they 
make it possible to reduce the specific cost of launch 
of payloads into circumterrestrial orbit. It is possible to 
define two principal directions in research on multiply 
reusable space transport systems with rocket engines for 
a vertical launch from the earth and a horizontal launch 
from a subsonic mother aircraft. 


Two-stage fully reusable systems of the Shuttle-II type 
require solutions of the complex problem of guaranteed 
engine-free return of the first stage to the launch site. 
It is possible that for this reason over the course of 
the last several years the United States has concentrated 
its research primarily on the concept of a single-stage 
carrier (SSTO) [3]. There is the widely known project 
for the Delta Clipper — a conical carrier for vertical 
takeoff and vertical rocket landing — developed by 
McDonnell Douglas, constituting one of the possible 
variants of construction of a single-stage carrier (Fig. 
l,a). A winged single-stage carrier for vertical takeoff 
and horizontal landing is another direction in research 
on promising single-stage carriers [3] (Fig. 1,b). 


The multipurpose aerospace system (MAKS), based on 
a heavy subsonic An-225 aircraft (horizontal takeoff and 
horizontal landing), will ensure a new quality due to the 
mobility of the launch. The MAKS system provides for 
three variants of the second stage positioned on top of 
the An-225 fuselage [7]. Variant | includes a manned or 
unmanned orbital aircraft outfitted with rocket engines 
and a jettisonable external fuel tank (Fig. 1,c). Variant 2 
inciudes a launch unit, used a single time, atop which is 
placed a payload (space vehicle with a space tug). This 
is a cargo variant of the system, known as the MAKS- 
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T. Variant 3 is a multiply reusable orbital aircraft with 
liquid-propellant sustainer rocket engines and internal 
fuel tanks. Variants | and 2 may become the first 
Stage in developing a highly efficient aerospace launch 
system because both these variants are based on actually 
existing technologies. As a result, in the first stage 
of construction of the aerospace system the technical 
risk is reduced and the time and cost of designing 
and producing system components are reduced. Variant 
3 requires the perfecting of new technologies related 
to the problem of integration of multiply reusable 
heat shielding, structure and cryogenic insulation. This 
variant is regarded as the second stage in a multipurpose 
aerospace system. The technical risk in the second stage 
must be reduced by use of the experience accumulated 
in the first stage. This variant also is known under the 
name An-224/HOTOL-Interim (Fig. 1,d). Some aspects 
of this concept were examined in the course of an 
international investigation which was carried out by the 
Organizations British Aerospace (Great Britain), NPO 
Molniya, Central Aerot ydrodynamics Institute, Central 
Institute of Aviation Engine Building (Russia) and 
ANTK Antonov (Ukraine). 


The principal results given in this article relate to inves- 
tigations of a winged single-stage Carrier with a vertical 
takeoff and a horizontal landing (SSTO concept), vari- 
ants | and 2 (MAKS concept) and variant 3 (HOTOL- 
Interim concept). The launch characteristics of all vari- 
ants of multiply reusable space transport systems are 
improved in the case of use of the three-component 
liquid-propellant rocket engine RD-701, developed at 
the NPO Energomash (Russia). It provides for the con- 
sumption of three fuel components in the first mode 
(kerosene, liquid hydrogen and liquid oxygen) and two 
components in the second mode (liquid hydrogen and 
liquid oxygen). 


We will examine the requirements on multiply reusable 
space systems. One of the early systems was developed 
in 1965 by the Mikoyan Design Bureau. It was the two- 
stage Spiral aerospace system based on a supersonic 
mother aircraft. A great number of concepts for multiply 
reusable space systems are now being studied in differ- 
ent countries: NASP, Delta Clipper (United States), HO- 
TOL (Great Britain), Zenger (Germany), MAKS (CIS), 
Ariane-5/Hermes (France), H2/Hope (Japan) and others. 
Such a diversity makes it necessary to formulate some 
fundamental requirements on multiply reusable launch 
systems. These requirements can be divided into three 
parts. 


Mandatory requirements: a substantial decrease in the 
cost of launch and return of payloads; ecologic clean- 
ness; reliability; safety; high operational characteristics; 
maximum use of experience with aerospace technology. 
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Special requirements: routine servicing of orbital sta- 
tions and other orbital objects; launch into equatorial or- 
bits for countries located in the middle latitudes; possi- 
bility of use of the orbital stage as a satellite for the solu- 
tion of various problems in orbit; realization of manned 
and unmanned variants; high landing accuracy. 


General requirements: good prospects for improvement 
and development; use of individual system components 
for other branches of the economy; assurance of capa- 
bility for handling not less than 70-80 percent of all 
future payloads. 














Fig. 1. 


Three-component liquid-propellant jet engines: Re- 
search has shown that for the launch of payloads of 
considerable mass in the case of an aerial launch from 
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existing subsonic mother aircraft it is necessary to use 
liquid hydrogen LH, and liquid oxygen LO,. A short- 
coming related to the great volume of the hydrogen 
tank can be partially offset by the three-component sus- 
tainer engine. The RD-701 engine ensures a 40-percent 
decrease in the working mass of hydrogen fuel and has 
a low weight-to-thrust ratio about 0.009. In the first 
mode the RD-701 has a specific thrust in a vacuum 
of 415 s (with a 6-percent LH, content in the total 
per-second fuel consumption and 461 s in the second 
mode (with a 14-percent LH, content) with a thrust 
reduced to the level 40 percent. Despite some decrease 
in the mean specific thrust along the trajectory, the 
RD-701 ensures a substantial gain in MAKS system 
payload mass due to a reduction in the mass of the 
external fuel tank and orbital aircraft. The gain in 
payload mass in the manned MAKS variant (with a 
takeoff mass of the second stage 275 tons) is estimated 
at 1.8 tons in comparison with the two-component 
engine of the Space Shuttle sustainer type (SSME) with 
a specific thrust in a vacuum of 455 s. 


The HOTOL-Interim variant is more sensitive to the 
use of three-component liquid-propellant engines: the 
gain in payload mass is about 3 tons in comparison 
with engines using LH,/LO, components having a 
specific thrust in a vacuum of 463 s. It must be 
noted that the "kerosene, liquid oxygen" components do 
not ensure launch of any appreciable payload quantity 
in the HOTOL-Interim variant and in MAKS system 
variant 1. Thus, the RD-701 is the optimum choice for 
aerospace systems with restricted takeoff mass. Similar 
conclusions concerning the advantages of the RD-701 
for the An-225/HOTOL-Interim system were drawn in 
[8]. 


Although the restriction on the takeoff mass in the 
concept of a single- stage carrier is not as rigorous 
as in an aerospace system, the use of three-component 
liquid-propellant engines gives a substantial effect. The 
RD-701 engine has a high thrust level during takeoff 
(first mode), which is very important for reducing the 
gravitational losses of characteristic velocity during a 
vertical launch. The second mode with a low thrust level 
makes it possible to dispense with strong throttling at 
the end of the trajectory. Thus, a high thrust during 
takeoff and a high specific thrust in the final part of 
the launch, being characteristic features of the sustainer 
engines of the first and second stages of an ordinary 
booster, are ensured in the very same three-component 
sustainer engine of a single-stage carrier. 


For a vertically launched winged single-stage carrier 
with a fixed takeoff mass of 550 tons data were obtained 
characterizing the advantages of the three-component 
liquid-propellant RD-701 rocket engine in comparison 





This report may contain copyrighted maternal. Copying and dissemination 
is prohibited without permission of the copyright owners. 

















So SPACE SCIENCE AND ENGINEERING 


with the Space Shuttle sustainer engine (SSME); the 
mass launchable into an elliptical orbit is 70,720 and 
73,460 kg; the SSTO dry mass is accordingly 50,880 
and 61,980 kg and the payload mass in circular orbit is 
10,400 and 1,940 kg. Calculations were made for the 
case of use of an extensible nozzle at altitudes greater 
than 9 km in each sustainer engine variant. Reference 
orbit: altitude 200 km, inclination 51°. 


A considerable gain in payload mass (about 8.5 tons) is 
related to a decrease in SSTO dry mass by 18 percent. 
If a study is made of the problem with a fixed payload 
mass the effect in a decrease in SSTO dry mass will be 
more significant (according to [3], up to 22 percent). 


The results give basis for assuming that the success 
of future multiply reusable space transport systems is 
related to the use of three-component engines and other 
key technologies. This conclusion can be applied to the 
single-stage vertically launched SSTO and to aerospace 
systems of the MAKS and An-225/HOTOL-Interim 


type. 


Use of aerodynamic load-carrying properties for im- 
proving launch characteristics. The subsonic mother air- 
craft creates optimum initial conditions for launch of 
the second stage of the aerospace system as a result of 
performance of a prelaunch maneuver in the longitudi- 
nal plane. After separation the orbital aircraft performs 
flight with optimum longitudinal control, taking into ac- 
count the influence of aerodynamic lift [9]. This part 
of the trajectory and launch can be called the aerody- 
namic turn. In contrast to aerospace systems, an ordi- 
nary booster does not use aerodynamic lift and executes 
a so-called gravitational turn. Nevertheless, if the single- 
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stage SSTO carrier is a multiply reusable winged vehi- 
cle, it can execute a gentle launch after a short vertical 
segment. In this case gravitational losses are reduced 
and the consumption of fuel necessary for launch of 
a single-stage carrier into circumterrestrial orbit is re- 


Table 1 gives data on the influence of aerodynamic 
load-carrying properties on the launch characteristics 
of three variants of multiply reusable space transport 
systems (MSTS). Two possible launch methods are 
examined for a vertically launched single-stage carrier. 
The characteristics are given for a reference orbit with 
an inclination 50° and an altitude 200 km. The following 
annotations are used in Table 1: h,, V,, 4, G,, altitude, 
velocity, trajectory angle and mass of the second stage 
of the aerospace system at the point of separation or 
SSTO at the point of launch from the Earth; AV,, AV, 
, AV.., AV, losses of characteristic velocity caused 
by gravitation, aerodynamic braking, presence of an 
angle between the thrust vector and the velocity vector 
and atmospheric pressure; V,, is the ideal characteristic 
velocity; G_, is the mass of the orbital aircraft at the 
point of sustainer engine shutdown; G,, is payload m: 43 
in a Circular reference orbit. 


It must be noted that the takeoff mass of the second 
stage of the MAKS and HOTOL-Interim variants was 
assumed equal to 275 tons. The orbital aircraft rocket 
engines are fired before separation and in the prelaunch 
maneuver segment of the An-225; the second MAKS 
Stage consumes 5 tons of rocket fuel, but in the HOTOL 
variant — about 2 tons. Thus, at the separation point the 
G, values are 270 and 273 tons respectively. 
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The results presented in Table | show the good contri- 
bution of the mother aircraft to launch energetics. How- 
ever, the An-225 mother aircraft plays more than the 
role of a first stage. It is a moving launch platform 
making possible approach to the equatorial plane and 
considerable reduction of the fuel necessary for rocket- 
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the case of launch of satellites into a geostationary orbit 
by means of the transport variant MAKS-T. High op- 
erating qualities are a second fundamental advantage of 
the aerospace system. The mobility of the launch from 
the mother aircraft makes it possible to minimize the to- 
tal time expended on flight and docking with an orbital 

































































dynamic rotation of the orbital plane, for example, in object. 
Table 1 

MSTS | MAKS | HOTOL ~=—s | SsTO 

Initial conditions: 
h,> m 8,180 7,600 0 
V,, mvs 196 206 0 
O,, degrees 36 24 90 
G,, tons 270 273 $50 

Velocity losses, m/s 
AV, -700 -646 -979 -1224 
Av... -312 -271 -267 -93 
AV... -114 -87 -304 -286 
AV, -42 -45 -69 -62 
Vi» avs 8505 8419 8996 9070 
G,,, tons 37.9 38.7 68.7 67.7 
G,,, tons 8.49.5 8.5 7.$ 

Crew 2.0 unmanned 











The aerodynamic characteristics of the orbital aircraft 
exert a substantial influence on the launch trajectory and 
on the working fuel supply. The total mass put into an 
elliptical orbit is dependent on the aerodynamic quality 
[9]. As a result, the HOTOL-Interim variant has a higher 
launchable mass G__,, but the MAKS, due to the effect 
of the jettisonable fuel tank, ensures a higher maximum 
payload G,, . 


A vertically launched winged single-stage carrier has a 
possibility for improvement in launch characteristics by 
use of an aerodynamic turn. The gain in payload mass 
in this case is about 1 ton or 13 percent in the case of 
a stipulated takeoff mass of 550 tons. 


Range of attainable orbital inclinations and altitudes. 
The relative payload mass in a circumterrestrial refer- 
ence orbit is characterized by a ratio G,/G,., where 


G,,. and G,,. respectively are the payload mass and the 
total mass of the orbital aircraft in a circular reference 
orbit with an altitude h° = 200 km and an inclination i 
= 51° (the asterisk denotes parameters applying to the 
indicated reference orbit with an altitude 200 km, but 
with other altitudes the payload masses G,, are denoted 
without an asterisk). The payload mass in higher orbits 
with altitudes h > h* is dependent on the dimensionless 
parameter G/G, ., where G,, is the dry mass of the 
orbital aircraft, including the structure, fuel remnants in 
the orbital maneuvering engines and jet control system 
and crew. The values of the dimensional parameters G 
a/G,,. and G,/G,,. are 0.320 and 0.648 respectively for 
a manned orbital aircraft of the MAKS system with a 
crew of two persons, but for the unmanned HOTOL — 
0.143 and 0.825 and for the SSTO — 0.112 and 0.856. 
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Table 2 
MSTS MAKS HOTOL SSTO 
G,, kg 26,370 38,033 67,081 
Gy kg 8,440 $,440 7,540 
pe 2 10°, km! 7.975 17.846 22.786 
-v), kg/cm 6.730 9.706 17.119 














ON fd ave is =1-¢«-G,_/G,., where « = 
AV. BJ: & she characteristic velocity reserve 
expended on maneuvering in a circular reference orbit 
and on the braking impulse for descent from an altitude 
r anh Pepe m/s*; J, = 319 s is the 

of the orbital maneuvering engines; « = 
0.032 is a dimensionless parameter which is identical 
for all orbital aircraft variants. 


For the case of orbits with a great altitude h > h° it is 
possible to obtain a linear approximation for a payload 
mass G, /G,,. = | - « - GJG,,. - Ah, where Oh =h 
-h’, A= = cig , © = const. For all variants of orbital 
aircraft it can be assumed that 3 = 0.0002552 km. ,. 


The change in the relative payload mass with altitude 
is vie = | - ph, where p = AKG ,4G,,.). The 
t of change in dimensional values of payload 
mass relative to orbital altitude h is dependent on the 
parameters 3 and G,.; v = 6G,/bh = - JG... 


Table 2 gives the values of the # and wv parameters 
characterizing the influence of orbital altitude on pay- 
load mass. These parameters were obtained using pre- 
cise mathematical simulation. 

Figure 2,a shows the precise dependencies of relative 
payload mass on orbital altitude obtained as a result of 
trajectory and weight calculations (1 — MAKS, 2 — 
HOTOL, 3 — SSTO). 


Due to the jettisonable external fuel tank the MAKS 
orbital aircraft may attain circular orbital altitudes up 
to 1,500 km. The small maximum altitude in the SSTO 
variant is attributable to the greater fraction of structural 
mass in the total mass of an orbital aircraft launched 
into a low circumterrestrial orbit. Thus, the presence of 


an external fuel tank, which is jettisoned into remote 
regions of the ocean, is a good solution for a first- 
Stage aerospace system. The tank does not constitute 
any ecological danger and the range of working orbits 
is considerably expanded. 


+ a change in orbital inclination i from 0 to 100° the 

values are reduced for the variants: MAKS — from 
‘dfs 10 5.399 he. HOTOL — from 7,242 to 2,239 kg 
and SSTO — from 10,480 to 2,227 kg. 


The dependencies of the relative payload mass a 
on orbital inclination i are shown in Fig. 2,b (G,,. 

the payload mass for a stipulated inclination i, A. Ie 

G,,« is the payload mass for an equatorial orbit with the 

same altitude 200 km). 


An increase in inclination results in greater losses of 
payload mass in the SSTO variant due to a greater 
total mass of the orbital aircraft, launched into a low 
circumterrestrial orbit, and the small fraction of payload 
mass in the total launchable mass. 


The transport variant of the MAKS-T aerospace system 
ensures delivery of a payload with a mass of 18 tons into 
a low circumterrestrial orbit with an inclination 51° and 
up to 20 tons into a low equatorial orbit. In the case of 
use of a space tug with cryogenic fuel components the 
MAKS.-T will make it possible to put 5 tons of payload 
into a geostationary orbit. 

Economic efficacy is a highly important criterion in 
a comparative analysis of different concepts of multi- 
ply reusable space transport systems. Expenditures on 
designing, production, testing and certification are the 
principal cost components. 





poe tele 9, 
is 


material. Copying and dissemination 


| prohibited without permission of the copyright owners. 








FBIS-UST-95-051 
18 December 1995 














SPACE SCIENCE AND ENGINEERING $3 

ar 

q 

Sd 
if 

i r = , 

as SS 

fev (2) 


Fig. 3. 1) MAKS (Russia), 2) Delta Clipper (USA), 3) HOTOL-Interim (Great Britain), 4) Zenger 
(Germany), 5) SSTO (R) (United States), 6) NASP (United States), 7) TSTO (United States); I) decrease in 
of 


specific cost 


payload launch by a factor of 4, II) same by a factor of 2, t — number of years from onset 


of development work. 


Key: 1) million, 2) t, years 


The principal characteristics of projects for multiply 
reusable space launch systems developed in different 
countries are represented in Table 3. Figure 3 shows 
the dynamics of change in total expenditures and profit 
when constructing and using different multiply reusable 
space transport systems. 


It can be concluded on the basis of data in [10) 
and research materials from the AO NIIEPU that 
the cost of development of multiply reusable space 
transport systems exceeds the construction cost by 
a factor of approximately 50-100. It must be noted 
that at the present time it is difficult to estimate the 
cost of construction of such a launch system as the 
NASP. The approximate estimates are about 50 billion 
dollars. Accordingly, 1 kg of NASP dry mass will cost 


approximately 500,000 dollars. 


A two-stage horizontally launched system is most 
desirable in the case of additional use of an An-225 
mother aircraft for other branches of the economy. 
There are good prospects for improving the launch 
characteristics and economic indices of the aerospace 
system due to future use of a superheavy mother aircraft 


with a lift up to 400 tons. A design of such an 
aircraft, designated Hercules, has been developed. It has 
a triplane aerodynamic design in two fuselages [11]. 


A more detailed analysis of the cost must take a number 
of additional factors into account. The results make it 
possible to conclude that the MAKS pays for itself most 
rapidly. A project of the type An-225/HOTOL- Interim 
is the next best in this respect. 


The development of single- and two-stage launch sys- 
tems with promising ramjet engines with subsonic and 
supersonic combustion is a task for the more remote fu- 
ture. Different variants of multiply reusable space trans- 
port systems based solely on rocket engines are proba- 
ble candidates for introduction as new space transport 
systems with a low specific cost of payload launch into 
circumterrestrial orbit. 


The use of a three-component two-mode engine of the 
RD-701 type ensures a considerable gain in payload 
mass for different transport system variants. The most 
considerable decrease in dry mass when using the RD- 
701 is attained in the variant of a single-stage carrier. 
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ss saeco semana 
Table 3 
Project Parameters | MAKS (CIS) | MAKS-T Ap-228/ lenges NASP (USA) Artane-S =| Delta Clipper 
(CIS) HOTOL (Germany) (France) (USA) 
(CIS-Gt. 
Britain) 
Launch weight, tons 630 630 630 366 180-230 620 460 
Orbis (1,°/7H, tam $1/200 51/200 51/200 28/200 28/400 28/460 90/400 
Payload mass, tons 8.4-9.5 18 5.4 7.5 9-12 3 4.5 
Crew, persoas 2 of + 2 3 2 
unmanned 
% of wtal load 71-75 92 50 32 70-73 25 60 
quanuty 
Possitility of return + + + ’ + + 
to Earth 
Cost of launch per 1,000 1,000 2,000 4,500 ? 36,000 1,100 
payload kg, dollars 
Parallax, km 2,000 2,000 2,000 2,000 ? 0 0 
Lateral range during | 2,000 2,000 2,000 2,000 2,000 2,000 
descent, km 


























fully ensure a great weight efficiency, but the construc- 
tion of such systems is complicated by the problem 
of ensuring a guaranteed engine-free return of the first 
stage to the launch site. This variant requires additional 
research. A vertically launched single-stage carrier is 
simplest in operation. The launch characteristics of such 
a carrier can be improved by use of an aerodynamic turn. 


The MAKS project is preferable for implementation 
immediate future with minimum expenditures 
and money. The An-225 mother aircraft is a 
launch platform which will make it possible to 
equatorial plane and considerably reduce 
on a rocket-dynamic rotation of the 
example, in the case of launch of 
orbit by means of the 
moving launch from an aircraft makes 
minimize the total time necessary for 
an 
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An improvement in the launch characteristics and eco- 
nomic efficacy of the aerospace system is possible by 
future use of the superheavy Hercules mother aircraft 
with a lift up to 400 tons. 


A multiply reusable aerospace system may become a 
real step along the path of next-generation multiply 
reusable launch systems. 
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Determining Structure of Meteor Streams From 
Radar Observations 
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[Article by S. V. Kolomiyets and I. A. Milyutchenko, 
Kharkov State Technical University of Radio Engineer- 
ing and Electronics, Kharkov, Ukraine; the first para- 
graph is an abstract; UDC 523.68] 


{[FBIS Translated Text] A method for determining the s 
parameter in the distribution of meteor bodies by mass 
and stream density is examined on the basis of the 
results of radar measurements of the number of meteors 
during the period of activity of meteor streams. The 
described method is used in analyzing data from radar 
observations oj the meteor streams n-Aquarids and 
Orionids at Kharkov in 1986. A value of the parameter 
s = 1.65 was determined for the dates of the maximum 
activity of both streams. The dependencies of the density 
of the stream of meteor bodies and the parameter s on 
solar longitude and the time of day when the stream is 
observed also were obtained. The mean density of the 
stream of meteor bodies with M 2 10° g for n-Aquarids 
is 2.9 x 10" m? s'; for Orionids it is 2.5 x 10" m s'. 
In a study of the structure of meteor streams the problem 
always arises of separating a stream from the sporadic 
background. When making observations at Kharkov the 
registry is in the range of the faintest meteors 12™- 
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14" (the mean daily number of meteors is several 
tens of thousands) and the sporadic meteor background 
becomes very strong. Accordingly, in addition to a 
formal determination of the days of the stream and 
background, on the basis of data in the literature the 
background number of meteors N__, is determined by 
using data obtained on different days directly adjacent tc 
days of stream activity and on more remote days (wher 
such data are available). 


When measurements are made of the number of radio 
meteors stream activity is manifested in the fact that 
there is a change in the shape of the diurnal curve of 
the number of meteors in comparison with the results 
of observations of sporadic meteors only. The time of 
visibility of the stream radiant is computed on the basis 
of the location of the observation point and the antenna 
directional diagram. The greatest difference of the N(t) 
curve (N is the hourly number of registered meteors) 
for the dates of stream activity from the corresponding 
curves for the other cited dates should be observed 
precisely at these hours. The differences in the diurnal 
dependencies of the number of meteors for other time 
intervals characterize the range of random fluctuations 
of the sporadic background. 


The mean number of sporadic meteors during the period 
of stream activity is determined separately for each hour 
using the formula N, = N,, + Ni 2/2, where NL, 
, and Nv, , are the hourly numbers of meteors on the 
days preceding the period of stream activity and on the 
days which follow respectively. 


Specifically the dates of observations which can be 
adopted as background dates are established on the basis 
of an analysis of the diurnal curves of the number of 
meteors and the dispersion of the results obtained later. 


The number of meteors of the observed stream is 
determined by subtracting the background number N._, 
for each hour from the hourly number N,, of meteors 
registered during the time of stream activity [sb is 
for the combination stream + background]. Negative 
N,, values have no physical sense and simply indicate 
random fluctuations in the number of meteors relative 
to the mean level, adopted as the background. The 
mean daily hourly number of faint meteors (12™ - 14 
™) on days when a meteor stream occurs may almost 
not be discriminated against the sporadic background, 
in contrast to observations in the range 5™ - 7”. 


The number of meteors was measured at Kharkov at 12 
threshold levels differing from one another with respect 
to the amplitude of the registered signals by a factor of 
1.6. The measurement data at the first, most sensitive 
registry level are usually not included in the processing 
because the influence of different kinds of interference 
on the measurement results is observed at this level. 
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The results of number measurements make it possible 
to determine the parameter k, of the distribution of 


radio meteors by amplitude and the parameter s of the 
distribution of meteor bodies by mass (2). 


A theoretical formula for the distribution of the ampli- 
tudes of signals reflected from unsaturated trails was 
derived by T. Kaiser (8). It was shown that in the case 
when the reflection occurs from the point of the ion- 
ization maximum and the maximum amplitude is pro- 
portional to the linear electron concentration, the distri- 
bution of signal amplitudes conforms to the Pareto law 
with a parameter related to s by the relation s = c + 1. 
This same law also is derived in a case when the re- 
flecting points are uniformly distributed along the trail, 
but the shape of the ionization curve is not dependent 
on the mass and remains constant. 


For assessing c it is customary to convert to an integral 
distribution of amplitudes. Then the curve constructed 
on a logarithmic scale approaches a straight line with a 


negative gradient. 
k,=s-1. 


In this case the recorder is a statistical analyzer per- 
forming the quantization of amplitudes of the reflected 
signals and accumulating the numbers of events which 
exceed some threshold level. A k, assessment can be 
made by the least squares method or by a method which 
is based on the fact that with a changeover from the vari- 
able x to the new variable y = Inx the Pareto distribution 
is transformed into an exponential distribution. For such 
a distribution an assessment of the k, parameter is found 
using the maximum likelihood test: 


b= (25) 


This method is more convenient when automating mea- 
surements. When there is a detector with a logarithmic 
characteristic and a sufficiently broad dynamic range the 
problem of finding a k, assessment is essentially reduced 
in actuality to computation of the mean arithmetic value 
and is realized even without using a built-in computer 
operating in real time. For stream meteors the K, . pa- 
rameter is determined for each hour using the formula 
(1) 
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Ky we = Nak - Neca: aN 

where k, ,., is an amplitude distribution parameter 
for sporadic meteors (background); k, ,, is this same 
parameter for the totality of stream and background 
meteors. This formula is a special case of the formula 


Not = Doki: 
(a 


The k, , and k, .., values are calculated using the 
formula 


pd - T55y, 





4, © Shh, (1) 
De- (Zs) 


where y, = lg N_ is the logarithm of the integral number 
at the corresponding registry level. 


In the further processing of meteor stream measure- 
ments use is made of the dependencies s = {(k,, V), 
obtained by the statistical simulation method on a com- 
puter [3]. The s parameter for stream meteors is deter- 
mined using the dependence s = f(k,)l,, computed for 
the velocity V of meteors of this stream (Fig. 1). 


Simultaneously the s parameter and its diurnal variation 
are determined for the sporadic background and for 
the totality of the background and the meteor stream. 
Usually there is a great scatter of values of the s 
parameter for stream meteors at different hours, on 
different days and in different observation years, in 
contrast to the s(t) curve for the sporadic background. 
This is associated not so much with the structure 
of the streams as with measurement errors due to 
difficulties in discriminating the stream against a strong 
and fluctuating background. At the activity maximum of 
the streams the stream parameter s is usually less than 
for sporadic meteors. By averaging the s parameter for 
those hours when there is a maximum number of stream 
meteors the mean s value is determined for different 
observation days and a graph of the dependence of s 
on solar longitude A is determined. This dependence 
Can serve as a cross-sectional characteristic of stream 
structure. 


In order to determine the relative change in the spatial 
density of meteor bodies in the swarm cross section 
there is need only to measure the total number of 
meteors during those same hours of the day when the 
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stream is visible (or the mean hourly number). Then 
> cite af ceanemmiente an dition de give a (| Hasazo » (1) 





transverse profile of the stream under the condition that 
the instrument parameters remain stable; such profiles 
are constructed for different amplitude levels N, = f(A). 



































For the absolute density of the stream of SANEDIX (2) 
meteor bodies Q(M* with the mass M > M” we deter- T = 

mine the effecting collecting surface of the antenna S_,, 

in the plane perpendicular to the direction to the stream Baus6op amaxmsepycnbix 

radiant (in the echo plane) and the minimum mass M”° SMEANTYRRDIX yposuch (3) 
of the registered meteor bodies. S_,, is a function of time 

of day t because the horizontal coordinates of the stream | t 
radiant change continuously. The stream density in this Bus6op Noi), opearposssm 

case is determined using the formula | Nught), Nelt) (4) 
Q(M") = N(tVS, q(t). (2) 

if the number of meteors is known. S.,, and M” are 

determined when calculating the geometric and physical 
detectability factors. On the basis of the dispersion Paceer k,_(0), k,. 

of Q(M” values, computed for adjacent hours of the (5) 
very same day, it is possible to judge the measurement | 

errors and the correctness of allowance for the geometry 


of several hours of observations the spatial density 


ili . = 
factor of meteor detectability because during the course | Paceer f,_(f) 
of particles in the meteor swarm can be considered 







































































constant. 
Onpenexeune s,(f) no s(k,, V) (7) 
ae : 
| Pactet TCKYIUEX KOOPAEEAT 
PORNANTA DOTOKA OG, f, (8) 
3.0} 
23+ 
—— (9) 
at | | PacveT RSOTHOCTE BOTOKS 
wereopunm Tex c M 2 hf* (10) 
1S t 
M= 10° r, 
- lpm) -* O00) r jan 
0 0S 1.0 1.$ 
| A 
Fig. 1. Dependence of s parameter of distribution of C Konex +) (12) 
meteor bodies hy mass on parameter k_ of 
distribution of amplitudes of reflected signals. Fig. 2. Algorithm for processing results of 
Key: 1. s°! sete ~~ ~~teedioe 


The algorithm for processing the results of observations 
of the number of meteors is illustrated in Fig. 2. 
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Key: 1. Beginning; 2. Input of a priori astronomical 
data; 3. Choice of analyzed amplitude levels; 4. Choice 
of N._,(t), formation of N,,(t), N,,(t) files; 5. Compu- 
tation of k, ,(t), k, ,,(0); 6. Computation of k, ,,(0); 7. 
Determination of s,,(t) from s(k,, V); 8. Calculation of 
momentary coordinates of stream radiant a, t,; 9. Calcu- 
lation of S_{a, t,); 10. Calculation of density of stream 
of meteor with M M”; 11. reduction to mass M 
= 10° g, Q(M") — Q(M); 12. End 


The cited method was used in processing the results 
of radar observations of the number of meteors in 
1986 during the period of occurrence of the streams 
n-Aquarids and Orionids [5]. A computation method at 
two levels (second and fourth) was used for assessing 
k, , and k, _,. In this case (1) is written as 





k, =25x 
NG (IgNS - NG) — Ne (lgNy" - I8Ny") 
x —2 bn ba 2 ar 2 am, 


where the superscripts (2) and (4) denote the number 
of the level. [The subscripts are for “stream,” "stream- 
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background" and “background,” as used in the preceding 
formulas.) 


The following dates were adopted for the background: 
27 and 30 April, 14, 15 May — for »-Aquarids; 13, 
18, 25, 26 October — for Orionids. For the principal 
maxima — 6 May and 24 October respectively. The 
values of the coefficient k, .. were computed for stream 
days for each hour during which the radiant was located 
in the visibility region of the meteorological radar. 
The k, ,, value for a given day was determined by 
averaging of the weighted hourly values k, .. The s(k,, 
V) dependence with a known velocity of meteors in the 
Studied streams V = 64-66 km s"' was used to construct 
a file of values of the s parameter. Figure 3 shows 
the dependence of the s parameter on time and solar 
longitude A for meteors of the streams and sporadic 
background. 


The curves confirm the comment made above that there 
is a great scatter of s parameter values. For calculating 
the stream density of meteor bodies it is customary 
to use the minimum s value near the stream activity 
maximum. Using the graph in Fig. 3b for both streams 
we obtain s = 1.65. 


Table 1. Density of stream of meteor particles belonging to n-Aquarids 





















































Oeiruls | el al alianilsai g 
M2Me pam 12 {07% 13 | 2s |12 - |0% |i4 Jis j13 
M2M, 0.02 | 013/002 | 004/002 | - | 0012] 0.023/ aos | om 
Mame Qxiewie'| age 064 (047 [031 | 15 |072 [036 | 022 | 056 |a72 | 061 
M2M, 0.01 | 0.008 | 0.005 | 0.025) 0.012 | 0.006 | 0.006 | 0.002 | 0.012 | 0.01 
MMe op as {38 | - | [as - - }42 17 
MzM, 007 }006 | - | ama} - | - | - |eo7 |al2 
Mz + M0a/e 16 }13 [ir | 50 132 jase jos fur fis | 1.76 
M2M,| & 0.028 | 0.021 | G.018 | 0.083) 0.052 | 0.006 | 0.001 | 0.018 | 0.028 | 0.029 














is report may contain copyrighted maternal. Copying and dissemination 
ppt eter he Sai oft the ne | 














FBIS-UST-95-051 
18 December 1995 


The equatorial coordinates of the stream radiant — the 
hour angle t, and declination 6 — are computed us- 
ing well-known astronomical formulas for the radiant 
visibility region. The effective collecting surface of the 
meteorological radar S_.(t,, 4) determines the geometry 
factor for detectability of this radiant. The geometry fac- 
tor curve for the Kharkov radar complex [7], computed 
by the graph analysis method, was employed when us- 
ing formula (2) for calculating Q(M*). By knowing the 
minimum admissible registrable mass M* and the s pa- 
rameter it is possible to determine the density of the 
stream of meteor bodies with M M, (for example, M, 
= 10° g) using the formula 


The minimum linear electron density of meteor trails 
a, Was Calculated using the well-known L-C method 
cited in [4]. In this case allowance is made for the de- 
pendence of meteor altitude H on velocity. (According 
to [6], with V = 65 km s H = 97 km.) Then for the sec- 
ond amplitude level (6.4 j:V) and a transmitter power 
390 kW a,,, = 9 x 10° cm", which on the scale of radio 
meteor masses [4] corresponds to the maximum admis- 
sible registrable mass M” = 1.8 x 10° g. 


The results of calculations of the stream density of 
meteor bodies assigned to the »-Aquarids and Orionids 
for M* and M, = !0° g are given in Tables 1, 2. In 
the tables for several \ values there are no data due to 
strong interference. Curves of the dependence Q(A) for 
the studied streams are shown in Fig. 3b. 
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Table 2. Density of stream of meteor particles belonging to Orionids 
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Optimum Synthesis of Multiray Optical ferent, but then the spatial spectral components of the 
Interferometer For Making Observations Through signal, measured by pairs of subapertures with identical 
inhomogeneous Atmosphere base distances, are “separated” on the secondary (exit) 
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[Article by L. Ye. Kopilovich, Radio Physics and 
Electronics Institute, Ukrainian National Academy of 
Sciences, Kharkov; the first paragraph is an abstract; 
UDC 520.274) 


(FBIS Translated Text) Problems in synthesis of the 
entrance aperture of an optical multiray interferometer 
ensuring single and double coverage of a maximum 
range of spatial frequencies with a stipulated number 
of subapertures and a corresponding exit nonexcess 
aperture of minimum size are examined. 


The nonoptimality of a traditional optical telescope 
when making observations through an inhomogeneous 
atmosphere has motivated researchers to seek new 
image-forming methods more resistant to the destruc- 
tive influence of atmospheric phase distortions. For this 
purpose it is possible to use an optical multiray interfer- 


(entrance) aperture which is excess in the sense that the 
base distances between the subapertures are not all dif- 


aperture plane in such a way that the exit aperture is al- 
ready nonexcess with the same number of subapertures 
in it as in the entrance aperture. This makes it possi- 
ble to extend the "phase closure” method [18, 8], which 
earlier was employed when using nonexcess aperture 
masks and therefore could be suitable only for retriev- 
ing the images of bright sources, to the case of weak 
sources when the influence of noise is extremely appre- 
ciable. Different aspects of use of a multiray interfer- 
ometer are discussed in [3]. 


However, until now no study has been made of the 
problems involved in optimizing an interferometer of 
such a type. This article is devoted to this subject. 
Here there are two optimization problems. First it is 
necessary to ensure the coverage of a maximum range 
of spatial frequencies (that is, registry of all frequencies 
in this region of the u, v plane) with a stipulated number 
of subapertures in the interferometer entrance aperture; 
second, it is necessary to construct a nonexcess mask 
of minimum size with the same number of subapertures 
on the exit aperture plane. These problems are solved 
in the article by the methods developed in (2, 14] and 
(1, 13) respectively. 


Bases for rectangular grids. We will introduce the con- 
cept of a base for a rectangular grid (2, 14]. For this 
purpose on the aperture plane we will construct an in- 
tegral grid and within it we designate a rectangular grid 
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measuring (n, + 1) x (n, + 1) (with side lengths n, and 
n, and on this grid we place K point elements (subaper- 
tures). The interval between two elements (points of in- 
tersection, subapertures) is called their vector (without 
allowance for its direction). If all possible intervals be- 
tween grid points of intersection are encountered among 
these elements, then we say that these elements form a 
bese for the grid. 


The problem of optimizing the interferometer entrance 
aperture thus essentially involves construction of a base 
from a stipulated number of elements which would 
cover a grid of maximum size, that is, have a minimum 
excess (which ensures interferometer transfer of the 
maximum range of spatial frequencies). The excess of a 
base of K elements for a grid measuring (n, + 1) x (n, + 
1) is determined by the ratio of the number of intervals 
between the base elements to the number of different 
intervals present in this grid, that is, by the quantity 


a” = K(K — ')/[2(2aym, + m, + m)) © K"/(4aym) 


where X >> |, n, + n, << nn,. The following quantity 
(2, 14) will be used as a characteristic of base excess 


a=Kivayn, - @ 


and from the mathematical point of view is more 
convenient for research. 


The simplest method for obtaining the base for a 
rectangular grid involves multiplication of the two bases 
for the segments which are the sides of this grid. 
Such one-dimensional bases (sets of K, whole numbers, 
among the differences between which are encountered 
all the numbers |, 2,....n, (5]) having a minimum or 
smali excess (characterized by the quantity 


a, '‘o K,/¥m) 


found in [15] for ali K, s 11 and in [16] — for K, < 19. 
The one-dimensional bases with a large number of 
elements can be obtained in the form [5, 15] 


fet... .ky @ 


PM +d}, Jody... Kyi 
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where ‘y, is the one-dimensional base for a small segment 
(we will call it the initial base) and d, is a cyclic 
difference set (CDS) with the parameters V,, k,, A, = 
1. In a general case the CDS with the parameters V , 
k,, A, is determined in the following way [9, 10): this 
is such a set of k, whole numbers on a segment of the 
length V, that any whole number 0 < wv < V, has A, 
different representations in the form 


vad, — d,mod V,. @) 
A base of the type (2) covers a segment of the length 


Balm + I-11, (4) 


where r, = d.., - d,,,- As 7, in (2) it is most common to 
use the nonexcess base 


{0, 1, 4, 6}, n, = 6. (5) 


However, since CDS with A, = | exist only when k, = 
q + | (where q is the power of a simple number) [9], 
for constructing the base (2) with a stipulated number 
of elements as the initial bases it also is necessary to 
use other bases from [15] or [16]. 


Optimization of a base in the form (2) essentially 
involves choice of an initial base with a minimum 
excess (with a stipulated number of its elements), and 
also, as can be seen from (4), choice of the most 
compact CDS. The procedure for obtaining a CDS with 
stipulated parameters was set forth, for example, in [12); 
CDS tables with k, s 100 are given in (10). A CDS with 
a minimum r, can be determined using the procedure 
described in (2, 14] (for k, < 48 tables of such CDS 
with A = 1 were given in (6, pp 169-171]). With an 
optimum initial base and CDS, using formula (2) we 
obtain a base for a segment of maximum length (4). 


Multiplication of two bases of the type (2) gives a base 
in the form 


{(n¥s +4,0V,+¢ ‘)) (6) 
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for a rectangular grid. In (6) taking 7, = 4°,, we obtain 
a base for a square grid with the side length (4). For 
example, we obtain an optimum two-dimensional base 
of 96 elements, in (2) taking the initial base (5) and a 
CDS with the V, = 553, k, = 24, A, = 1 from 
{6}. We find from (4) that R, = 3,445 and the base (6) 
for a square grid in this case consists of 9,216 elements, 
and for it a = 2.675. 


It also was demonstrated in (2, 14) that by using the 
two-dimensional initial base +,, 6, and a two-dimensional 
difference set (TDDS) (a, b, with the parameters v,, v,, 
k, A, (being a generalization of the CDS) it is possible 
to construct a base for a grid with the side lengths 


R, = (n, + lve - 7,-1.8 = 
(m + I - ry - 1, am 


where 
Ye Gane ~ Gams ls Bang = Day (8) 


This base consists of K = kK, elements (where K, is 
the number of elements in the initial base) and has the 
form 


(Om +e Ime wr” 


Bases of the form (9) have, generally speaking, a lesser 
excess than the base (6); however, when A, = | they 
cannot be constructed for strictly square grids and it is 
necessary to select the characteristics of the initial bases 
and TTDS parameters in such a way as to compensate 
approximately for the difference between R, and R.,. 


It must be noted that a small part of the elements of the 
bases (6) and (9) goes beyond the limits of the grids 
covered by them and therefore the size of the integral 
grid on the aperture plane must be somewhat greater 


than for the grid proper. 


Generalization: multiple bases. The statement of the 
problem formulated in the preceding section can be gen- 
eralized in the following way. We will assume that it is 
necessary io ensure total coverage of a stipulated range 
of spatial frequencies when one of the elements mal- 
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functions, which is especially significant with respect to 
the reliability of a space-based interferometer. In this 
case coverage can be ensured by means of a double 
base. In general such a set of point elements that any 
interval between the grid points of intersection is en- 
countered at least A times among the intervals between 
the elements of this set can be called a A-multiple base 
for a rectangular grid. A double base is obtained when 
A = 2. 


When constructing multiple bases with a stipulated 
number of elements the optimization problem also 
arises, and specifically, coverage of a maximum range 
of spatial frequencies. For a well-known system of 
radioastronomical arrays, the Mills cross [17] (see Fig. 
1), constituting a double base, the quantity (1) is equal 
to 


a = An, + n)/Va,n; 2 4, (10) 


attaining a minimum in the case of square grids and 
increasing with an increase in the n,/n, ratio. By means 
of methods examined in the preceding section it is 
possible to obtain double bases with an excess less than 
for Mills crosses covering the corresponding regions of 
spatial frequencies. The next section is devoted to this 
problem. 


Construction of double bases. Double bases can be 
obtained both by multiplying single and double linear 
bases and on the basis of two-dimensional difference 
sets (TDDS) and two-dimensional initial bases. 
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Fig. 1. Diagram of Mills cross with central element. 
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t of the length n, = 3(K, - 


4) When K = 42 the of the bases which we found 


doubly cover a segment of the length n, = 25 and has 
the form = 
( 


{0, 2, 3, 4, 6, 12, 13, 19, 20, 24, 25, 27), 


but when K, = 13 a segment of the length n, = 30 is 
covered, but the base itself has the form 


{0, 2, 3, 4, 6, 13, 17, 20, 24, 25, 29, 30,32). (13) 
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If any of the bases (11)-(13) is multiplied by the single 
base (5), for the determined double two-dimensional 
base 


a= K,KJVim <4, 


that is, less than for the Mills cross. However, by 
selecting pairs of segments covered by these bases 
and the single bases from [15] it is possible to obtain 
double bases for nearly square grids (Table 1). The last 
column gives the gain in the number of elements (M) 
in comparison with the corresponding Mills cross. It is 
obvious that the gain will be considerably greater for 
more extended grids. 





Tabie 1. Characteristics of bases of type (6) for small grids 



































LS a (n, + 1) x (m, + 1) M 
32 3.771 13%7 4 
35 3.889 10% 10 1 
4x 3.843 13x14 2 
54 3.868 16x 14 2 
70 3.889 19x 19 2 
ld 3.920 22 x 25 2 
96 3.918 26 x 25 2 _ 





3. Two variants also are possible when constructing 
double bases on the basis of TDDS: using a double two- 
dimensional initial base and a TDDS with A, = 1 or a 
single two-dimensional initial base and a TDDS with A, 
= 2. In the first case as the initial base it is possible to 
take one of the bases found in #1, but TDDS with A, = 
1 are obtained from the CDS with A, = | with the same 
k and with V = v.v,, (v,, v,) = 1, using the formulas 


a, © d, mod v,, b, = d, mod v, 


(14) 
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By optimizing the TDDS by the method described in (2, 


SPACE SCIENCE AND ENGINEERING 65 


them, as well as the parameters of the corresponding 


14) it is possible to construct double bases for a number § TDDS and the initial bases, are given in Table 2. 


of rectangular grids. The characteristics of several of 





Table 2. Characteristics of bases of type (9) and parameters of corresponding TDDS 


























K a (R, + 1) x (R, +1) (v,, v3 &) (m, + 1) x(a, +1) M 
320 3.856 83 x 85 (13, 7;10 7213 12 
432 3.846 120 x 107 (19, 7;12) 7x16 13 
2,176 3.759 567 x 593 (93, 49;68) 7213 140 
2,252 3.795 591 x 588 (97, 39;62) 7% 16 122 











With respect to a TDDS with A, = 2, only two types of 
such sets are known: with k = 6 and v, = v, = 4 with 
this same k and v, = 8, v, = 2 [20]. These TTDS cannot 
be obtained from the CDS. A TDDS with minimum r, 
and 1, in the first case has the form: 


{O.O.. 0. .L BAB ogy 


If the initial base is taken in the form of the product 
of two bases (5), using formula (9) we obtain a double 
base for a square grid with a side length R = 25 and 
with a = 3.840, represented in Fig. 2. In the second case 
the TDDS with minimum r,, r, has the form: 


(00. 0). 63.0), 4.0.15.0.0.0.13.%. a6 


and the initial base is taken in the form of the product 
of the base (5) and the base from [15] with K, = 8: 


(0. $. i3. 19, 22, 24.32.39}, a= 24. «17) 


As a result we obtain a double base for a grid with side 
lengths R, = 50, R, = 48 and with a = 3.918. 
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Fig. 2. Double base for square grid with side length 
R = 25 (enclosed by solid line), constructed on basis 
of two-dimensional difference set. 
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Optimization of exit aperture. The nonexcess configu- 
rations (NC) of an interferometer exit aperture with a 


stipulated number of subapertures K positioned in a K 
x K grid S, can be computed by the method developed 
in [1, 13}. As a point of departure we use the set [11] 


A. @modp):. 129. 1.....9- 1. (18) 


where p is a simple number. The elements of the set 
(18) form a NC in the S, grid, the same as the elements 
of any set obtained from it with cyclic shifts along the 
grid axes or with multiplication of the coordinates of the 
elements of the set (18) by (mod p) for pairs of numbers 
reciprocally simple with p. It was demonstrated in [1, 
13} how by operations of this type the table obtained 
from (18) (if the value | is assigned to the set elements 
and a 0 is assigned to the remaining grid points of 
intersection) can be reduced to such a form that by 
successive deletion of pairs of lower rows and right- 
hand columns from it both the order of the table and 
the number of set elements in it is reduced by unity, 
that is, a NC of n elements on S, is obtained with n = 
p-1,p-2,...p-¢+ 1, where 


c= max(l,, - 4.4), (19) 


and |, is a series of i*mod p values written in an 
increasing order. In actual practice, proceeding on the 
basis of the closest simple number p > K by such a 
method it is usually possible to obtain an NC of K 
elements in the S, grid. In the example considered above 
the number of NC elements must be equal to K = 9,216 
and therefore in formula (18) we assume that p = 9,221 
(for example, see the table in [7]). In order to obtain the 
sought-for NC it is necessary to reduce the order of the 
matrix (18) by 5, which is known to be less than the 


quantity (19). 


Summary. The methods for optimizing optical interfer- 
ometer configuration developed in this article evidently 
also will find application when developing large mu:!ti- 
clement radio telescopes operating in the millimeter and 
submillimeter ranges, in the construction of which it is 
necessary to reckon with the high cost of each structural 


component. 


The author expresses appreciation to Yu. V. Ko- 
rniyenko, who drew his attention to the problem of 
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synthesis of optical interferometers, as well as making 
a number of comments on the text. 


The work was done with the financial support of the 
State Fund for Fundamental Research, Ukrainian State 
Committee for Science and Technology (Project No 2.3/ 
474). 
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Eko-Psi-95 Experiment in Mir Complex Described 


964D0110A Moscow SEGODNYA in Russian 
2 Nov 95 p9 


[Article by Viktor Gritsenko} 


{FBIS Translated Text] A few days ago the State Center 
of Medical-Biological Problems (formerly the Institute 
of Medical-Biological Problems) began the Eko-Psi-95 
experiment (ecology and psychology). Three volunteers 
will face total isolation for 90 days in an experimental 
complex on Earth which recreates the Mir manned 


orbital complex. 


The main goal of this research is to study psychological 
compatibility, mechanics of conflict situations, and 
methods of getting out of them. As the scientific director 
of the experiment, Vadim Gushin, explained, “they 
will study the dynamics of psychological comfort and 
productivity.” 
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The experiment was prompted by the fact that while 
working in space, despite the existing special psycho- 
logical selection process, the crew cannot always man- 
age to avoid problems. In 1976 the flight of the mili- 
tary orbital station Almaz-3 (officially called Salyut-5) 
was concluded prematurely because of the psychologi- 
cal incompatibility of military cosmonauts Boris Voly- 
nov and Vitaliy Zholobov. Or a very recent example. 
Viadimir Dezhurov and Gennadiy Strekalov were part 
of this fall's eighteenth main expedition with American 
Norman Tagard on the Mir orbital compl: .. The cos- 
monauts were fined 15 percent of the to' sum of the 
contract for “violation of seniority in the _rew." 


In order to avoid such complications, the State Center of 
Medical-Biological Problems decided to "repeat" them 
on Earth. Moreover, the crew's problems will be created 
intentionally. On the twenty-second day of the “flight’ 
there are plans to “dock” with a ship whose crew will 
“visit.” However, at this time the automatic controls will 
"fail" and the "docking" will have to be done manually. 
The cosmonauts will spend 60 hours doing this, during 
which they will not have the opportunity to even 
sleep. When the “docking” finally occurs, three more 
researchers will arrive and live an entire week there. The 
specialists promise to provide the “cosmonauts” with 
many other problems. 


The State Center of Medical-Biological Problems has 
its own team of seven cosmonaut-researchers, of which 
only Valeriy Polyakov has been fortunate enough to fly 
in space, twice. Thus, the center is especially famous 
because those who have not flown in space participate 
in the ground-based experiments as a consolation here. 
In September 1994-January 1995 the HUBES experi- 
ment, lasting 135 days, was carried out under contract 
with the ESA. Three representatives of the center partic- 
ipated: cosmonaut-researchers Vladimir Karashtin and 
candidate of medical sciences Vasiliy Lukyanyuk, as 
well as candidate of medical sciences Igor Nichiporuk. 


The experienced 48-year-old Colonel Aleksandr 
Andryushkov is participating in the current experiment. 
He is a correspondent of KRASNAYA ZVEZDA who 
underwent training at the Mission Control Center in 
1990-1992 in preparation for space flight as part of the 
team from the Journalists Union of the USSR. How- 
ever, he was not the first Soviet cosmonaut-journalist 
(he was preceded by a Japanese journalist). Also 
participating in the experiment are two students of the 
Moscow Aviation Institute, Aleksandr Ivyanskiy and 
Yaroslav Balakhontsev. At the end of the experiment 
each of them will receive $1,500. 
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Mir Leak Described 


964D0110B Moscow SEGODNYA in Russian 
2 Nov 95 pl 


[Article by Mikhail Chernyshov] 


(FBIS Translated Text) Tuesday evening Mir teleme- 
try recorded a drop in temperature and a small leak 
of heat transfer fluid in the station’s heating system. It 
seems the incident caused a mild panic among those in 
charge of space matters. The next day a hastily created 
commission consisting of representatives of Energiya, 
the Mission Control Center, and Military Space Forces 
were at work. The situation was complicated by the fact 
that the pipelines for the station's "central heating” are 
behind panels and instrument shields. The specific loca- 
tion of the heat transfer fluid leak was not determined 
over the course of the day. The specialists stated that 
nothing strange had occurred; only a small leak in the 
heat regulation system. The vapors of the leaking heat 
transfer fluid, like automobile antifreeze, poses no dan- 
ger to the crew's respiration. All remember, however, 
the "ice age" which occurred on Salyut-7, when a sim- 
ilar “small” flaw turned the internal sections of the or- 
bital complex into an ice cave. In addition, the batteries 
discharged and the lights went out. Only desperate mea- 
sures taken by cosmonauts Vladimir Dzhanibekov and 
Viktor Savinyy, resulted in a “thawing” of the station, 
returning it to life. 


The current situation, according to the assurances of the 
Mission Control Center, is much simpler. Mir is not 
threatened with the onset of an ice age. "The heat trans- 
fer fluid leak is small, and deviations of temperature 
from the norm are insignificant.” Nonetheless, all cur- 
rent and proposed operations to repair the station have 
for some reason been kept secret and withheld from the 
press. The attempt to hide something that will come out 
anyway is completely in the spirit of the traditions of do- 
mestic cosmonautics. However, there is a small compli- 
cation: in addition to our countrymen, Yuriy Gidzenko 
and Sergey Avdeyev, who could be persuaded to keep 
quiet, ESA astronaut Thomas Righter is on board Mir. A 
Western researcher, by terms of his contract, is obliged 
to promptly report to his home agency all failures which 
occur on the station. There has been no talk of reexami- 
nation of the program of operations or emergency evac- 
uation of the crew from on board Mir yet. The specialists 
promise to finally put the situation in order today. 


Russian Space Industry Described 
964D0122A Moscow KRASNAYA ZVEZDA 
in Russian 11 Nov 95 p 5 


(Interview with Sergey Zhiltsov, department head at the 
State Space Science and Production Center imeni M.V. 
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Khrunichev, by KRASNAYA ZVEZDA correspondent 
Mikhail Rebrov; place and date not given: “Iridium- 
Eurasia; New Company Promises Russian New Quality 
of Space Communication"; first two paragraphs are 
boldface KRASNAYA ZVEZDA introduction] 


(FBIS Translated Text) Today it is difficult to imag- 
ine our business and everyday lives without radio re- 
ceivers, televisions, telephones, and computers. Such 
technology feeds us information, connects us with the 
outside world, and makes leisure time interesting. 


Even back in the dawn of their development, hu- 
mankind had ways to send “meaningful signals" 
across great distances by using the beating of drums, 
torches, campfires, and smoke. Space has brought a 
new level of quality to the task. As they say, it seems 
like only “ yesterday" we dreamed of a time when we 
would be able to call a person anywhere on earth by 
simply dialing a number. We supposed that commu- 
nication with that person would automatically be es- 
tablished regardiess of whether he was in the middle 
of the ocean, in the center of a big city, or crossing 
the desert. It was from such a historical digression 
that I began my conversation with Sergey Zhiltsov, 
department head at the State Space Scientific and 
Production Center imeni M.V. Khrunichev—an en- 
terprise that is among other things working on de- 
signing space communications. 


[Rebrov] Sergey Aleksandrovich, in our stormy era a 
man deprived of communication is deprived of prac- 
tically everything. Information is a resource. An inex- 
haustible and important resource. Not that long ago or- 
dinary relay satellites were fantastic, whereas today we 
talk about completely new concepts on a global scale. 


[Zhiltsov] I would put it more precisely. We are not 
simply talking about them but are actually developing 
them. You are correct: Each day, common-use commu- 
nication occupies a bigger place in society. Today vir- 
tually all people have the opportunity to use telephone 
communication. But life is being dictated by new con- 
ditions. And supporting business activity now requires 
more than constantly being on the telephone. Today 
the telephone services market offers a wide selection 
of equipment needed by modern businessmen. This in- 
cludes data transmission, paging, facsimile communica- 
tion, and confidential communication. A significant shift 
in clients’ preferences in favor of mobile communica- 
tions, namely, a preference for pocket devices, has been 
observed. 


[Rebrov}] As far as I know, progress in the field of 
telecommunications and communication, as well as 
developments in rocket and space technology, enabled 
the U.S. firm Motorola, Inc., to begin developing the 
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IRIDIUM satellite system of mobile communications 
back at the end of the eighties. The system, which 
is based on the use of low-orbiting spacecraft, would 
enable subscribers to use an ordinary portable telephone 
to establish communications regardless of their location. 


[Zhiltsov] Such a project is in fact being implemented 
at the present time. It is intended first and foremost 
to provide communication services to governmental 
and business organizations and to services requiring 
inumediate communication in emergencies. I would add 
that this is extremely important to populations living in 
hard-to-reach regions and in territories without ground 
systems of collective communication. I am not speaking 
about those members of population in middle and high 
income brackets, representatives of business circles, 
travelers, and so forth. 


[Rebrov} Let's clarify something. Is this system an 
international system? 


[Zhiltsov] Yes. In 1993, an international consortium 
called Iridium, Inc., was formed to implement the 
IRIDIUM project. The work is continuing. After the 
system has been put into operation, Iridium, Inc., 
will own its space segment. The consortium includes 
leading firms in the field of telecommunications, as well 
as developers and manufacturers of communications 
equipment and rocket and space technology in the 
United States, Canada, Russia, Italy, Japan, China, and 
other countries. 


Extra capital and “brains” from such well-known firms 
as Lockheed Corporation, Siemens, Telespacio, Mo- 
torola, McDonnell Douglas, Stet, and Vebak [translit- 
eration] have been brought in to work on the project. 
Russia is represented on the list by the leading rocket 
and space company State Space Scientific and Produc- 
tion Center imeni M.V. Khrunichev (GKNTPTs]. 


{Rebrov] What is the relative extent of your involve- 
ment? In other words, what do you have to invest and 
what will you receive? 


{Zhiltsov} As an investor, the space center must, first 
and foremost, create a financial base of operations in 
Russia, purchase the required equipment and construct 
ground interface stations, and set up a network of 
system service providers. In addition, in January 1993, a 
contract was signed with the U.S. company Motorola for 
three launches of the rocket Proton to put 21 IRIDIUM 
' communication satellites into orbit (seven satellites each 
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time). The first launch is slated for the end of 1996. Our 
center is developing and manufacturing carrier rocket 
containers to hold seven satellites, as well as a system 
to separate them from the Proton’s final stage and put 
them into orbit. This is a commercial project. 

[Rebrov] What guarantees are there that it will be 
implemented? 


[Zhiltsov) A number of decisions of the State Com- 
mission on Electronic Communication and the "Memo- 
randum Concerning Mutual Understanding and the Ba- 
sic Principles of Cooperation in Implementation of the 
IRIDIUM System on Russian Federation Territory” be- 
tween the Russian Federation Ministry of Communica- 
tions and the Space Center imeni M.V. Khrunichev have 
permitted successful implementation of the project on 
Russian Federation Federation territory. 


[Rebrov] Is everything going “smoothly” as the space 
specialists say? 


[Zhiltsov) In general, yes. The amount of work sched- 
uled, short time frame, and plans to begin regular op- 
eration of the system in 1998 have all necessitated the 
creation of the new company Iridium-Eurasia, which 
should solve problems regarding the IRIDIUM system 
in Russian territory and provide the communication sys- 
tem’s services within the territories of a number of gov- 
ernments as well as in those served by the base interface 
Station that is being created. And that company has al- 
ready been created. In the future, the State Center imeni 
Khrunichev will grant the company the right to operate 
the interface station and provide the communications 
services of the IRIDIUM system, as well as equip the 
interface station. I don’t think that we will have a prob- 
lem. 


[Rebrov] Sergey Aleksandrovich, tell us a little about 
the partners. 


[Zhiltsov] The GKNTPTs’ partner is the closed joint- 
stock company Zond syvaz, which is the operator of a 
digital satellite network linking a number of Russia's 
cities. 

(Rebrov] And what's next? 


[Zhiltsov}] There are plans to create a second system 
interface station, preferably in the Asian part of Russia. 
This will make it possible to increase the number of 
subscribers served, reduce rates, and increase the quality 
of services provided. The preferred use of the system 
is to provide electronic communications services in 
rural, sparsely populated, and difficult-to-reach regions, 
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including regions of the Far North, Siberia, Far East, 
the Chukotski Peninsula, the Kamchatka Peninsula, and 
Sakhalin, most of whose population centers are more 
than 100 km from telephone stations. 


The IRIDIUM system must become a link in the world’s 
infrastructure of mobile communications networks, in- 
cations systems. We have had much of this prepared for 
many years now—both the satellite and the ground seg- 
ment. We are using existing communications lines and 
do not need to create anything new. And we are already 
inviting all those operating telephone communications 
networks to join us. 

At the beginning of our conversation, I called your 
attention to the fact that we have already advanced 
beyond the “conversation” stage and have gotten down 
to business. The Ministry of Communications, as the 
country’s leading expert, is supporting our project. It 
therefore appears highly certain that Russia will have 
modern mobile communications, which will in turn open 
fundamentally new opportunities for the system's users. 
[Rebrov] And now one final question, Sergey Aleksan- 
drovich: How soon will the investments in the system 
be recouped? 


(Zhiltsov} Ill tell you. The system's cost recovery 
period is estimated at a year and a half. After that, there 
will be a profit. 


FBIS-UST-95-051 
18 December 1995 


Russian Engines for American Rockets 


964D0122B Moscow KOMMERSANT DAILY 
in Russian 4 Oct 95 No 183, p 9 


(Unattributed article: “Ministry of Defense Does Not 
Oppose Equipping American Rockets With Russian 
Engines”) 

(FBIS Translated Text] Yevgeniy Gritsenko, general di- 
rector of the Samara NTK Dvigateli NK [NK Engines], 
denied the numerous media reports in late September 
stating that the Russian Federation Ministry of Defense 
is blocking a joint Russian-American program to use 
NK-33 rocket engines. In his words, it was precisely 
in September that the Ministry of Defense approved the 
collaboration of the Samara NTK Dvigateli NK with the 
American company Aerojet. In July, two NK-33 engines 
were sent to the United States, where final preparations 
are being made for firing tests on Aerojet engine test 
beds in Sacramento (California). The NK-33, which was 
developed more than 25 years ago within the framework 
of the program to achieve manned flights to the Moon, 
is slated for use on new-generation American carrier 
rockets and to modernize existing rockets. If the Samara 
NTK Dvigateli NK program is implemented, significant 
amounts of capital could be obtained from the sale of 
NK-33 engines that have already been manufactured and 
from license fees. That capital could then be invested 
in the development and production of promising aircraft 
engines and commercial gas turbine units. 


> 
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Degradation of the Interphase Si-SiO, Boundary 
Under Field and Radiation Effects 

964D0047A St. Petersburg PISMA V ZHURNAL 
TEKHNICHESKOY FIZIK] in Russian 26 May 95 
Vol 21 No 10, (manuscript received 7 Nov 94; after 
revision 16 Feb 95) pp 1-5 


[Article by I. V. Klimov, Yu. M. Listopadov, A. I. 
Nazarov] 


(FBIS Abstract) It is important to understand the reasons 
for the degradation of metal-dielectric-semiconductor 
transistors due to “hot carriers” (i.c., X-rays, avalanche 
hole and tunnel electron injection) to predict perfor- 
mance, optimize operating modes and improve pro- 
duction methods. This paper studies the structure of 

ine silicon-silicon dioxide-silicon with three 
types of oxides (vapor, pyrogenic, and dry) grown on 
an n-type silicon substrate. The radiation and avalanche- 
injected charge is formed in the same way in all oxide 
types. Study of the electric field strength showed that 
as the avalanche duration increases, the captured charge 
decreases, either due to leakage of the positive bulk 
charge from the trap into the electric field, or due to 
additional electron injection leading to recombination 
of electrons with the captured positive charge. A fig- 
ure shows the dependence of surface state density on 
surface potential for various effects. Tunnel electron in- 
jection and avalanche hole injection with an electron 
component, as well as emission at a positive voltage at 
the gate causes a significant increase in surface states at 
the silicon-silicon dioxide interface. The greatest den- 
sity is observed in samples with dry oxide. A study of 
the degradation of field metal-dielectric-semiconductor 
transistors after homogeneous avalanche injection of 
electrons showed an increase in the fluctuations of sur- 
face potential during injection. After a critical charge 
(10 pulses) the potential relief equalized. The most sta- 
ble were transistors with a subgate dielectric manufac- 
tured using moist oxide. There is a clear relation be- 
tween the density of surface states and the method of 
manufacturing the oxide layer when there is an electric 
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field in the dielectric and charges of both signs in the in- 
terphase carrier boundary. The effects of “hot carriers” 
are determined so thal one can predict the degradation 
changes in metal oxide semiconductor structures. Fig- 
ures 2; references 5: 2 Russian, 3 Western. 


Growth of 3D Polytype 4H and 6H Silicon Carbide 


Crystals 

964D0047B St. Petersburg PISMA V ZHURNAL 
TEKHNICHESKOY FIZIK] in Russian 26 May 95 
Vol 21 No 10, (manuscript received 20 Dec 94) pp 
20-24 


[Article by A. Yu. Maksimov, A. A. Maltsev, N. 
K. Yushin, I. S. Barash; A. F. loffe Physicotechnical 
Institute, St. Petersburg] 


(FBIS Abstract] This paper describes an experiment to 
try to grow bulk silicon carbide crystals of polytypes 
6H and 4H, as well as an attempt to obtain substrates 
(monocrystalline silicon carbide substrates) with a di- 
ameter of greater than an inch. The growth rate of poly- 
type 6H monocrystals depends on temperature and ar- 
gon pressure: too fast a growth rate and quality suffers; 
too slow and insufficient crystal height is obtained dur- 
ing the experiment. The optimal growth rate was 0.8 
mm/hr at 2230°C and an argon pressure of 7 millibars. 
Boules of the 6H polytype 25 mm in diameter were ob- 
tained. The 4H polytype crystals were grown by placing 
a thick epitaxial layer of polytype 4H on a 6H sub- 
strate with an excess of carbon in the atmosphere and 
with the addition of a group IV element (Ge, Se, Pb). 
Later the substrate was removed. X-ray analysis demon- 
strated a sufficiently good crystal quality. The boules 
were used to manufacture substrates 25 mm in diameter 
(or a square shape 30x30 mm) 0.3-0.5 mm thick. The 
monocrystal plates of polytypes 4H and 6H had a poly- 
type homogeneity of better than 80 percent, a dopant 
concentration between 10’ and 10"* cm’, a pore con- 
centration of 10? cm* and a dislocation density of 10° 
cm~. Figures 3; references 5: 1 Russian, 4 Western. 
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Diffusion Waveguides in Polycarbonate 
964D0047C St. Petersburg PISMA V ZHURNAL 
TEKHNICHESKOY FIZIK] in Russian 26 May 95 
Vol 21 No 10, (manuscript received 14 Feb 95) pp 
25-28 


{Article by A. V. Tomov, Institute of Applied Optics, 
Academy of Sciences of Belarus, Mogilev) 


sion time is from 0.5 to 16 minutes. The waveguides 
have a smooth, virtually linear function of the index of 
refraction from the surface to the substrate. There is a 
drastic change at the waveguide-substrate interface. In 
the best samples optical losses did not exceed 2.8 dB/ 
cm (A = 0.633 4m). A mathematical model is proposed 
which describes the index of refraction profile, which is 
of particular importance in the study 0 nonlinear phe- 
nomena in these waveguides. Figures 2; references 9: 4 
Russian, 5 Western. 


Effective Heating of Electrons and lons by Low 
Hybrid Waves in an FT-2 Tokamak 


964D0047D St. Petersburg PISMA V ZHURNAL 
TEKHNICHESKOY FIZIKI in Russian 26 May 95 


Vol 21 No 10, (manuscript received 16 Feb 95) pp 34-39 


[Article by V. N. Budnikov, V. V. Dyachenko, L. A. 
Yesipov, Ye. R. Its, M. A. Irzak, S. I. Lashkul, K. 
A. Podushnikova, A. Yu. Stepanov; A. F. loffe Physi- 
cotechnical Institute, Russian Academy of Sciences, St. 
Petersburg) 


(FBIS Abstract} In low hybrid heating in an FT-2 
tokamak wave absorption and ion heating is associated 
with the development of decay parametric instability of 
heating waves in the central plasma region. Generation 
of fast ions begins at the center at a minimum energy 
of eight times the ion temperature. As plasma density 
increases, the area of decay shifts to the periphery 
and local density increases slightly. The energy of the 
ions gradually increases due to the development of the 
longwave spectrum of the parametrically excited waves. 
Most ions are heated in this transition region. This is 
a detailed study of the heating of ions and electrons 
in the transition region of plasma densities. Installation 
and discharge parameters are given. The effect of low 
hybrid waves on plasma parameters and the dependence 
on density are shown in a figure. The highest efficiency 
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of heating is reached in a very narrow range of 
densities. Generation of fast ions begins at a density only 
somewhat lower than the low hybrid density. The central 
plasma region satisfies the conditions of low hybrid 
resonance and the heating can formally be considered 
classical. However, some facts point to a parametric 
mechanism, for example more slower ions than faster 
ions appear at lower densities, and the fact that the 
beginning of generation of fast ions occurs at the same 
time as a mechanism which leads to effective heating 
of electrons. It is stated that generation of fast ions and 
heating of electrons are associated with the development 
of parametric instability of the pump wave. A change 
in discharge parameters when a high-frequency pulse 
is used may Cause a transition from a parametric to a 
Classical mechanism of interaction. Figures 2; references 
5 (Western). 


Possibility of Establishing the Neutron Generation 
Mechanism in the Mechanical Effect on Deuterated 
Matter 


964D0047E St. Petersburg PISMA V ZHURNAL 
TEKHNICHESKOY FIZIKI in Russian 26 May 95 
Vol 21 No 10, (manuscript received 15 Feb 95) pp 48-52 


[Article by Ye. G. Fateyev, Institute of Applied Me- 
chanics, Urals Division, Russian Academy of Sciences, 
Izhevsk] 


{[FBIS Abstract] In 1986 it was found that neutrons are 
emitted when there is a powerful mechanical effect on 
deuterated substances, but the neutron emission mech- 
anism has remained a mystery. This paper discusses 
the possibility of establishing the mechanism of nev- 
tron emission which arises at the time of a rheologi- 
cal explosion or a rheological explosion-like process in 
impact loading and vibration disintegration of deuter- 
ated materials. It is more likely that the cause may be 
found using this type of process because a stage-by- 
Stage analysis can be done. A classical mechanism of 
neutron generation is proposed which is associated with 
acceleration processes in the plasma discharges between 
the edges of microfractures. However, this type of dis- 
charge is only sporadic and one cannot unambiguously 
identify them with neutron emission bursts. Mechani- 
cal effects and low-frequency electric fields are used to 
determine the origin of neutron emissions. This paper 
discusses the method, but does not present results of 
the application of the model. Figure |; references 18: 
15 Russian, 3 Western. 
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High-Temperature Compound Electric Conductivity of Heated Superdense Plasma 
YBa,Cu,Se, with T, = 371 K 964D0047G St. Petersburg PISMA V ZHURNAL 


964D0047F St. Petersburg PISMA V ZHURNAL 
TEKHNICHESKOY FIZIK] in Russian 26 May 95 
Vol 21 No 10, (manuscript received 7 Feb 95) pp 67-71 


{Article by V. D.Shabetnik, S. Yu. Butuzov, V. I. Plak- 
siy, Central Scientific Research Experimental Institute 
of the Ministry of Defense of the Russian Federation, 
Mytishchi)} 

{[FBIS Abstract} A compound with the formula 
YBa,Cu,Se,, (where x is disruption of the stochiome- 
try of Se, and varies from 0 to 0.5) was created to try to 
obtain a substance with a high temperature of transition 
to the superconducting state. X-ray analysis was used 
to study the compound. Specific features of its structure 
are discussed. The device used to prove the presence 
of superconducting properties in the compound is 
described. Superconductivity is signaled by a change in 
the sign of magnetic susceptibility and a drastic change 
in the resistance of the material. A figure shows the 
magnetic susceptibility and resistance for YBa,Cu,Se, 
and this is contrasted with the same information for 
YB,Cu,O,, a known superconductor. It is found that 
the transition temperature for the former is 371 K (93 
K for the latter). The higher temperature is due to the 
presence of a selenide ion Se, in the crystalline lattice. 
The volume of material in the superconducting phase 
is 12 percent. Figures 4; references 2 (Russian). 


TEKHNICHESKOY FIZIKI in Russian 26 May 95 
Vol 21 No 10, (manuscript received 18 Jan 95) pp 89-92 


[Article by V. V. Pogosov, I. T. Yakubov, Institute 
of High Temperatures, Russian Academy of Sciences, 
Moscow] 


([FBIS Abstract] A superdense heated plasma arises 
when large energy fluxes affect matter, for example, 
irradiation of metal surfaces by strong subpicosecond 
laser pulses. The metal plasma does not expand, and 
retains its initial density. It cannot be described using 
common approximations. The plasma consists of ions 
of charge Ze and electrons with the parameters of ion- 
ion, electron-ion, and electron-electron interactions. The 
ion subsystem is described by a cellular model, each 
ion in a cell screened by Z electrons. The distribution 
of free electrons, which move unhindered from cell to 
cell, is found to be nearly uniform. The cell potential 
is defined and the effect of scattering on the ion 
core is considered. The time of flight through the 
cell is calculated classically. An equation is derived 
for the coefficient of electric conductivity. Comparison 
with experimental results shows that the cellular model 
developed here provides a good description, but can be 
refined. Figures 2; table 1; references 8: 2 Russian, 6 
Western. 
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Estimation of the Thermodynamic Nonideality of 
the Detonation of an Emulsion Explosive 
957A1198A Moscow KHIMICHESKAYA FIZIKA 

in Russian Nov 94 Vol 13 No 11, pp 79-87 


[Article by V.V. Odintsov, V.I. Pepekin, and B.N. Ku- 
tuzov, Chemical Physics Institute imeni N.N. Semenov, 
Russian Academy of Sciences, Moscow; first paragraph 
is KHIMICHESKAYA FIZIKA abstract; manuscript re- 
ceived $ Jul 94; UDC 534.222.2+536.7] 


(FBIS Translated Text) An attempt has been made to es- 
timate the thermodynamic component of the nonideality 
of detonation in an ammonium nitrate-based mixture of 
emulsion explosive with glass microspheres as compared 
with the hydrodynamic component that is manifested in 
the effect of charge diameter on detonation rate. To ob- 
tain emulsion explosive and glass microsphere mixture 
detonation rates correlating with experimental data, the 
assumption was made that the reaction of products at the 
Chapman-Jouget point is incomplete. In other words, it 
was assumed that the products consist of three chemical 
subsystems and that each has its own chemical equi- 
librium established within it. The assumption of incom- 
plete reaction of the emulsion explosive and glass mi- 
crosphere mixture detonation products at the Chapman- 
Jouget point because of one of two factors (either the 
diffusion of the decomposition products of the emulsion 
matrix’s components is insufficient or else the rate of the 
chemical reaction is insufficient on account of low tem- 
peratures) results in calculated and experimental rela- 
tionships between the detonation rate and initial density 
that were consistent with one another. The order of the 
degree of thermodynamic nonideality numerically char- 
acterizing the degree of reaction of the products may be 
approximated by a linear dependence of the degree of 
reaction on the initial density with a slope coefficient 
on the order of 0.1-0.15 percent x (kg/m'y'. Accord- 
ing to the results of calculations performed by using a 
BKW-RDX model [BKW being an acronym for Becker- 
Kistyakovskiy-Wilson equations of state], the diffusion 
factor should be the main reason for the increase in non- 
ideality (decrease in the degree of reaction) as the con- 
centration of hollow microspheres in the mixture with the 
emulsion matrix is increased. When a BKW-RR model is 
used, on the other hand, the effect of chemical kinetics 
canno! not be ignored. 


Emulsion explosives have come to be used extensively 
in the extraction industry thanks to their rather high 
detonation parameters, water resistance, and explosion- 
proof nature. The mixture’s sensitivity is increased by 
adding hollow glass microspheres to the system, which 
reduces the density of the emulsion explosive and glass 
microsphere mixture. 
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Studying the process has made it possible to establish 
the characteristic features of an emulsion explosive’s 
detonation as compared with other classes of explosives. 
It has, for example, been noted' that the density depen- 
dence of detonation rate indicates a behavior that is sig- 
nificantly different from that of powerful monomolec- 
ular explosives. The critical diameter of detonation in- 
creases exponentially as the emulsion explosive’s den- 
sity increases. It has also been noted? that the detonation 
rate of an emulsion explosive peaks at densities below 
its peak density. The deviation of experimental detona- 
tion rate values extrapolated for an infinite charge di- 
ameter from the values predicted on the basis of calcu- 
lations performed by using the Becker-Kistyakovskiy- 
Wilson [BKW] program or the Kemlet [transliteration] 
method will become greater as the initial density of the 
emulsion explosive and glass microsphere mixture is de- 
creased. 


Meyder?® [transliteration] noted similar detonation rate 
behavior for mixtures of solid ammonium nitrate with 
combustible components. The decomposition and reac- 
tion of the ammonium nitrate, part of which behaves 
as an inert component, has been hypothesized to be the 
cause. The experimental values of the detonation rate in 
Amatex-40 and Amatex-20 were successfully obtained 
in calculations performed under the assumption that 50 
percent of the ammonium nitrate in the mixture reacts 
while 50 percent remains inert. The experimentally ob- 
tained values of the detonation rate in amatol are repro- 
duced in calculations if 19 percent of the ammonium 
nitrate reacts while the rest remains inert. 


In this study, an attempt has been made to estimate the 
thermodynamic component of the nonideality of deto- 
nation in an emulsion explosive and glass microsphere 
mixture as opposed to the hydrodynamic component that 
is manifested in the effect that charge diameter has on 
detonation rate. 


The study is based on the results of a published 
experimental study' of the detonation of a emulsion 
explosive and glass microsphere mixture in which 
the starting mixture was obtained by adding unsorted 
hollow glass microspheres with a true density of 150 
kg/m’ and average diameter of 70 jm to an emulsion 
with the following composition (percent by weight): 
ammonium nitrate, 77; water, 16; oil, 6; and emulsifier, 
1. The density of the matrix was specified as 1,353 kg/ 
m’, and the particle size of the emulsion ranged from | 
to 5 ym. 


Table 1 presents the composition of the emulsion ex- 
plosive and glass microsphere mixture and parameters 
of its components that were assumed in the thermody- 
namic calculations. 
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\ cmpenent Caemical Formaia Percent by Weight Enthalpy, kJ/cg 
Amsmoasum astratc H,N,O, 77 -4,563 
Waser H,0 16 -15,860 
Saturated hydrocarbon (liquid) (CH,), 6 -2,088 
Oleic acid C,,H,,0, 1 -2,841 
Giaas SiO, above 100 percent -15,160 














The calculations were performed by using Becker- 
Kistyakovskiy-Wilson equations of state with different 
sets of coefficients and the respective hypotheses re- 
garding condensed carbon: 


— BKW with an enthalpy of graphite formation of 15 
kcal/mol’; 

— BKW-RDX? with the possibility of the formation of 
disperse diamond and graphite; 


ee @<,- = 
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Figure 1. Dependence of the detonation rate of an 
mixture 
lines: 1, BKW;; 2, 








— BKW-TNT? with the possibility of the formation of 
disperse diamond and graphite’, 


— BKW-RDX” with the possibility of the formation of 
disperse diamond and graphite’, 


— BKW-R with an enthalpy of graphite formation of 
12 kcal/mol; 


— BKW-RR with the possibility of the formation of 
disperse diamond and graphite.’ 


In all of the models, the graphite was assumed to be 
compressible and its behavior was described in terms of 
a Kovena (transliteration) equation of state, whereas the 
glass and diamond were assumed to be incompressible. 


Table 2 and Figure | present the results of equilibrium 
calculations of the detonation of an emulsion explosive 
and glass microsphere mixture within the framework of 
the aforesaid models. The points in Figure | represent 
experimentally obtained/ detonation rate data that have 
been extrapolated for an infinitely large charge diame- 
ter, which will for brevity’s sake henceforth be termed 
the ideal detonation rate. The experimental value in Ta- 
ble 2, on the other hand, was obtained by a least squares 
approximation of the values of the ideal detonation rate’ : 


D, = 5.383p, - 152, 


where the detonation rate is expressed in m/s and the 
density is expressed in kg/m’. 
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Table 2. Detonation Parameters of the Emulsion Matrix With o Density of 1,353 kg/m’ 


























Medd ee — ee Pressure, GPa Temperature, K Density, kg/m’ Speed of Sound, m/ 
s 

BKW 7,584 18.41 974 1,772 $,791 

BKW-RDX 7,134 16.04 976 1,764 $,473 

BKW-TNT 6,547 14.53 1,254 1,805 4,907 

BKW-RDX" 6,992 17.73 1,126 1,773 $,335 

BKW-R 7,078 15.75 1,719 1,763 5,434 

BKW-RR 6,613 13.95 1,781 1,770 5.054 

Experiment! 7,131 




















A comparison of the lines in Figure 1 reveals that the 
density dependence of the detonation rate of an emul- 
sion explosive and glass microsphere mixture differs for 
the different models. The slopes of the dependence are 
approximately identical for all of the models with the 
exception of the BKW-RR model, in which case the 
slope is noticeably smaller. The remaining models may 
be divided into three groups based on absolute value of 
the detonation rate: the BKW model yields the greatest 
value; the BKW-RDX, BKW-RDX", and BKW-R mod- 
els yield mid-range values; and the BKW-TNT model 
yields the lowest value. The detonation rate values cal- 
culated on the basis of the BKW-RR model go from 
values characteristic for the middle group given low 
densities to the value calculated in accordance with the 
BKW-TNT model given near-peak densities. The slope 
of all of the calculated dependences is less than that 
of the experimental dependence in the case of the ideal 
detonation rate. 


To obtain calculated rates of detonation of an emulsion 
explosive and glass microsphere mixture that correlate 
with data obtained in experiments, we assumed that the 
reaction of products at the Chapman-Jouget point was 
incomplete, which is just what Meyder assumed for 
mixtures of solid ammonium nitrate. Unlike Meyder, 
who considered part of the ammonium nitrate to be 
inert, we hypothesized that the products consist of 
three chemical subsystems and that cach has its own 
The first subsystem consists of products of the chemical 
decumposition of the aqueous solution of ammonium ni- 
trate. The second subsystem consists of products of the 
Chemical decomposition of the mixture of hydrocarbons. 
The third subsystem consists of products of the chemical 
interaction of all of the emulsion matrix's components. 
The ratio of components constituting each of the subsys- 


tems corresponds to the ratio of the same components in 
the starting emulsion matrix. In all remaining respects, 
the detonation products are considered to be a unified 
system, i.c., the state variables are one and the same for 
the entire system. 


The said assumptions were all made on the basis of the 
diffusion theory of incomplete mixing of the products 
of the decomposition of the emulsion’s components and 
the incompleteness of the chemical reaction, which is in 
turn the result of the products’ rather low temperatures. 


,;*, FF fF FF F * FE 
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Figure 2. Calculated dependence of the detonation 
rate of an emulsion explosive and giass microsphere 
mixture on the degree of reaction of the products 
given an initial density of 1,000 kg/m’. (Y-axis 
represents detonation rate in m/s.) 


It is proposed that the weight percent of the third subsys- 
tem, i.c., the weight percent of those products that have 
been completely mixed and that are in a state of chemi- 
cal equilibrium, be considered the quantity characteriz- 
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ing the thermodynamic component of the nonideality of 
detonation. This quantity may henceforth be tentatively 
termed the degree of the products’ reaction. The results 
of calculations of the detonation of an emulsion explo- 
sive and glass microsphere mixture presented in Figure 
1 and Table 2 correspond to complete equilibrium in 
the products or a 100-percent degree of reaction. 


It follows from the data presented in Table 2 that not 
Only the detonation rate but also the other detonation 
parameters calculated within the framework of the vari- 
ous models differ markedly from one another. The very 
low detonation temperatures calculated in accordance 
with th BKW, BKW-RDX, BKW-TNT, and BKW- 
RDX"* models are . The detonation tempera- 
tures in the BKW-R and BKW-RR models correspond 
to the freezing temperatures of the detonation products 
of powerful explosives during the process of isentropic 
expansion.‘ 


In the first model, the calculated and experimental 
values of the detonation raie of an emulsion explosive 
(without glass microspheres) virtually coincide, which 
makes it possible to assume that the detonation of the 
emulsion matrix is ideal and to calculate nonideality 
from the said value. In the calculations based on the 
BKW-RR model, it was not the absolute value of 
the detonation rate that was used as the criterion of 
thermodynamic nonideality of detonation but rather the 
decrease in detonation rate that occurred as the amount 
of glass microspheres in the mixture was increased and 
the mixture’s density was thereby decreased. 


Figures 2 through 4 present the results of calculations of 
thermodynamically nonideal detonation of an emulsion 
explosive and glass microsphere mixture by using the 
BKW-RDX model and assuming a density of 1,000 kg/ 
m’. The very low temperature of the detonation products 
is noteworthy. Along with the othur detonation parame- 
ters, Figure 4 presents the dependence of the amount of 
chemical heat released in the products at the Chapman- 
Jouget point on the degree of reaction. The amount 
of chemical heat was calculated as the difference be- 
tween the values of the standard enthalpy of formation 
of the starting system and the detonation products at the 
Chapman-Jouget point given a temperature of 298.15 K 
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The sharp bend in the curve is the result of release of 
the heat of the phase transition of SiO, as the tempera- 
ture of the products changes. The sharp bend does not 
affect the relationships of the other detonation parame- 
ters, however. 








$ 
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Figure 3. Calculated dependence of the change in 

detonation temperature of an emuision explosive 

and glass microsphere mixture on the degree of 

product reaction given an initial density of 1,000 
kg/m’. 
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depending on the structure of the free carbon formed in 
a0 pune. The Ease capesmemateas (made by the 
least squares method) of the dependence of the degree 
of reaction of the products at the Chapman-Jouget point 
On initial density are expressed in terms of the following 


For BKW-RDX: 

diamond leg, x = 0.7554p, — 0.8; 
graphite leg, x = 0.1046p, — 39.9; 
For BKW-RR: 

diamond leg, x = 0.1651, — 123.5; 


graphite leg, x = 0.1239, — 75.9; 


where x is the degree of reaction expressed as a 
percentage and p, is the initial density of the emulsion 
Sunes afd ts cdieonasiame ett is tae. 
It may thus be concluded that the thermodynamic 
nonideality of the detonation of an emulsion explosive 
and glass microsphere mixture increases in a more 
or less linear manner as the density of the mixture 
decreases and is characterized by a rate of change in 
degree of product reaction at the Chapman-Jouget point 


on the order of 0.1 to 0.15 percent/(kg/m’). 








nonideal detonation of an emulsion ex- 
mixture and the resul- 
tant values obtained for the degree of reaction. For the 
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sake of comparison, the same figures also contain the 
results of calculations performed for the case of com- 
pletely ideal detonation. Except for temperature (Figure 
10) and chemical heat (Figure 11), the detonation pa- 
rameters as a function of density of the starting mixture 
behave similarly for all of the models considered. 
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Figure 9. Dependence of the detonation rate of an 
emulsion explosive and giass microsphere mixture 
on initial density. Calculated lines: BKW-RDX, 
100-percent reaction; 2, BKW-RR, 100-percent 
reaction; 3, BKW-RDX, partial reaction; 4, 
BKW-RR, partial reaction. Points represent the 
extrapolation of experimental data for an infinitely 
large charge diameter.' (Y-axis represents 








detonation rate in m/s.) 
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Figure 10. Dependence of the detonation 
temperature of an emulsion explosive and gizas 
mixture on initial density. Calculated 
lines: BKW-RDX, 100-percent reaction; 2, 
BKW-RR, 100-percent reaction; 3, BKW-RDX, 
partial reaction; 4, BKW-RR, partial reaction. 
(X-axis represents initial density in kg/m’.) 
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In the calculations based on the BKW-RDX model 
the temperature of ideal detonation remains virtually 
unchanged and amounts to approximately 1,000 K, 
whereas using the BKW-RR model results in a decrease 
from a temperature of 1,800 K in the case of ideal 
detonation of the emulsion matrix to 1,100 K in the case 
of the maximum degree of nonideality. The amounts of 
chemical heat calculated for ideal detonation obtained 
by using the two models are virtually identical, whereas 
the values calculated for thermodynamically nonideal 
detonation differ significantly. This is because of the 
difference in the products’ temperature behavior: Given 
the low temperature of the BKW-RDX model, silicon 
oxide undergoes two phase transitions, which in turn 
affects the dependence of chemical heat on initial 
density. 


«ha /lue 
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Figure 11. Dependence of the chemical heat 
released in the products of the detonation of an 


emulsion explosive and giass mixture 
on initial density. Calculated lines: 1, BK W-RDX, 
100-percent reaction; 2, BKW-RR, 100-percent 
reaction; 3, BKW-RDX, partial reaction; 4, 
BKW-RR partial reaction. (X-axis represents initial 
density in kg °. 


The assumption that the reaction of products of the det- 
onation of an emulsion explosive and glass microsphere 
mixture at the Chapman-Jouget is incomplete, possibly 
because of insufficient diffusion of the products of the 
decomposition of the emulsion matrix's components or 
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possibility because of the insufficient rate of occurrence 
of the chemical reactions on account of low tempera- 
tures, thus made it possible to achieve calculated depen- 
dences of the detonation rate on initial density that were 
comparable to the dependences obtained in experiments. 
The order of the degree of thermodynamic nonideality, 
which is a numerical characterization of the degree of 
reaction of the products, may be estimated by a linear 
dependence of the degree of reaction on the initial den- 
sity with a slope coefficient on the order of 0.1 to 0.15 
percent x (kg/m’)'. 


According to calculations performed by using the BKW- 
RDX model, the diffusion factor should be the main 
cause for the increase in nonideality (decrease in de- 
gree of reaction) as the concentration of hollow micro- 
spheres in the mixture with the emulsion matrix is in- 
creased. The reason for this is that the temperature of 
the products of thermodynamically nonideal detonation 
remains virtually unchanged (see Figure 10). Using thc 
BKW-RR model makes it possible to consider both the 
diffusion factor and the effect of chemical kinetics as 
result of the significant temperature drop of thermody 
namically nonideal detonation. 
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Institute, Russian Academy of Sciences; first paragraph 
is KHIMICHESKAYA FIZIKA abstract; manuscript 
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(FBIS Translated Text] A study examined the lasing 
mechanism of a pulse laser at the transition between 
components of the fine structure of a chlorine atom, i.e., 
*P.,7P,, (A = 11.33 ym) during photodissociation of 
iodine monochloride [IC1]. The dependence of output 
power on the following was studied: energy and spectral 
composition of the pumping, pressure of the ICI and 
its chemical purity, and the number of pumping pulses 
before the mixture is replaced. The important role of 
processes of quenching of the Cl° on the photolysis 
products, i.c., the I atoms and especially the Cl° atoms 
that impede chemical purification of the lower laser 
level under conditions of high-power lamp pumping 
because of its preferential drain in reactions with the 
native ICi molecules, was noted. A quasi-continuous 
lasing model, including one with solar pumping, was 
shown to be possible in principle, and the conditions 
required for achieving it were discussed. 


Introduction 


In the past few years, the following two types of pho- 
todissociative gas lasers have aroused special interest: 


1. Lasers with solar pumping (space or ground based) 
capable of converting solar energy directly into the 
energy of laser radiation. The success of such lasers is 
contingent upon having an active medium that is capabie 
of operating with a high coefficient of converting solar 
energy in a closed cycle in a continuous or periodic 
pulse mode. 

2. Lasers operating in the mid-infrared range and lasing 
in an atmospheric window (A = 11-13 ym) for purposes 
of transmitting energy to satellites. probing deep space, 
and conducting space communications or else in a range 
(A = 2-5 ym) where interference on the part of the 
earth's thermal radiation with laser communication is 
minimal provided that the lasers are located outside the 
atmosphere. 


Three active media for gas lasers with solar pumping 
have been proposed and studied in the literature: the fa- 
miliar iodine atom laser, which has a wavelength of A 
= 1.315 jm and uses photodissociation of various per- 
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fluoroalky! iodides'*; the atomic bromine laser, which 
operates at a wavelength of A = 2.714 um with photol- 
ysis of IBr’*; and the molecular CO, laser, which oper- 
ates at A = 10.6 xm with photoinduction of Br,/CO, by 
E-V transfer of energy from an excited bromine atom 
to carbon dioxide.** None of these active media fully 
satisfies the aforementioned criteria. The search for new 
laser media for gas lasers with solar pumping is there- 
fore both important and timely. 


From that standpoint, this study examines the character- 
istics of a photodissociative ICI laser. 


lodine monochloride was first used as an active medium 
for a photodissociative laser with high-power lamp 
pumping in our laboratory.’* The wide and intensive 
band of continuous absorption of ICI between 360 and 
600 nm (¢ = 111 I/mol/cm at 470 nm), which is the 
region of maximum solar radiation at the level of the 
earth, and the high quantum yield of chlorine atoms 
in the excited state Cl(?P,,.) during photodissociation of 
ICI (7. ~ 0.55 with wideband pumping*) make the said 
laser attractive for use with solar pumping. 
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Figure 1. Potential energy curves and absorption 
spectrum of ICI. Key: 1, E, 10°/cm; 2, «, Vmol/cm. 
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Figure 1 presents the potential energy curves and 
absorption spectrum of ICI. In the first absorption 
band, the probability of transition to the state B’I1,+ is 
significantly higher than the probability of transition to 
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the state A°Il, and 0°.’ ICI is therefore initially excited 
into the state B°Il,+ 


ICI (X",) + Av © IC] (B°M,»), 


360 << 600 nm (la) 


which may dissociate adiabatically 


IC) (B'M,+) © 1 + CLCPy), 


hs 548 nm (1b) 


yielding a chlorine atom in an excited state, or nonadi- 
abatically 


1Ct (B'N,+) © IC (B0*) #1 + CI @P,). (1c) 


resulting in unexcited atoms. The domination of channel 
(1b) explains the inversion of the populations at the tran- 
sition between levels of the fine structure Cl?P,,-*P,,.) 
(A = 11.33 ym) in the very act of photodissociation. 


The rapid secondary reaction of the unexcited chlorine 
(Cl) atoms with the native IC] molecules*"'."? 


Ci (?P,,) + ICi Ch, +], 
k, = 1:10" c/s (2) 


occurs at a rate that is more than 30 times faster than 
the total rate at which Cl° is deactivated (here and 
henceforth X° # X(*P,,), X = Cl, Br, 1)" 


Cl WP,) + ICl~ Cl, +], i 2a) 


Ch (*P,,) + ICI, 


2 33-10-%o ig OP! 
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which additionally increases and maintains the inversion 
population at the laser transition thanks to selective 
depopulation of the lower laser level.” The prerequisites 
for a laser with a continuous operating mode thus exist. 


The radiative life time Cl(}P,,) is very long (87.9 
seconds'®); however, the cross section of stimulated 
radiation is only half that in the analogous transition 
of the iodine atom used in the familiar laser with A = 
1.315 ym.'* At room temperature (K = 3.2 x 10°),'” the 
equilibrium constant 


IC] 2 Wal, + Val, 


facilitates the regeneration of ICI molecules after pho- 
tolysis and thus permits multiple lasing based on one 
mixture. 


To clarify the mechanism of lasing, a study was 
conducted to determine the dependence of the laser's 
Output power on the energy and spectral composition 
of light pumping, the pressure and chemical purity of 
the ICI, and the number of pumping pulses before the 
mixture was replaced. The effect of the addition of 
buffer gases to the working mixture was also examined. 
In addition, the method of laser atomic resonance 
spectroscopy was used to investigate the time profiles 
of inversion in the active medium. 


1. Experimental Unit 


Figure 2 is a diagram of the experimental unit. A quartz 
laser cell with a length of 120 cm and inner diameter of 
2.8 cm and with NaCl Brewster windows had a 0.5-cm- 
thick outer “sleeve” that could be filled with gaseous or 
liquid materials when necessary to filter the pumping 
light. A quartz pulse lamp with xenon filling (pressure, 
Xe ImF 15 torr; length, 120 cm) that was, together with 
the laser tube, located inside an aluminum foil reflector 
served as the light pumping source. The energy of the 
discharge loop ranged from SO to 1,350 J, and the light 
pulse had a duration of approximately 6.5 js at its half- 
height. An optical cavity with a base of approximately 
2 m was formed by mirrors with a dielectric coating: 
an “opaque” mirror with a radius of curvature of 2 m 
and a flat mirror with a reflection factor of 95 or 90 
percent at 11.33 ym. The energy of the laser radiation 
was measured with an IMO-2N calorimeter. The laser 
pulse’s time characteristics were recorded by using a 
Ge:Au photoresistor (that had been cooled to 77 K) with 
a set of calibrated attenuators. An FEK-19 photoelement 
was used to control the lamp’s light. 
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Figure 2. Diagram of the experimental unit: 1, 
sounding laser; 2, study laser (amplifier); 3 and 3’, 
mirrors of the multipass circuit; 4, FEK-19 
phototubes; 5, Ge:Au photoreceivers; and 6, 
oscillographs with memory. 


All of the experiments were conducted under static 
conditions. The pressure measurements were taken with 
reference vacuum gauges and a type MAS-EZ pressure 
transducer with electrical output (upper limit, 30 torr) 
that was connected directly to a measuring cell. The 
iodine chloride was prepared as follows: Cl, with a 
grade of chemically pure was dried over P,O, three 
times and subjected to vacuum distillation at 77 K. It 
was then mixed in a stoichiometric ratio (1:1) with L.. 
After the reaction had been completed, the composition 
of the gaseous phase over the ICi was analyzed and 
periodically checked based on its absorption spectra’® 
on a Specord UV-VIS spectrophotometer. The typical 
composition of the gas had a purity of approximately 
99 percent. 


2. Experiment Results 
2.1. Characteristics of a Cl Laser 


Even with pumping energies significantly above the 
threshold level and long pumping pulse durations (up 
to 15-20 ys), the lasing had the form of a pulse with 
Ti, ~ 5 us and a delay relative to the beginning of 
light initiation of 1-2 js. In view of the aforementioned 
effect of “eating away” of the lower laser level in the 
active medium of a chiorine laser (provided that the 
lower value of the constant &, is accurate), such short 
pulses were somewhat unexpected because in 
level laser with collision resettlement of the 

level, the duration of the lasing pulse is, as is 
known, dictated by the duration of the pumping 
. In a study‘ that examined the characteristics of 
laser operating with photodissociation of IBr and 


se 


i 
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having a kinetic diagram identical to that of a Cl laser, 
lengthening of the pumping pulse did not lead to a 
Significant increase in the duration of lasing either. 
This indicates that in neither laser's kinetic diagram 
was consideration given to additional processes that 
either depopulate the upper laser level or slow down 
resettlement of the lower laser levels. 


The pumping threshold amounted to 100 +/- 20 J, 
which is just 1.5 times higher than that of a bromine 
laser in the equivalent cavity.“* This corresponds to 
approximately 3,000 solar constants and, in principle, 
indicates the possibility of obtaining Cl lasing with solar 
pumping by using solar energy concentrators, which as 
is well known have a theoretical limit of approximately 
20,000 solar constants.‘ As the pressure of the IC! rises, 
the threshold lasing energy begins to increase when 
(ICl} 2 7 torr (Figure 3). Consequently, operation under 
elevated pressures will only be possible given elevated 
degrees of solar energy concentration. 
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Figure 3. Dependence of the lasing energy of a Cl 


laser on pumping energies given different 
concentrations of ICI: 1.9 (1), 3.8 (2), 7.5 (3), and 


25.3 (4) mm Hg. The x-axis represents pumping 
energy in joules, and the y-axis represents the 
lasing energy in millijoules. 


The dependence of lasing energy on energy stored in 
the discharge loop (proportional to the light pumping 
energy) presented in Figure 3 is linear throughout the 
entire range of pumping energies ‘EF. = 0.1 to 1.35 kJ) 
and ICI pressures (1 to 25 torr) studied. For a bromine 
laser, the increase in output power virtually ceased after 
E_ & $00 kJ,“*"* which the authors of another study‘ 
attributed to either heating of the active medium or 
quenching of the excited Br° atoms on the iodine atoms. 
The absence of lasing energy saturation in an IC! laser 
may accordingly be explained as follows: a) the smaller 
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cross section of absorption of IC] compared with IBr 
and, consequently, the lesser amounts of heating; b) the 
stronger temperature of the rate constants 
of the interaction of CI° and Cl atoms with ICI molecules 
as compared with the analogous dependences of the rate 
constants of the reaction of Br° and B atoms with [Br 
molecules; and c) the stronger quenching of Br’ atoms 
on the photolysis products of IBr (I, Br, and Br° atoms 
and oscillation-excited molecules IBr(1))) as compared 
with the g of Cl’ atoms on the photolysis 
products of ICi (1, Cl, and Cl atoms and oscillation- 
excited molecules ICI(v)): 


Cr+1+C+i, (4) 
T+Cr-Aacrn+c, 
r+a-a+a, (6) 
+10 (9-0, +1 
“A+. ( 


Unfortunately, the literature data regarding spin-orbital 
relaxation of Br° and Cl° atoms on atoms and radicals 
are extremely scant.*'*?' The temperature dependences 
of the rate constants of the interaction of bromine atoms 
in ground and excited states with parent IBr molecules 
that were obtained by Boriyev et al.” confirm that in 
an IBr laser, the effect of temperature is extremely 
important: k/k* is approximately equal to 35 at 300 K 
versus only approximately 10 at 350 K. The temperature 
of the rate constants of the reactions of Cl 
and Cl° atoms with ICI are lacking in the literature. 


The effect of inert gas additives on a Cl laser’s output 
characteristics was studied to clarify the role of heating 
of the active medium. It turned out that the addition 
of He, Ar, and Xe does not result in an increase in 
the intensity and/or duration of the laser's lasing pulse. 
Moreover, a monotonic, albeit rather weak decrease in 


transition line CiC?P,, *P,,) in the said atoms’*). As is 
thus evident in Figure 4, for 9 torr, IC], and E_ = 726 
J, adding 30 torr of He and Ar, as well as 40 torr of 
Xe, reduced E __. by approximately half. Oxygen had a 
slightly stronger effect on lasing: When approximately 
15 torr of O, was added, the laser's lasing energy (E_._] 
dropped by half. 
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Figure 4. Dependence of the lasing energy of a Cl 
laser on concentration of the different additives to 
the active medium [M]: M = Ar (1), He (2), Xe (3), 
and O, (4); [ICI] = 3 x 10"/cm’, E = 726 J. The 
x-axis represents the of a given 
additive, i.e., [M){ICI], and the y-axis represents 
the resultant lasing energy divided by the initial 


lasing energy, ie. E/E...’ 


Besides demonstrating that heating has no effect on the 
mechanism of lasing, the said results make it possible 
to estimate the upper bound of the rate constant of the 
quenching of Ci(?P,,) atoms on O, and noble gases. 
The estimates are consistent with data published in two 
studies'?*; however, they contradict the results of other 
published studies. In one of those studies,” a rate 
constant of 2.1 +/- 0.5 x 10°'' cm*/s was obtained by the 
method of vacuum ultraviolet adsorption spectroscopy 
for the process of quenching of Cl" on O,. That value 
was confirmed in a later study by the same team of 
researchers.» Clearly, with such a rate constant, adding 
even approximately | torr of O, with a lasing time 
scale of 10 ys should lead to a significant decrease in 
lasing energy. (Such behavior actually occurred when 
strong relaxants were used as additives: hydrogen and 
a number of multiatomic gases, ic., CF,Cl, CCl,, 
CF, Br, CF,I, i-C,F,1, CH,I, CD,I, and SF,, completely 
quenched the lasing energy even when added in very 
small concentrations (<1 torr), in total consistency with 
the values of the deactivation rate constants.*) Our 
estimate of the rate constant of the quenching of CI° 
on O, was ks 3x 10" cm’/s. 
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A strong increase in rale constant in the series from 
He to Xe (from 3.8 x 10° cm*/s in the case of He to 
1.8 x 10°" cm/s in the case of Xe) was discovered in 
one study” for the process of relaxation of an excited 
chlorine atom in collisions with atoms of inert gases. 
In another study,”’ the asymptotic method of plotting 
quasi-molecular terms was used to calculate the rate 
constants of spin-orbital relaxation, i., FCP,,) and 
Ci@P,,.), during collisions with He, Ne, Ar, Kr, and 
Xe, and an inverse relationship was obtained: for the 
case of quenching of Cl", the rate constant decreases 
monotonically from 4.2 x 10 to 1.4 x 10” cm/s 
in the series from He to Xe. As is evident from the 
curves presented in Figure 4, Ar and Xe are, at the 
very least, not more effective quenching partners for 
Ci?P,,) than He is. From our measurements and in 
view of experimental data on the quenching of excited 
chlorine atoms on Ar atoms'*#' and the results of a 
theoretical study,”’ it is obvious that the values presented 
elsewhere” for the rate constants of the relaxation of Cl° 
on atoms of heavy noble gases are inflated by several 
orders of magnitude. 


It is difficult to determine the reasons for so great a 
discrepancy between the data obtained in some pub- 
lished studies** and our studies*’ and other published 
data."*°2* This is especially true when consideration 
is given to the fact that the rate constants presented in 
the one group of studies*” for relaxation of CI° on O, 
and noble gases are highly inflated whereas the rate con- 
stants presented in the same studies** for the cases of 
collision of Ci?P,,) atoms with molecules of H,,?!?°2* 
CO,’ CH, HCI” Cl,,"? and H,O” are low. Suffice 
it to present the opinion of the one team of researchers”® 
that the time profile of the concentration, i.e., (CI°(s)], 
with lamp photolysis of CCI, (the source of CI° atoms in 
the aforementioned group of studies) is rather com- 
plicated and its decrease is not described in terms of a 
simple exponent as was proposed in the studies of Hu- 
sein’s team.*** 


_. on the concen- 
tration of ICI, ic., EL. = E({ICl}), that are presented 
in Figure 5 are distinguished by the presence of sat- 
uration at an iodine chloride pressure above 10 torr. 
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or the aforementioned processes (described in equations 
(3) through (7)) of spin-orbital relaxation Ci?P_,.) of the 
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state on ICI and its photolysis products: iodine, chlorine 
atoms and oscillation-excited molecules ICI(v). 


‘, whe 
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Figure 5. Dependence of the lasing energy of a Cl 
laser on the concentration of IC] at the 
different pumping energies [E___: 600 (1), 864 (2), 
1,037 (3), and 1,350 (4) J. The x-axis represents the 
concentration of IC] in mm Hg, and the y-axis 
represents the lasing energy in millijoules. 


The maximum energy of lasing obtained in a cavity 
with a 10% output mirror amounted to 4.5 +/- 0.5 
mJ with E ~ 1 kJ and an iodine chloride pressure 
of approximately 18 torr. This is approximately one- 
third the lasing energy of a Br laser in an equivalent 
Cavity in the optimum mode, which allowing for the 
difference in energy of the radiated quantum of hy, / 
hv, = 4, indicates that both lasers have similar quantum 
efficiencies. 


The absorption spectrum of ICI consists of two bands of 
identical intensity that are located in the visible (band 
1) and ultraviolet (band 2) regions of the spectrum.'° 
Experiments involving filtration of a pulse lamp's 
radiation were performed to determine the contribution 
of the ultraviolet absorption band of IC] (A__ = 248 
am, o,., = 4.4 x 10°'* cm’)"* to the formation of CiCZP,, 
,) excited atoms. For that purpose, gaseous Br, under a 
pressure of approximately 200 torr was let into the laser 
cell's “sleeve,” which in turn weakened the luminous 
flux in the region A > 360 nm by a factor of 7-10 
virtually without touching the ultraviolet portion of the 
lamp’s radiation. No laser radiation was observed when 
E__ = 1.3 kJ and when the pressure of the iodine 
chloride amounted to about 9 torr. It is thus possible to 
draw the qualitative conclusion that during photolysis 
of ICI in the ultraviolet region of the spectrum with 
A < 360 nm, the main atoms formed are chlorine 
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atoms in the ground state *P,, and +,,. < 0.33. It is 
thus the photodissociation of ICI in the visible band 
of absorption of light with 360 < A < 548 nm that is 
decisive to the formation of Cl° atoms. Consequently, 
with solar pumping of ICI, it is necessary either to cut 
off the ultraviolet portion of the radiation or, better yet, 
transform it into the visible region. Using a solution of 
the dye 490 RT (1.8-naphthylene-1',2’-benzimidazole) 
in benzene not only enabled us to eliminate photolysis 
of ICI in the ultraviolet region of the spectrum but also, 
thanks to the dye’s luminescence, enabled us to increase 
the density of the luminous flux in the visible region 
of the spectrum and significantly improve the laser's 
Operating stability. 
Photodissociative Br lasers are known to be sensitive 
to molecular bromine impurities in IBr on account of 
the process of parasitic photodissociation of the Br, that 
the lower laser level.*'* An analogous 
effect may also be important for a Cl laser. Studies 
were therefore conducted to determine the effect that 
additions of Cl, 
its lasing . It turned out that a chlorine laser is 
much more in the face of such impurities, and 
its lasing is completely cut off only when [C1,):(ICl] = 
1. Using the aforementioned dye 490 RT in the laser 
“sleeve” improves the stability of the laser's output 
characteristics when — to Cl, impurities by an 
additional factor of 1.5-2. 


In cases of repeated photolysis (up to 50 flashes) 
without replacement of the active substance and with 
an interval of 1 minute between pulses, the intensity 
and shape of the laser signal and its delay time remained 
virtually unchanged from pulse to pulse. This confirms 
the effective regeneration of the starting IC] molecules 
in recombination processes and indicates that an ICI 
laser may operate in a sealed version. 


2.2. Time Profile of an Inversion 


Studies of the time behavior of inversion of the popula- 
tions at the Ci?P, .-*P,) transition in the active medium 
of a Cl laser wi different pumping levels are very 
important to assess the possibility of creating a quasi- 
continuous Ci laser with solar pumping. To take such 
measurements, the experimental unit was modified in 
accordance with a “generator-amplifier” circuit." An 
analogous photodissociative laser model was added to it 
(Figure 2) that served as a source of sounding laser ra- 
diation at the transition between levels of the fine struc- 
ture of the chlorine atoms, and the study Cl laser was 
transformed into an amplifier (its resonator mirrors were 
removed and replaced with mirrors permitting multi- 
ple'* passage of the sounding radiation through the ac- 
tive medium). The sounding laser was triggered with 
a delay relative to the pulse of the amplifier lamp [7,,, 
anging from 0 to 1,000 ys. The sounding laser's lasing 


to a Cl laser's active medium have on 
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signal was a single pulse with a duration of 1.5 jus. The 
gain [x] of the laser’s signals as it passed through the 
cell four times was found by using the relationship x 
= In(LJ,°/1,1,°, where I, and I, are respectively the laser 
signal's measured intensities at the cell's entrance and 
exit and /,° and /,° are the analogous intensities for an 
empty cell. 


The time dependence of the gain factor of a C1 laser, i.e., 
x1), was measured under the following conditions: 


a) with an identical pumping energy, 'c., FE = | 
kJ (the degree of photodissociation of ICI was kept 
constant at a = 3 percent) but with the IC! under 


different levels of pressure (0.5, 1.0, and 3.0 torr); 


b) with an identical ICI pressure (3 torr) and different 

ing levels, ic., E = 1 kJ (a = 3 percent) and 
E_ = 270J (a = 0. cent). Figure 6 presents the 

ts. It turned out that the decrease in the curves 
plotted for ys) with one 2:id the same pumping level 
but different ICI pressures is described by an identical 
overall constant [k,“"] of 5 x 10? cm’/s. At the same 
time, the drop in the curves given identical values 
of [ICl} but different values of a is described by 
Significantly different values of k,*", i.e., (k,*"), = 5 x 
102 cms (E__ = 1 kJ, a = 3 percent) and (k,*), = 
8 x 10 cms TE = 270 J, a = 0.9 percent). The 
dependence of the value of k, on the degree of the 
active substance’s collapse may be clearly traced, with 
the k, decreasing and approximating the published” 
value as a decreases. It is significant that channels 
(IV) through (VII) of the quenching of Cl(?P,,.) on 
the photolysis products of ICI play an important role 
even with degrees of collapse as low as a 2 | percent 
and must be taken into consideration in the kinetic 
diagram describing lamp photolysis of ICI with high- 
power pumping. 
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Figure 6. Dependence of intensification of a C] 
laser’s resonance signal during photolysis of IC1: 
(ICl] = 1 x 10°/em’; E__ = 1,000 (1) and 270 (2) J. 
The x-axis represents time in microseconds, and the 

y-axis indicates ,“ 
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All of the photolysis products of ICI, i.e., atoms of 
Cl, Cl", and I, are easily observed in our method 
when the active medium is sounded with | and Cl 
lasers. Of course, the formation of oscillation-excited 
ICK(X, v) molecules is in principle possible when IC! is 
excited with white light. As is clear from the absorption 
spectrum of ICI that is presented in Figure 1, however, 
the relative concentration of such “hot” molecules barely 
exceeds 3-5 percent of the total number of dissociated 
molecules of ICI. It is therefore not likely that process 
(VII) plays an important role in the kinetic diagram 
of photolysis of ICI. (Therein lies one of the main 
differences between IC] and IBr. In the case of wideband 
photolysis of molecules of IBr, the quantum yield of the 
photodissociation may amount to approximately 0.25- 
0.5, which means that the concentrations of atoms and 
"hot" molecules, i.e., IBr(X, v) are comparable.”) 


We will note that the concentration of I atoms increases 
on account of process (II), reaching a stationary state al 
the level 2a[ICi}, whereas the concentration of chlorine 
atoms in the ground state [C1] falls symbasically to zero 
(k, = 1.5 x 10" cm/s). Using a Ci laser to sound a 
measuring cell containing IC] under a pressure of 3 torr 
mixed with CCl, under a pressure of | torr, which is an 
effective quencher of the state Cl(?R,.), demonstrated 
that there are virtually no Cl atoms left after about 5 
us following the beginning of the light pulse. At our 
request, the authors of one study*' measured the value of 
the constant k, by the method of electron paramagnetic 
resonance in a flow and obtained a value of 4.7 +/- 
1.5 x 10°? cm’/s. It may therefore be hypothesized that 
the contribution of process (VI) to quenching of the Cl° 
atoms is, in all likelihood, also small. 


In addition, given the very different kinetics at short 
times, the rear “tails” of the amplification profiles x1) 
at times ¢ 2 12 ys virtually coincide for degrees of 
collapse of IC] of 3 and 0.9 percent (see Figure 6) 
despite the fact that the concentration of iodine atoms 
under the said conditions differs by more than a factor 
of 3. This indicates that quenching of the Cl@P,,,) atoms 
on I atoms (process (I'V)) is less important than is their 
relaxation in process (V). The profile of x(t) calculated 
by the Runge-Kutta method from a kinetic diagram that 
includes processes (1) through (VI) is consistent with 
the results of an experiment performed with values of 
k, = 5 +/- 2x 10" and k, = 10 +/- 5 x 10°" cm’/s, 
respectively. We will note that while the value obtained 
for k, by variation falls within the range of published 
values for the rate constants of the process of relaxation 
of X(?P.) atoms on I atoms,'*° the value obtained for k, 
appears to be anomalously high at first glance. However, 
the literature already mentions just such an anomaly for 
the quenching of I?P,.) on Cl atoms.”-” 
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Final clarification of the contribution of processes dic- 
tating the kinetics of inversion during photolysis of ICI 
requires experiments with independent sources of Cl’, 
Cl, and I atoms. The results that have been obtained 
nevertheless demonstrate that during high-power lamp 
photolysis of IC] (when a > | percent), additional chan- 
nels of quenching of Ci@P,,.) on photolysis products 
appear and sharply increase the overall constant of its 
decomposition so that k,*“/k, comes to equal to approx- 
imately 1/2. Using reaction (II) to selectively purify the 
lower laser level and thereby create a quasi-continuous 
Cl laser, including a Cl laser with solar pumping, be- 
comes problematic. Evidently, a continuous lasing mode 
can only be achieved in the case of "weak" pumping 
with degrees of dissociation [a] of markedly less than 
1 percent, wherein quenching of the Cl° atoms on the 
photolysis products is not yet significant. 


The existing discrepancy between the values of the rate 
constant of the process of deactivation of Cl° molecules 
by ICI molecules, i.e., k,, measured in published stud- 
ies'*'> may be explained by proceeding from the results 
obtained. In one of the two studies,'® high-power lamp 
photolysis of IC] (E. = 1.2 kJ, a = 5 percent) was 
used, and the measured value of k, = 3 x 10" cm/s 
was evidently due to the effect of secondary processes 
of quenching on the photolysis products. In the other 
study,'? laser photolysis of IC] (a < 0.5 percent) was 
used, the likelihood of a contribution on the part of the 
secondary processes was low, and the value k, = 3.3 
x 10'? cm’/s obtained was obviously closer to reality. 
The reason why the value obtained for k, in a recent 
publication’ was inflated by nearly an order of mag- 
nitude is more difficult to determine inasmuch as the 
authors also used laser photolysis. It can only be spec- 
ulated that either the degree of photodissociation of ICI 
in the said study'* was higher or else the accumulation 
procedures used (the said procedures were used" to im- 
prove the signal:noise ratio) made it impossible to avoid 
the effect of secondary processes of quenching on the 
photolysis products. 


3. Conclusion 


The output characteristics of a pulse photodissociative 
Cl laser and the time behavior of the inversion in its 
active medium have thus been studied. The fundamental 
possibility of implementing a quasi-continuous lasing 
mode, including one with solar pumping, has been 
demonstrated. An analysis of the reasons for the rapid 
decomposition of the inversion indicated that processes 
of quenching of Cl° atoms on the photolysis products, 
i.c., atoms of I and especially Cl, play an important 
role in the said laser's lasing mechanism. Chemical 
purification of the lower laser level by its preferential 
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therefore only be implemented with low levels of light _ KVANTOVAYA ELEKTRONIKA, Vol 12, No 8, 
pumping. 1985, p 1914. 
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Footnote 


“It should be noted that the spread between the values 
provided in the literature for the total rate constant of 
processes (3a) and (3b) is very large. Donovan and Hu- 
sain’? give a value that is two orders of magnitude larger 
than given elsewhere,'* and a recent publication" 
presented an intermediate value, ic., k, = (4 +/- 1) x 
10°? cm’/s. One of the purposes of the present study 
was therefore to clarify the reasons for so great a dis- 


crepancy. 
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Using the Electrical Response of a Semiconductor 
Compound To Record the Chemical 
Transformation Behind the Front of a Shock Wave 


957A1206B Moscow KHIMICHESKAYA FIZIKA 
in Russian Dec 94 Vol 13 No 12, pp 175-176 


[Article by S.S. Nabatov and A.V. Lebedev, Chem- 
ical Physics Institute, Russian Academy of Sciences, 
Chernogolovka; manuscript received 12 Jul 94; UDC 
$39.63+537.226] 


(FBIS Translated Text} A method of recording exother- 
mal reactions in shock wave conservation experiments 
based on the phenomenon of thermoelectromotive force 
in semiconductors has been proposed.' In essence, the 
method entails measuring the difference in potentials 
between the faces of a semiconductor element that is 
placed into a flat storage tube together with the study 
sample. it was noted that such a formulation can only 
be used to study processes in an unloaded state. The 
method was tested by way of the example of the reac- 
tion occurring during the synthesis of stannous mono- 
sulfide (SnS} from the elements. A thermoelement was 
made from SnS, which is a p-type semiconductor with a 
thermoelectromotive force coefficient of +550 V7 un- 
der normal conditions. It turned out that with a loading 
pressure of 16 GPa, the reaction is initiated within sev- 
eral milliseconds after the shock wave eniers the study 
specimen. There are, however, grounds for hypothesiz- 
ing that at higher pressures, the synthesis may occur di- 
rectly in a shock-compressed state. As has been shown 
elsewhere,’ an SnS thermoelement generates electrical 
signals with a characteristic amplitude of 200-300 mV 
when a shock wave passes through it. This commu- 
nication proposes that the said signals be used to ob- 
serve the reaction occurring in the system Sn + S during 
shock compression. The experiment consisted of sound- 
ing a shock-compressed mixture with a secondary shock 
wave. It was hypothesized that if a semiconductor com- 
pound forms behind the front of the first shock wave, 


the two loading amplitudes were insufficient 
reaction in a shock-compressed 
state, which is consistent with findings reported else- 
where.' Additional experiments with higher pressures 
planned to clarify whether it is possible to effect a 
synthesis reaction behind the front of a shock wave. 
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Estimating the Parameters of a New Class of 
Detonation of Emulsion Explosives 


9S7AI1199A Moscow KHIMICHESKAYA FIZIKA 
in Russian Dec 94 Vol 13 No 12, pp 131-144 
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Russian Academy of Sciences, Moscow; first paragraph 
is KHIMICHESKAYA FIZIKA abstract; manuscript re- 
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(FBIS Translated Text] This article is devoted to quan- 
titative estimation of the parameters characterizing the 
detonation of mixtures of an emulsion explosive and gun- 
powder based on thermodynamic calculations of ideal 
and nonideal detonation. Two models of the equation of 
State of gaseous detonation products, i.e., the Becker- 
Kistyakovskiy-Wilson [BKW]-RDX and BKW-RR mod- 
els, were used for the calculations. Three thermody- 
namic models of detonation were examined: 1) ideal 
detonation, i.e., complete chemical equilibrium of the 
detonation products at the Chapman-Jouget point; 2) 
Slightly less than ideal detonation, i.e., the detonation 
products consist of two chemical subsystems (decom- 
position products of the emulsion matrix and gunpow- 
der) such that a state of equilibrium exists within each 
subsystem and the two systems do not interact with one 
another; and 3) far-from-ideal detonation, i.e., the det- 
onation products consist of three chemical subsystems 
(decomposition products of the aqueous ammonium ni- 
trate solution, hydrocarbon mixture, and gunpowder) 
such that a state of equilibrium exists within each sub- 
system and the three systems do not interact with one 
another. It is demonstrated that estimates of detonation 
parameters based on thermodynamic calculations and 
general discussions of the mechanism of emulsion explo- 
sives’ sensitivity make it possible to boost the said explo- 
sives’ power characteristics by 10-20 percent or more by 
replacing the hollow glass microspheres or gas-forming 
additives with disperse explosives, especially gunpow- 
der, the problem of whose recycling in view of conver- 
sion of the defense complex still remains far from solved. 


It has been hypothesized elsewhere' that the power 
characteristics of emulsion explosives that are now 
in existence and used in practice may be increased 
significantly by adding gunpowder to them. This study 
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is devoted to quantitative estimation of the detonation 
parameters of mixtures of emulsion explosives and 
gunpowder. 

Emulsion explosives possess good water resistance and 
are fairly explosion-proof, which is to say that they have 
a low sensitivity to initiation of detonation. While low 
sensitivity is a positive feature during manufacturing 
and transport, it becomes a drawback when the explo- 
sives are put to direct use inasmuch as the traditional 
method of boosting their sensitivity by adding hollow 
microspheres or a gas bubble-forming chemical reagent 
to the emulsion matrix results in a decrease in charge 
density and, accordingly, in a decrease in the efficiency 
of the detonation products of emulsion explosives. 


Adding disperse gunpowder to the emulsion matrix in- 
stead of microspheres or gas-forming chemical additives 
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appears to be a promising way of overcoming the afore- 
said drawback in using emulsion explosives. When that 
is done, a smaller reduction or possibly even an increase 
in detonation parameters as the amount of additive is in- 
creased may be anticipated instead of the rather sharp 
drop in detonation parameters that occurs when hollow 
glass microspheres are added to the emulsion matrix.’ 


Table 1 presents the composition of an emulsion matrix 
identical to the one considered elsewhere.? The compo- 
sition of the gunpowder whose detonation parameters 
were studied in another study’ (gunpowder No. 3) and 
that was considered as an additive to commercial explo- 
Sives in yet another study’ are also presented in Table 
1 along with the enthalpies of the various components. 























Table 1. 
Component Chemical Formula Weight Percent Eathalpy, kj/ag _ 
Ammonium nitrate H,N,0, 77 -4,563 
Water H,O 16 15,860 
Saturated hydrocarbon (liquid) (CH,), 6 -2,088 
Oleic acid C,gH,,0, 1 -2,841 
Guapowder CoHy3 13% 15 Pec more than 100 percent -2,272 














The emulsion matrix’s density was determined to equal 
1,353 kg/m’.* For the sake of generalization, two 
gunpowder density values, i.e., 1,600 and 1,353 kg/m’, 
have been cons'dered. The latter value was selected to 
clarify the effect of adding gunpowder apart from the 
effect of the initial charge density. 


Two different cases of determining charge density were 
thus considered: 1) an increase in charge density from 
1,353 kg/m’ by adding gunpowder with a density of 
1,600 kg/m’ and 2) a charge density of 1,353 kg/m’ that 
remains constant regardless of the amount of gunpowder 
in the mixture. Figure | 1s a graphic representation of the 
dependence of charge density on the amount of additive 
for each of the two aforesaid cases. The dependence of 
the density of the mixture of an emulsion matrix and 
hollow glass microspheres is presented for comparison’s 
sake. In Figure 2 the relative increase in initial charge 
density when gunpowder is used as compared with when 
hollow glass microspheres are used to boost emulsion 


explosives’ sensitivity is presented as a function of the 
amount of additive. 


The Automated System of Thermodynamic Calculations 
and Algorithms [ASTRAL] was used to perform the 
thermodynamic calculations of emulsion explosives’ 
detonation based on a method described elsewhere.’ A 
Becker-Kistyakovskiy-Wilson [BKW] equation served 
as the equation of state of the detonation products’ 
gaseous phase. The possibility of formation of disperse 
diamond and graphite was assumed for the condensed 
carbon.° A Kovena [transliteration] equation of state was 
used for the graphite, and the diamond was considered 


incompressible. 


Two sets of coefficients were used for the BKW equa- 
tion: BKW-RDX and BKW-RR. Selection of the two 
sets of coefficients for calculating emulsion explosives’ 
detonation is discussed elsewhere.’ 
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Table 2. 
Exploatve Density, kg/m? Model Detonation Rate, m/s Pressare, GPa Temperature, K 
Guapowder 1,353 BKW-RR 6,706 15.68 3,375 
1,353 BKW-RDX 6,568 15.44 2,834 
Emulsion mani 1,353 BKW-RR 6,613 13.48 1,781 
1,353 BKW-RDX 7,134 16.04 976 









































Figure 3 is a graph representing the calculated equilib- 
rium dependence of the detonation rate on the amount 
of gunpowder in a mixture with an emulsion matrix. 
The results of calculations based on the BKW-RDX and 
BKW-RR models differ not only from the standpoint of 
absolute values but also with respect to the nature of 
the dependences. In fact, when the initial charge den- 
sity remains unchanged (1,353 kg/m’), the detonation 
rate calculated on the basis of the BKW-RDX equa- 
tion moderates as the amount of wder in the mix- 
ture is increased, whereas the BKW-RR equation yields 
a slight increase in detonation rate. This is a conse- 
quence of the fact that the calculated detonation rates 
a. ped ny are slightly different for the BKW-RDX 
~ % BKW-RR 8, whereas the detonation rate of 
Ae matinee RDX mode! is signifi iy Wher taan leet calociened 

is significantly hi calculated 

<o are placed by using the BKW-RR equation (Table 2). Moreover, 








a AS 4 A 








of the emulsion matrix detonation rate calculated on the 

5; 2, gumpowder, basis of the BKW-RDX model is noticeably higher than 

: that calculated for gunpowder with an identical initial 
charge density. 
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Figure 5S. Dependence of the temperature of the 
products of an ideal detonation on the weight 
percent of gunpowder in a mixture with an 
emulsion matrix. 1, BKW-RDX model and 
gunpowder density of 1,600 kg/m’; 2, BKW-RDX 
model and gunpowder density of 1,353 kg/m’; 3, 
BKW-RR model and gunpowder density of 1,600 
kg/m’; 4, BKW-RR model and gunpowder density 
1,353 kg/m’. Values on the y-axis represent 
quantities in K. 


The advantage of using gunpowder instead of hollow 
microspheres in a mixture with an emulsion matrix 
may be graphically illustrated by the ratio of the 
detonation rates of the said mixtures given an identical 
weight percent of additive in the mixture. Figure 6 
is a graph representing the dependence of the ratio 
of the detonation rates of mixtures of emulsion with 
gunpowder and with microspheres for the two equation 
of state models and two methods oi determining the 
charge density of a mixture with gunpowder. Previously 
calculated’ values obtained by using the BKW-RDX 
and BKW-RR models and assuming partial reaction 
of the detonation products by matching the slopes of 
the calculated and experimeatally obtained dependences 
of the detonation rate on initial charge density were 
used as values for the detonation raie of mixtures of 
an emulsion matrix and microspheres. As follows from 
the graph in Figure 6, when the amount of additive 
to the emulsion matrix is increased, the ratio of the 
detonation rates increases sharply. The calculated rate of 
the detonation of an emulsion explosive with gunpowder 
is thus approximately twice that of the detonation rate 


of a mixture with microspheres when the weight percent 
of additive is assumed to be approximately 10 perceat. 
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Figure 6. Dependence of the ratio of the rate of ideal detonation of a mixture with gunpowder to rate of 


ideal detonation of a mixture with 


on the weight percent of additive. Linus: 1, BKW-RDX 


microspheres 
model and gunpowder density of 1,600 kg/m’; 2, BKW-RDX model and gunpowder density of 1,353 kg/m’; 
ee a Oe 
1 B, 


Those results were, however, obtained in the case where 
completely ideal detonation was assumed for the mix- 
tures with gunpowder and thermodynamically nonideal 
detonation with incompletely reacted products at the 
Chapman-Jouget point was assumed for the mixtures 
with microspheres. For that reason, to examine the effect 
of the possibility of thermodyramic nonideality of the 
detonation of mixtures of an explosive emulsion with 
gunpowder, we will consider two additional cases. The 
first is the worst possible case, and it may serve as the 
lower bound of our estimate of the said mixtures’ deto- 
nation parameters. 


The first case of thermodynamic nonideality of detona- 
tion may be modeled by assuming that the detonation 
products of the emulsion matrix and gunpowder do not 
engage in a chemical reaction with one another. When 
such is the case, the entire system is assumed to con- 
sist of the following three chemical subsystems whose 
products are in a state of chemical equilibrium: 1) a so- 
bution of ammonium nitrate in water; 2) a mixture of 


The calculation results demonstrate that the absence of 
chemical interaction between the decomposition prod- 
ucts of the emulsion matrix and gunpowder have virtu- 


ally no effect on the parameters of the mixture’s deto- 
nation. This is expressed in the indistinyuishability of 
the lines depicting the dependence of the detonation pa- 
rameters of mixtures of emulsion matrix and gunpowder 
in Figures 3 through 5 for completely ideal and ther- 
modynamically nonideal (two subsystems) detonation. 
The dependence of the ratio of the detonation rate of a 
mixture of emulsion «nd gunpowder to the detonation 
rate of a mixture of emulsion and microspheres on the 
weight percent of additive in the case of thermodynam- 
ically nonideal detonation virtually ¢«.»cides with the 
corresponding dependence for completely ideal detona- 
tion. 


A different pattern is observed in the case of thermody- 
namic nonideality and the formation of three chemical 
subsystems in the detonation products. The said case 
corresponds to the cace of thermodynamically nonideal 
detonation with a zero degree of product reaction that 
was examined elsewhere.’ For brevity’s sake, we will 
refer to the system in question (three subsystem) as ther- 
modvnamically nonideal detonation with a zero degree 
of reaction. Because the decomposition products of the 
oxidizing subsystem (solution of ammonium nitrate in 
water) do not react with the decomposition products of 
the burning subsystem (hydrocarbon mixture), a sub- 
stantial portion of the chemical component will not be 
transformed ito the thermal and elastic components of 
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the internal energy to the Chapman-Jouget point. As a 
result, the detonation parameters turn out to be signifi- 
cantly lower than in the cases considered previously. 


Figure 7 is a graph represenung the dependences of 


the rate of thermodynamically nonideal detonation of 


a mixture of emulsion matrix and gunpowder on the 
amount of additive as calculated by using the two 
equation of state models and two methods of specifying 
initial charge density in the case of a zero degree of 
product reaction. The dependences differ from the cases 
of ideal and slightly less than ideal detonation. The 
detonation rate values for a zero degree of reaction 
that are obtained with both models are significantly 
lower than that for a pure emulsion. As the amount of 
gunpowder in the mixture increases, the detonation rate 
increases and approaches the detonation rate of ideal 
and slightly less than ide. detonation. 


Ti.¢ possibility of obtaining higher detonation parameter 
values when gunpowder is added to the emulsion matrix 


instead of hollow microspheres, even in the case of 


detonation with a zero degree of product reaction, 1s 
illustrated in Figure 8, which graphs the dependences 
of the ratio of the detonation rates of mixtures with 
gunpowder and a zero degree of product reaction to the 
detonation rates of mixtures with microspheres on the 
amount of additive in the emulsion mixture. 
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Figure 7. Dependence of the rate of far.’rom-ideal 
(with a zero degree of product reaction) detonation 
on the weight percent of gunpowder in a mixture 
with an emulsion matrix. Lines: 1, BKW-RDX 
model and gunpowder density of 1,600 kg/m’; 2, 
BKW-RDX model and gunpowder density of 1,353 
kg/m’; 3, BKW-RR model and gunpowder density 
of 1,600 kg/m’; 4, BKW-RR model and gunpowder 
density 1,353 kg/m’. Values on the y-axis represent 
quantities in m/s. 
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Figure 8. Dependence of the ratio of the rate of a far-from-ideal detonation (with a zero degree of product 
reaction) of a mixture with gunpowder to the rate of detonation of a mixture with microspheres, on the 
weight percent of additive. Lines: 1, BKW-RDX model and gunpowder density of 1,600 kg/m’; 2, 
BKW-RDX model and gunpowder density of 1,353 kg/m’; 3, BKW-RR model and gunpowder density of 
1,600 kg/m’; 4, BKW-RR model and gunpowder density 1,353 kg/m’. 
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The detonation rate of pure emulsion in the case of 
a zero degree of reaction is less than the rate of 
ideal detonation, which is what a comparison with 
experimental data’ indicates that detonation of pure 
emulsion is, and amounts to more than 80 percent of the 
rate of ideal detonation, When, however, the amount of 
additive in the mixture is increased to levels as low 
as 1.5-2.5 percent (by weight), the rate of far-from- 
ideal detonation of a mixture with gunpowder becomes 
comparable to the rate of somewhat nonideal detonation 
of a mixture with microspheres. 


When the amount of additive is increased further, the 
ratio of the rates increases rapidly. Al a concentration of 
additive equal to 9 percent (by weight), the detonation 
rate of a mixture with gunpowder and a zero degree of 
product reaction exceeds the rate of somewhat nonideal 
detonation of a mixture with microspheres by 50-60 


percent. 


it follows from experimental data‘ that the critical di- 
ameter of detonation of a mixture of emulsion and mi- 
crospheres, which characterizes the charge’s sensitiv- 
ity, reaches 25 mm and virtually stops decreasing at 
densities below 1,200 kg/m’, which corresponds to ap- 
proximately a concentration of microspheres in the mix- 
ture equal to 1.6 weight percent. The experimentally 
obtained* dependences of detonation rate on the initial 
density peak at approximately the same value. It may 
be assumed that the values of the rate of somewhat 
nonideal detonation of a mixture of emulsion and mi- 
crospheres calculated on the basis of the BKW-RDX 
and BKW-RR models with the given charge density are 
optimum for actual use from the standpoint of the re- 
lationship of the detonation rate and critical diameter. 
The specified detonation rate values calculated in ac- 
cordance with the two equation of state models may 
therefore be used as reference values for determining 
the effectiveness of adding gunpowder instead of mi- 
crospheres. Specifically, detonation rates that have been 
calculated for mixtures of emulsion and gunpowder by 
using the respective models for the two methods of de- 
termining charge density may be referenced to the said 
values. 


The graphs presented in Figures 9 and 10 represent the 
ratios of the detonation rates of mixtures of emulsion 
matrix with gunpowder to reference values of the deto- 
nation rate of a mixture of emulsion matrix and micro- 
spheres as a function of the amount of gunpowder in 
the mixture for the two different methods of determin- 
ing mixture density. 
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Figure 9. Dependence of the relative rat: of 
detonation on a mixture’s weight percent of 
gunpowder with a density of 1,353 kg/m’. Lines: 1, 
ideal detonation, BKW-RDX model; 2, ideal 
detonation, BKW-RR model; 3, far-from-ideal 
detonation, BKW-RDX model; 4, far-from-ideal 
detonation, BKW-RR model. 


It is evident from Figure 9 that when the density of the 
mixture equals the density of the emulsion matrix, it 
may be expected that adding gunpowder will increase 
the detonation rate by 10-15 percent over the reference 
values. The anticipated detonation rate of a mixture of 
emulsion and gunpowder should lie in the regron whose 
upper bound is the dependence for an ideal detonation 
rate and whose lower bound is the the dependence for 
detonation with a zero degree of product reacuon. We 
will note the different natures of the dependences for 
ideal detonation rate calculated in accordance with the 
different equation of state models. In the case of the 
BKW-RR model, the relative rate of ideal detonation 
increases slightly as the amount of gunpowder in the 
mixture is increased, whereas it decreases in the case 
of the BKW-RDX model. The values obtained for the 
relative rate of detonation with a zero degree of reacuion 
increase with both models, but to different degrees. In 
any case, the rate of detonation with a zero degree 
of product reaction surpasses the reference detonation 
rate of a mixture of emulsion and microspheres given 
a concentration (by weight) of gunpowder in a mixture 
with an emulsion matrix of more than 12 percent. 
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Figure 10. Dependence of the relative rate of 
detonation on a mixture'’s weight percent of 

gunpowder with a density of 1,600 kg/m’. Lines: 1, 
ideal detonation, BKW-RDX model; 2, ideal 

detonation, BKW-RR model; 3, far-from-ideal 
detonation, BAW-RD\ model; 4, far-from-ideal 
detonation, BKW-RR model. 


In the case in which the density of the mixture increases 
as its content (by weight) of gunpowder with a density 
of 1,600 kg/m’ ts increased (Figure 10), the increase in 
the rate of ideal detonation when gunpowder is added 
reaches values 15-20 percent above the reference values 
for a mixture with microspheres. The nature of the de 
pendences of the relative rates of ideal and maximally 
possible nonideal detonation ts identical for both equa- 
tion of state models: In all cases, the detonation rate 
increases as the concentration (by weight) of gunpow- 
der in the mixture is increased. The only difference ts 
in the dependences’ slopes. The detonation rates, even 
with a zero degree of product reaction, are more than 10 
percent higher than the reterence values for the given 
concentration (by weight) of gunpowder in the mixture 


It has been hypothesized** that detonation of an emul- 
sion matrix is ideal from a thermodynamic standpoint, 
e., the detonation products are completely intermixed 
and react by the Chapman-Jougct point. Increasing the 
amount of hollow microspheres increases the mixture ’s 
sensitivity but results in an increase in the thermody 
namic nonideality of detonation, as modeled> by an in- 
crease in the degree of chemical reacuon of the de 
composition products of the aqueous ammonium nitrate 
soluuon and hydrocarbon mixture 


Under the effect of the leading shock wave, the hollow 
glass microspheres collapse, thereby creating “hot spots 
in whose vicinity chemical reactions in the emulsion 
matrix itself are initiated and accelerated in the inal 
Stages. The creation of significant local temperature 
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spikes, Leo, hot spots, is, however, accompanied Dy a 
Jecrease in the total amount of energy per anit mass 

volume 
displaced trom part 


replaced by inert glass 


Of mixture Decause enerpy-releasing matter is 
of the mixture’s bulk and partially 
which if turn results ina 


Gecrease in the degree of completion of the chemica 


reacuions im the late stages of the reachon process 
itler dispersion oof the hot spots. As a result. the 
increase in temperature of the products in the chenuc al 


reacuon zone as they move to the Chapman-Joupgel 
point decreases significantly as a result of the lower 
erage Concentration of energy in the system. whi 
In turn restrains the increase in the detonation products 
temperature. The result is a closed cirele, or in other 
words, a system with positive feedback: The decrease 


in the rate of energy release results in a decrease 
temperature rise that in turn leads to a decrease in rate 
of energy release that in turn leads to a decrease 
temperature rise, and so on 


Increasing an emulsion explosive’s sensiuvity by the 
iddition of hollow glass microspheres thus leads to a 
decrease in the degree of reacuon of the detonation 
products 


As has been noted,’ achieving detonauon conditions 
In gunpowder requires powerful initiation and a rather 
large charge diameter. Systematic experimental data 
regarding the sensitivity of gunpowder to shock wave 
Inithauion and its critical diameter are not to be tound 
in the literature, however. Consequently, we can only 
make a few general remarks that must be verified by 
systematic experimental studies 


With regard to the possibility of the development ot 
a detonation mode of propagation of the chemical 
decomposition of gunpowder under the effect of a shock 
wave passing from the emulsion matrix, we can only 
state that the rate and pressure of detonation of a pure 
emulsion are comparable with the detonation rate and 
pressure of gunpowder (see Table 2) and that it may 
therefore be anticipated that the strength of a shock 
wave 1s entirely adequate to initiate detonation in a 
gunpowder particle 


The critical diameter of detonation, even if it exceeds 
the dimensions of a gunpowder particle, is not directly 
applicable to the case under consideration inasmuch as 
the particle is surrounded by a dease medium subjected 
to shock wave compression 


The proposed mechanism of increasing the sensitivity of 
an emulsion explosive by adding gunpowder particles 1s 
based on the very same principle as in the case ¢ 
microspheres or gas bubbles. The temperature of the 
products of the gunpowder’s detonation 1s significantly 
higher (by 1,S00-2.000 K for different equation of state 
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models) than the temperature of the detonation products 
of pure emulsion, In that sense, the gunpowder parucles 
will act as “hot spots.” The products of the gunpowder’s 
detonauon act as “hot spots” not only in close-lying 
regions, as in the case of microspheres, but also in 
significantly more remote regions inasmuch as al least 
as much energy is stored in gunpowder molecules as 1s 
stored in the explosive emulsion itself. 


Even with a very cautious estimate, it may therefore be 
anticipated that the degree of reaction of the detonation 
products of a mixture of emulsion matrix and gunpow- 
der particles will be at least equal to that occurring upon 
detonauon of 4 mixture of emulsion and microspheres. 
A more realistic estimate is the hypothesis of the ther- 
modynamic ideality of the detonation of a mixture with 
gunpowder or nonideality caused by partial mixing or 
total nonmixing of the detonation products of the emul- 
sion and gunpowder as a result of the dimensions of the 
latter's microparticles. In the latter case, the aforemen- 
tioned scenario involving the formation of two chem- 
ical subsystems in the detonation products may serve 
as a mode! of nonideality of detonation. The results of 
calculations performed for the said scenario practically 
coincide with the results of calculations performed for 
ideal detonation. 


Estimates of detonation parameters based on thermody- 
namic calculations and general discussions of the mech- 
anism of the emulsion explosives’ sensitivity thus make 
it possible to conclude that the strength characteristics 
of such explosives can realistically be increased by 10- 
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20 percent or more (depending on the reference values) 
by replacing the sensitivity-increasing hollow glass mi- 
crospheres and gas-forming additives with disperse ex. 
plosives, especially gunpowder, the problem of whose 
recycling in view of conversion of the defense complex 
sull remains far from solved 
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New Hypothesis on the Role and Importance of the 
Electric Field in Oceans 


W4D0070A Khabarovsk TIKHOOKEANSKAYA 
GEOLOGIYA in Russian Vol 14 No 1, May-Jun 95 
(manuscript received 10 May 94) pp 121-123 


[Annotauon and Article by E. A. Alftan, commentary 
by V_ A. Akhmadulin, St. Petersburg; UDC $50.37] 


{[FBIS Translated Text} 


Annotation 


This arucle examines the reasons for the generation 
of a telluric electric current in oceans. It is shown 
that it flows not only in the water, but also to a 
significant extent below the ocean floor and continents. 
This explains the formation of uplifts and troughs on 
the surface of the oceans above those areas of the 
ocean floor which have a high electric conductivity. 
Even small shifts in the Earth's crust or the covering of 
underground channels of juvenile waters may lead to a 
drastic change in the vertical component of this current, 
which in a number of cases leads to the formation of 
tsunamis and earthquakes, the generation of infrasonic 
sound and, as a result of the latter, to catastrophes 
involving ships and aircraft. 


The materials of the article also confirm and supplement 
the little known hypothesis of V. S. Borkhsentus on the 
oceanic Origin of the earth's magnetuc field. 


{Text of article] 


In many places, the surface of the ocean has smooth 
deflections of tens of meters, that 1s, unique upifts and 
troughs. Usually they are explained by a hypothesis 
which presumes the presence of gravitational anoma- 
lies. These deflections of the ocean can be explained 
by the effect on the water of electrostatic fields in the 
depths of the ocean and below it. The existence of such 
fields 1s evidenced by the presence of electric telluric 
currents in oceans. Water, like other dielectrics, is at- 
tracted to a zone of greater electrostatic field strength. 
Consequently, if there are local regions in the ocean in 
which the vertical component of the electrostatic field 
reaches a significant value, the surface of the water 
should deflect downward or upward in these zones. Pos- 
sible orders of magnitude of the electrostatic field are 
indicated in the calculations of V. §. Borkhsenius.’ He 
showed that when water evaporates and condenses the 
configuration of molecules changes, as does their elec- 
tric charge. Consequently, between heated and cooled 
zones of the ocean there occurs a difference of potentials 
and an electric current flows. Since dunng evaporation 
and condensation each molecule changes in charge by 
| eV, the total difference of potentials that is formed is 
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sufficient to explain not only the ongin of tellune cur 
rents, but according to the hypothesis of V. S. Borkhse- 
nius, the Earth's magnet field as well 


The meridional direction of the tellurie current in the 
Atlanuc confirms the presence of differences in potential 
between heated equatorial and cooled arctic regions ol 
the oceans. A sample calculation using data presented 
in Rei. 4 shows that on average a current on the order 
of 1,000,000 A flows in a band of ocean 100 km wide 


According to the hypothesis of V. S. Borkasenius the 
Earth's magnetic field is due to an electric current 
flowing in the oceans in 4 direction of parallels between 
the day side of the Earth heated by the sun and the 
cool night side. Actually, on planets which do not have 
oceans one observes only a weak magnetic field. The 
essential objection to Borkhsenius’ hypothesis ts that 
several continents lie in the path of the oceanic current 
in the equatorial direction. The path of the current 
around them in the ocean would significantly weaken 
it. Borkhsenius died without finishing and popularizing 
his hypothesis. 


To solve our problem we seek the vertical component of 
the electrostatic field while at the same time refining and 
supplementing Borkhsenius’ hypothesis. Actually, 1s 11 
possible that the electric current which arises due to the 
differences in potential at the surface of the ocean flow 
only along the water? Obviously, it should branch and 
flow in the depths of the Earth under the ocean floor 
and under continents, where the electric resistance of 
the rock is significantly less than the electric resistance 
of ocean water. This is shown schematically in Fig. 
1. Current | flows through the ocean from zone A of 
the heated ocean to zone B of the cooled ocean. The 
electric conductivity of ocean water is about 2-5 O 
'm'. In the Gutenberg layer’ (under the substrate) the 
temperature rises to 1500-1800°C, that 1s, the melung 
point of minerals. Their electric conductivity reaches 
1000 O'm'. In deeper layers of the mantle and tn the 
Earth's core the electric conductivity rises by several 
orders of magnitude. The basalt rocks of the ocean floor 
are a good insulator. Thus, according to the notation in 
Fig. 1, the electric resistance of the lower mantle Rn 
and the core Rec 1s significantly less than the electric 
resistance of the ocean Ro, but between them are layers 
of greater electric resistance, Rks, Rs, Rg. However, 
in those p'aces where this layer is disrupted or very 
thin, a significant part of the electric current 1s directed 
from the ocean downward or from beneath the ocean 
floor upward. Consequently, a current flows along the 
lower mantle (current I,) and the Earth's core (current 
I.) from point A to point B parallel to current I in the 
ocean. Currents I, and I, obviously also flow beneath 
continents between the day and night sides of the Earth, 
participating in the formation of the Earth's magnetic 
field. 
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Figure 1. Schematic of electric resistances of a portion of the ocean and layers of the Earth beneath it. 


Current | moves through the ocean from zone A of the heated ocean to zone B of the cooled part of the ocean 
The electric resistance of the lower part of the mantle Rm and core Re are significantly less than the resistance 
of the ocean Ro, but between them are layers of large resistance Rks, Rs, and Rg. The calculations show that in 
zones where the integrity of the earth's crust is disrupted (the emergence of juvenile water, magma, etc.) the electric 
resistance Rkv is significantly less than Rk, and thus a substantial portion of current | is directed from the ocean 
into the depths and flows parallel to the current in the ocean from point A to point B along the lower mantle (current 
L,) and the Earth's core (current I,). Key: a. ocean; b. Earth’s crust; c. lithosphere d. upper layer (substrate) of the 


upper mantle; e. Golitsin layer. 


In the Bermuda Triangle there are several zones of 
high electric conductivity. First, there is the Puerto 
Rican depression. Above it the surface of the ocean 1s 
depressed downward by 15-25 m. The juvenile waters 
have a high electric conductivity and their suboceanic 
bed links the ocean with the depths of the Earth. 
These waters probably emerge in the Sargasso Sea (the 
depression of the surface is about 60 m). An excellent 
conductor of electric current is the lava of underwater 
volcanos, for example, at the Middle Atlantic fault. 


Since, as noted above, a current on the order of 
1,000,000 A flows through every 100 km band of ocean, 
and in regions where the integrity of the ocean floor 
is disrupted the current may be directed from zones 
with a width of several hundred kilometers, the value of 
the vertical component in the areas of these disruptions 
also may be on the order of millions of amperes. This 
type of current obviously, is associated with significant 
electrostatic fields. 


Since the electric current in these zones 1s constant, or 
changes very slowly, deflections of the ocean surface 
retain their size and physical properties. The value of 
the electric current flowing under the ocean or from 
beneath the floor of the ocean may change drastically, 
for example, as a result of small shifts in the layers 
of rock under the ocean floor in a region of high 
eleciric conductivity. This may lead to the generation 
Of infrasonic sound, to the formation of long waves 
(tsunamis) and earthquakes 


An electric current in the stream of juventie waters helps 
to dissolve rock, eroding the underground channe!. The 
separated chunks of rock may partially or completely 
cover its underground bed for some time. Partial co 
erage may lead to the formation of an infrasonie liquid 
or vapor gush below the ocean. If the frequency of the 
infrasonic sound coincides with the frequency of the al 
pha rhythm of the brain, then this may lead to insanity 
and death of people in the affected zone. The infra- 
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sonic sound may froth the water at a depth of tens or 
hundreds of meters. The frothed water is lighter than 
common water and a ship in this region will collapse 
under the water, 


Moving from the ocean into the atmosphere and creating 
an alternaung decrease and increase in air density, 
the infrasonic sound leads to an alternating (with the 
frequency of the infrasonic sound) decrease and increase 
in lift, affecting the (wing) planes of aircraft. This, 
obviously, 1» che reason for the known, but unexplained 
phenomenon which occurs when even in calm weather 
an aircraft over the ocean suddenly begins to shake, 
which sometimes leads to a catastrophe. The drastic 
change in the vertical component of the electric current 
due to the movement of rock or covering of an 


underground part of the beds of juvenile waters in itself 


may generate an infrasonic sound pulse, because the 
rock of the ocean floor has piezoproperties. This may 
lead to the aforementioned phenomena, as well as to 
earthquakes. 


If the current-conducting channel below the ocean floor 
is covered for a long time, the deflecuon upward 
of the ocean surface begins to overflow, generating 
long waves (tsunamis). The nature of these tsunamis 


was unexplained because they arise in the absence of 


cyclones and other known causes. 


Using the invention’ of the author of this article, use of 


the laws described above makes it possible to determine 
the points of origin of these tsunamis and the moment 
of their formation This makes it possible to promptly 
warn the population of the coastal zones of the ocean 
and save tens of thousands of lives. 
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Discussion by V. A. Akhmadulin 


The proposed hypothesis is based on the litte known 
ideas of V. S. Borkhsenius on the orgin of electric 
Charges in the evaporation and condensation of water, 
which themselves require proof. The author of the article 
should have presented in more detail the essence of this 
hypothesis and presented quantitative calculations, as 
well as shown how such fantastic figures (1,000,000 
amperes per 100 km or 10 A/m) may be obtained 
from the initial data in reality. In addition, it should be 
pointed out that the latitudinal currents of Borkhsentus, 
which flow trom the heated day side to the cool 
night side, interact with meridional currents which 
arise due to differences in temperature between heated 
equatorial and cooled arctic regions of the ocean 
In addition it was important to esumate the currents 
which are induced in conducting sea water by the 
Earth's external electromagnetic field. It 1s completely 
unclear what conducting branches in the depths of 
the Earth the author ts talking about, because the 
conductivity of sea water is much higher than the 
conductivity of sedimentary rock, not to mention the 
magmatic and metamorphic rock comprising the upper 
part of the Earth's crust of oceans and continents 
Practically, current conducting zones, in the opinion of 
the author, should penetrate to the lower mantle, but 
the arguments in support of this view are not presented, 
although the experience of geoelectricians in the study 
of conducting objects in the Earth's crust and upper 
mantle exclude the presence of such vertical channels 
There is no basis to link this with the lava of underwater 
volcanos of the Middle Atlantic fault. The effect of 
these channels would undoubtedly have been recorded 
by geoelectricians. It is also unclear what stream of 
juvenile water the author is talking about, implying that 
under the ocean floor there may exist underground rivers 
of the type which are formed in karst cavities. All of 
this requires proof because tt contradicts available data 
and existing ideas 





This report may contain copynghted matenal. Copying and dissemination 
18s prohibited without permission of the copynght owners 

















FBIS-UST-95-051 
18 December 1995 


Shares Sold in Seven New Oil Companies 


964D0121A Moscow KOMMERSANT DAILY 
in Russian 2 Nov 95 pS 


[Arucle by Olga Gorokhova under the “Oil Company 
Stock Auctions” rubric: “Seven New Oil Holdings 
Sold” | 


[FBIS Translated Text} The lethargy of the financial 
market has finally been shaken by an event that 
investors could not even have dreamt of half a year 
ago. The sale of shares in seven new oil companies 
began straight off at the all-Russian financial auction 
conducted by the Federal Fund Corporation (FFK) 
on behalf of the Russian Federal Property Fund 
({[RFFD. Because of a number of factors, these 
valuable pieces of paper may turn out to be an 
extremely profitable acquisition. 


At the beginning of August of this _ year, 
Goskomimushchestvo [Russian Federation State 
Committee on the Administration of State Property] 
issued a directive confirming the plan for selling 
shares designated state property. The list contained 
136 enterprises, including the following newly formed 
holdings: Eastern Oil Company (VNK), Eastern Siberia 
Oil and Gas Company (Vostsibneftegaz), Siberia-Ural 
Petrochemical and Gas Chemistry Company (SIBUR), 
Orenburg Onl Joint-Stock Company (ONAKQO), and 
Siberian-Far Eastern Oil Company (SIDANKQ). The 
RFFI later added the Sibneft Oil Company, NORSI-Oy! 
Qil Company, ard Tyumen Oil Company (TNK) to 
the list. The said enterprises’ shares have not been 
circulated on the secondary market, which makes 
potential investors’ task much more difficult inasmuch 
as investors generally use market prices as a guide 
when submitting orders. 


At the same time, rather large blocks of shares—from 
10 to 1S percent of each holding’s authorized capital— 
have been put forward. As a result, shares worth more 
than 241 billion rubles [R] (a figure based on the shares’ 
initial prices) have entered the market at the same ume. 
Clearly, so large a “surge” of shares is inconsistent 
with the drop in demand that is currently evident in 
the market. 


It should also be borne in mind that the auction results 
will not be tallied until after the elections, the outcome 
of which 1s unclear as of today. The political uncertainty 
is not very likely to inspire participation in trading even 
by Russian investors who are accustomed to everything, 
and nervous foreigners are most likely to be scared away 
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altogether. All of this would seem to indicate that the oil 
shares offered for sale may turn out to be underpriced. 
However, that is precisely what makes them a good way 
of investing capital. 


True, investors should consider the fact that each of 
the companies (except for Vostsibneftegaz) has retained 
ownership of 34 percent of its shares (in addition to the 
government-owned block of shares) and that, in accor- 
dance with the privatization plan, no more shares will 
be offered through financial auctions but will instead 
only be able to be acquired at investment competitions 
planned for 1996. And that is generally more expensive 
and more troublesome than at auctions. As far as fur- 
ther privatization of the holdings is concerned, all the 
companies except for SIDANKO are still quite young, 
and the sale of those blocks that have been deemed fed- 
eral property cannot begin until at least autumn 1998 
(the state-owned block of SIDANKO may be offered 
for sale as early as spring 1997). It will therefore be 
difficult to form a rather large block of shares. 


According to an unnamed source, large Russian banks 
and investment companies are now appearing. True, 
experience has demonstrated that most orders will 
probably be placed in the last few days of the order 
period in the second 10 days of November. The “Nashe 
otechestvo"” [Our Fatherland] IPK [Advanced Training 
jInstitute], an agent of the Federal Fund Corporation, has 
proposed that rather large Russian structures take part 
in the trades. This includes both outside oil companies 
and possibly those companies that are issuing the shares, 
which will buy their own shares. Private individuals will 
hardly want to take part in such a game. 


_It is interesting that in accordance with the Chubays 


directive dated 21 July 1995, a new method of deter- 
mining the initial price of a share is being used at the 
auctions. Previously the initial price simply could not 
exceed 20 times the face value. Now, however, book 
values are used to make the calculation. It is true that 
for all practical purposes these "fresh-out-of-the-oven"™ 
holdings have no such values, and data regarding their 
subsidiary structures are used. Thus the initial price of 
shares in SIDANKO is R120,000, and that of shares 
of ONAKO is R110,000. The initial prices of shares 
in the other companies are as follows: VNK, R80,000, 
SIBUR, R70,000; Sibneft, R115,000; TNK, R85,000, 
and NORSI-Oyl, R160,000. Only the initial price of a 
share of Vostsibneftegaz has remained at R20,000. 


{Boxed item] 
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Companies Whose Shares Are Being Offered at 
Financial Auctions: 


VNK: Authorized capital of approximately R3 billion 
Owns controlling blocks of nine Companies, including 
Tomskneft and the Tomsk NKF 


Vostsibneftegaz: Authorized capital of more than R442 
million. The holding includes the Yeniseynettegaz and 
Yeniseygeotizika 


SIBUR: Authorized capital of approximately RS14 
million. Owns controlling blocks of six companies, 
including the Scientific Research and Design Instutute 
of Gas Refining (NIPI gazpererabotki) and Sorbent. 


ONAKO: Authorized capital of approximately R2 bil- 
lion. Owns controlling blocks of six Compamies, includ- 
ing Orenburgneft and Orenburggeolgiya. 
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SIDANKO: A ithorized capital of approximately RY 4 
billion. Owns controlling blocks of 26 joint-stock ¢ 
pames, including Udmurtnett and the Angara Peu 
chemical Company (Angarskaya Netleahimicheskaya 
kompantiya), 


Sibneft. Authorized capital of approximately R?.22¢ 
billion. Owns controlling blocks of the joint-stock 
companies (AOs) Noyabrsknettegaz AQ and Omskiy 
NPZ [oil retinery] AO and a number of other companies 


TNK: Authorized capital of R7YO.882 million. Owns 
controlling blocks of the Tyumennettegaz AQ, Nizh 
nevartovskneflegaz, etc 


NORSI-Oyk: Authorized capital of R1.497 billion. Its 
mihorized capital includes blocks of shares in six 
iecovery and oil ref.aing companies, including the Nor 
AQ). 
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